fit_tidal_hsieh_1987: fit_tidal_hsieh_1987

Description Usage Arguments Examples

View source: R/tidal_hsieh_1987.R

Description

This function fits a empirical response function value to the Hsieh model. For high transmissivity wells (>1e-4 m2/sec) this will likely not give reasonable storativity values.

Usage

1
2
fit_tidal_hsieh_1987(frequency = 1.9322736/86400, gain, phase,
  radius_well, radius_casing = radius_well)

Arguments

frequency

frequency in cycles per second

gain

value to fit in m/ns

phase

phase shift relative to Earth tide

radius_well

radius of the well

radius_casing

radius of the casing

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
data('hsieh_1987_fig_2_3')
storativity <- 1e-07
transmissivity <- 1e-05
radius_well <- 0.05
frequency <- 1.9322736/86400
hsieh <- tidal_hsieh_1987(frequency, storativity, transmissivity, radius_well)

res <- fit_tidal_hsieh_1987(frequency = 1.9322736/86400,
                 Mod(hsieh$response),      # m / strain
                 Arg(hsieh$response),      # degrees
                 radius_well,
                 radius_casing = radius_well)
all.equal(res$par, log10(c(storativity, transmissivity)), tolerance = 0.05)

jkennel/aquifer documentation built on Dec. 25, 2019, 9:25 a.m.