tidal_hsieh_1987: tidal_hsieh_1987 Solution for estimating transmissivity and...

In jkennel/hydrorecipes: Hydrogeology steps

tidal_hsieh_1987 Solution for estimating transmissivity and storativity from earth tides.

Description

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])

jkennel/hydrorecipes documentation built on Feb. 23, 2025, 11:03 a.m.