slug_cbp | R Documentation |
This function calculates H/H0 for a well of finite diameter using equation 7 from Cooper, Bredehoeft and Papdopulos, 1967 Response of a Finite-Diameter Well to an Instantaneous Charge Water.
This function uses Stehfest method for the Inverse Laplace Transform.
slug_cbp(time, radius_casing = NULL, radius_screen = NULL, radius = NULL, storativity = 1e-05, transmissivity = 10, head_0 = 1, n = 12L)
time |
elapsed time since slug introduction |
radius_casing |
radius of the casing where the water level is |
radius_screen |
radius of the screened interval |
radius |
radius from borehole center to calculate the response |
storativity |
aquifer storativity |
transmissivity |
aquifer transmissivity |
head_0 |
height of introduced head change |
n |
coefficient for Stehfest algorithm (should be < 20) |
vector of H/H0
t <- 5.776 * 10^(seq(-4, 3, 0.01)) s <- slug_cbp(t, radius_screen = 7.6) plot(s~t, log = 'x', type = 'l')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.