tidal_hsieh_1987 Solution for estimating transmissivity and storativity from earth tides.

Hsieh, P.A., Bredehoeft, J.D. and Farr, J.M., 1987. Determination of aquifer transmissivity from Earth tide analysis. Water resources research, 23(10), pp.1824-1832.

Usage

tidal_hsieh_1987(
  frequency,
  storativity,
  transmissivity,
  radius_well,
  radius_casing = radius_well
)

Arguments

frequency: the frequency of the response
storativity: layer storativity (L/L)
transmissivity: aquifer transmissivity (L*L/t)
radius_well: well radius
radius_casing: casing radius

Value

tidal response

Examples

data('hsieh_1987_fig_2_3')
storativity <- 1e-07
transmissivity <- 1e-03
radius_well <- 0.05
frequency <- 10^seq(-5, 2, by = 0.05)
tau   <- 1 / frequency
hsieh <- tidal_hsieh_1987(frequency, storativity, transmissivity, radius_well)
plot(Mod(response)~dimensionless_frequency, hsieh,
 type='l',
 log = 'x',
 xlim = c(1, 1000))
points(response~dimensionless_frequency, hsieh_1987_fig_2_3[variable=='gain' & S == storativity])


plot(unwrap(Arg(response)) * 180/pi~dimensionless_frequency, hsieh,
 type='l',
 log = 'x',
 xlim = c(1, 1000),
 ylim = c(0, -90))
points(response~dimensionless_frequency, hsieh_1987_fig_2_3[variable=='phase' & S == storativity])