step_aquifer_theis: step_aquifer_theis

In jkennel/hydrorecipes: Hydrogeology steps

step_aquifer_theis

Description

Radial Flow (GRF) model with flow dimension equal to 2. This method defaults to a fast FFT convolution so many rates can be included, but requires a regular time series.

Usage

step_aquifer_theis(
  .rec,
  time,
  flow_rate,
  thickness = 1,
  radius = 100,
  specific_storage = 1e-06,
  hydraulic_conductivity = 1e-04,
  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 well (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 Theis model

References

Barker, J.A., A generalized radial flow model for hydraulic tests in fractured rock. Water Resour. Res., 24 (1988), pp. 1796-1804, 10.1029/WR024i010p01796

Theis, C.V., 1935: The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage, Transactions of the American Geophysical Union, 16th Annual Meeting, Part 2, pp. 519-524.

See Also

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

Examples

dat <- data.frame(x = as.numeric(1:100),
                  y = rep(0.01, 100))
formula <- as.formula(y~x)

frec = recipe(formula = formula, data = dat) |>
  step_aquifer_theis(time = x, flow_rate = y) |>
  prep() |>
  bake()

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