fit_tidal_hsieh_1987: fit_tidal_hsieh_1987

In jkennel/aquifer: Well solutions for aquifer tests

fit_tidal_hsieh_1987

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

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

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 July 31, 2022, 11:33 p.m.