14.07.2013 Views

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

96 Processing <strong>Modflow</strong><br />

MODFLOW < Well<br />

An injection or a pumping well is defined by using the Data Editor to assign two values to a<br />

model cell:<br />

3 -1<br />

- Rechare rate <strong>of</strong> the well (Q) [L T ] and<br />

- Parameter Number [-]<br />

Negative cell values for the “Rechare rate <strong>of</strong> the well” are used to indicate pumping wells,<br />

while positive cell values indicate injection wells. The parameter number is used to assign Q as<br />

a parameter for an automatic calibration by the inverse models PEST or UCODE, see PEST <<br />

Parameter List... or UCODE < Parameter List....<br />

The injection or pumping rate <strong>of</strong> a well is constant during a given stress period and is<br />

independent <strong>of</strong> both the cell area and the head in the cell. It is implicitly assumed by MODFLOW<br />

that a well penetrates the full thickness <strong>of</strong> the cell. MODFLOW can simulate wells that<br />

penetrates more than one model layer. In this case, the injection or pumping rate for each layer<br />

has to be specified. The total injection or pumping rate for a multilayer well is equal to the sum<br />

<strong>of</strong> those from the individual layers. For confined layers, the injection or pumping rate for each<br />

layer (Q ) can be approximately calculated by dividing the total rate (Q ) in proportion to the<br />

k total<br />

layer transmissivities (McDonald and Harbaugh, 1988):<br />

Q k ' Q total @ T k<br />

ET<br />

3.6.1 MODFLOW<br />

(3.26)<br />

where T is the transmissivity <strong>of</strong> layer k and GT is the sum <strong>of</strong> the transmissivities <strong>of</strong> all layers<br />

k<br />

penetrated by the multilayer well.<br />

An other possibility to simulate a multi-layer well is to set a very large vertical hydraulic<br />

conductivity (or vertical leakance), e.g. 1 m/s, to all cells <strong>of</strong> the well. The total pumping rate is<br />

then assigned to the lowerst cell <strong>of</strong> the well. For display purposes, a very small pumping rate<br />

-10 3<br />

(say, -1×10 m /s) can be assigned to other cells <strong>of</strong> the well. In this way, the exact extraction rate<br />

from each penetrated layer will be calculated by MODFLOW implicitly and the value can be<br />

obtained by using the Water Budget Calculator. See Chapter 2 for how to calculate subregional<br />

water budgets.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!