step_aquifer_wellbore_storage: step_aquifer_patch

In jkennel/hydrorecipes: Hydrogeology steps

step_aquifer_patch

Description

Papadopulos-Cooper 1967 solution for wellbore storage.

Usage

step_aquifer_wellbore_storage(
  .rec,
  time,
  flow_rate = 1,
  radius = 0.15,
  radius_casing = 0.15,
  radius_well = 0.15,
  thickness = 1,
  specific_storage = 1e-06,
  hydraulic_conductivity = 1e-04,
  n_terms = 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)
radius	distance from center of well
radius_casing	radius of casing in the interval over which the water level declines
radius_well	effective radius of well screen or open hole
thickness	the aquifer thickness (L)
specific_storage	specific storage of aquifer (L/L)
hydraulic_conductivity	the hydraulic conductivity (L/t)
role	character - the name of the role
...	additional arguments

Value

The drawdown using the Papadopulos-Cooper model

References

Papadopulos, I.S. and H.H. Cooper, 1967. Drawdown in a well of large diameter, Water Resources Research, vol. 3, no. 1, pp. 241-244.

See Also

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

Examples

dat <- data.frame(x = 10^seq(-5, 2, length.out = 100),
                  y = rep(0.01, 100))
formula <- as.formula(y~x)

frec = recipe(formula = formula, data = dat) |>
  step_aquifer_wellbore_storage(time = x,
                                flow_rate = 10.0,
                                hydraulic_conductivity = 10.0,
                                specific_storage = 1e-4) |>
  prep() |>
  bake()

jkennel/hydrorecipes documentation built on Dec. 24, 2024, 5:38 p.m.