step_aquifer_patch

barker_herbert 1982 solution for radial patches.

Usage

step_aquifer_patch(
  .rec,
  time,
  flow_rate = 0.01,
  thickness = 1,
  radius = 200,
  radius_patch = 100,
  specific_storage_inner = 1e-06,
  specific_storage_outer = 1e-05,
  hydraulic_conductivity_inner = 1e-04,
  hydraulic_conductivity_outer = 1e-06,
  n_stehfest = 12L,
  role = "predictor",
  ...
)

Arguments

.rec

the R6 recipe object.

time

the time for evaluation (t)

flow_rate

the flow rate from the well (L^3/t)

thickness

the aquifer thickness (L)

radius

the distance to the observation location (L)

radius_patch

the radius of the cylindrical patch (L)

specific_storage_inner

specific storage of inner patch (L/L)

specific_storage_outer

specific storage of outer patch (L/L)

hydraulic_conductivity_inner

the hydraulic conductivity of the inner patch (L/t)

hydraulic_conductivity_outer

the hydraulic conductivity of the outer patch (L/t)

n_stehfest

integer number of terms to use in Stehfest method (typically < 18)

role

character - the name of the role

...

additional arguments

Value

The drawdown using the Theis model

References

Barker, J.A., and R. Herbert, 1982: Pumping tests in patchy aquifers, Ground Water, vol. 20, No. 2, pp. 150-155.

Butler, J.J., 1988: Pumping tests in nonuniform aquifers – The radially symmetric case, Journal of Hydrology, Vol. 101, pp. 15-30.

See also

Other aquifer: step_aquifer_constant_drawdown(), step_aquifer_grf(), step_aquifer_leaky(), step_aquifer_theis(), step_aquifer_theis_aniso(), step_aquifer_wellbore_storage()

Examples

dat <- data.frame(time = as.numeric(1:100))
formula <- as.formula(time~.)

frec = recipe(formula = formula, data = dat) |>
  step_aquifer_patch(time = time, flow_rate = 0.01) |>
  plate()