GetWellConfig: Get Well Completion and Pumping Rate in Model Space

Description Usage Arguments Value Author(s) Examples

View source: R/GetWellConfig.R

Description

This function determines well completions and pumping rates in model space. The pumping rate is specified for each model cell intersecting a well's open interval(s) and calculated by multiplying the estimated pumping demand by the cell's transmissivity fraction. The transmissivity fraction is calculated by dividing the cell's aquifer transmissivity by the sum of all transmissivity values for cells belonging to the same well. The transmissivity fraction calculation assumes saturated conditions in the model cell.

Usage

1
GetWellConfig(rs.model, wells, well.col, rate.col = NULL, lay2.hk.tol = 0.01)

Arguments

rs.model

RasterStack. Composed of raster layers describing the model grid and hydraulic conductivity distribution: lay1.top, lay1.bot, lay2.bot, lay3.bot, lay1.hk, lay2.hk, and lay3.hk.

wells

SpatialPointsDataFrame. Average pumping rate for each well during various times.

well.col

character. Column name of the well identifier field.

rate.col

character. Vector of column names for the pumping rate fields.

lay2.hk.tol

numeric. Hydraulic conductivity tolerance for model cells in layer 2. Used to prevent pumping in the aquitard layer of the aquifer system. Pumping is prohibited in model layer 2 cells with hydraulic conductivity values less than lay2.hk.tol and a well opening isolated to layer 2; for these cases, pumping is allocated to the adjacent layer 1 cell.

Value

Returns an object of class data.frame with the following components:

...

unique identifier assigned to a well, its name is specified by well.col.

lay,row,col

layer, row, and column number of a model cell, respectively.

hk

hydraulic conductivity of the model cell, in meters per day.

thk

vertical length of the well opening (open borehole or screen) in the model cell, in meters. A value of zero indicates that the well opening is unknown or below the modeled bedrock surface.

frac

transmissivity fraction for a model cell, where transmissivity is defined as hk multiplied by thk.

...

pumping rate allocated to the model cell for each time period specified by rate.col, in cubic meters per day. The pumping rate is calculated by multiplying the pumping demand for a well (specified in wells) by frac.

Author(s)

J.C. Fisher, U.S. Geological Survey, Idaho Water Science Center

A.H. Wylie, Idaho Department of Water Resources

Examples

1
## Not run: # see Appendix D. Uncalibrated Groundwater-Flow Model

USGS-R/wrv documentation built on June 30, 2020, 11:07 p.m.