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

radius_well

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 Dec. 24, 2024, 5:38 p.m.