Slug test steps
Cooper Bredehoeft Papadopulos 1967
# check vs. CBP 1967 table 1
Tt_r2 <- rep(c(1.0, 2.15, 4.64), 4) * 10^(rep(c(-3, -2, -1, 0), each = 3))
a_01 <- c(0.9771, 0.9658, 0.9490, 0.9238, 0.8860, 0.8283,
0.7460, 0.6289, 0.4782, 0.3117, 0.1665, 0.07415)
a_0001 <- c(0.9969, 0.9949, 0.9914, 0.9853, 0.9744, 0.9545,
0.9183, 0.8538, 0.7436, 0.5729, 0.3543, 0.1554)
a_000001 <- c(0.9992, 0.9985, 0.9970, 0.9942, 0.9888, 0.9781,
0.9572, 0.9167, 0.8410, 0.7080, 0.5038, 0.2620)
dat <- list(x = Tt_r2)
formula = formula(x~.)
frec1 = recipe(formula = formula, data = dat) |>
step_slug_cbp(
time = x,
radius = 1.0,
radius_casing = 1.0,
radius_well = 1.0,
specific_storage = 1e-1,
hydraulic_conductivity = 1.0,
thickness = 1.0,
head_0 = 1.0,
n_terms = 12L
) |>
plate("dt")
plot(slug_cbp~x, frec1, type = "l")
points(a_01~dat$x, col = "red", pch = 20L)
# estimate transmissivity with bouwer rice
data(bouwer)
bw <- as.data.frame(bouwer)
rc <- 4/2/12 # radius of 2 inches
rw <- 8.25/2/12 # radius of screen
Le <- 10.0 # screen length
y0 <- 1.44 # initial drawdown
Lw <- 17.92 # height of water above screen bottom
wl <- 6.08 # static wl
H <- 18.92 # height of water level above base
t <- as.numeric(bouwer$datetime - bouwer$datetime[1L]) # elapsed time
y <- wl - bouwer$val # change in wl
bouwer_rice(t, y, rw, rc, Le, Lw, H) * 86400 # should be ~4.5
#> [1] 4.45846References
Bouwer, H., 1989. The Bouwer and Rice Slug Test—An Update a. Groundwater, 27(3), pp.304-309.
Cooper, H.H., J.D. Bredehoeft and S.S. Papadopulos, 1967. Response of a finite-diameter well to an instantaneous charge of water, Water Resources Research, vol. 3, no. 1, pp. 263-269.
sessionInfo()
#> R version 4.6.0 (2026-04-24)
#> Platform: x86_64-pc-linux-gnu
#> Running under: Ubuntu 24.04.4 LTS
#>
#> Matrix products: default
#> BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
#> LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
#>
#> locale:
#> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
#> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
#> [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
#> [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
#> [9] LC_ADDRESS=C LC_TELEPHONE=C
#> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
#>
#> time zone: Etc/UTC
#> tzcode source: system (glibc)
#>
#> attached base packages:
#> [1] stats graphics grDevices utils datasets methods base
#>
#> other attached packages:
#> [1] splines2_0.5.4 RcppRoll_0.3.2 tibble_3.3.1 bench_1.1.4
#> [5] data.table_1.18.4 gslnls_1.4.2 collapse_2.1.7 hydrorecipes_0.0.6
#> [9] Bessel_0.7-0 rmarkdown_2.31
#>
#> loaded via a namespace (and not attached):
#> [1] gtable_0.3.6 xfun_0.57 bslib_0.11.0
#> [4] ggplot2_4.0.3 htmlwidgets_1.6.4 recipes_1.3.2
#> [7] lattice_0.22-9 vctrs_0.7.3 tools_4.6.0
#> [10] generics_0.1.4 parallel_4.6.0 pkgconfig_2.0.3
#> [13] Matrix_1.7-5 RColorBrewer_1.1-3 S7_0.2.2
#> [16] lifecycle_1.0.5 compiler_4.6.0 farver_2.1.2
#> [19] codetools_0.2-20 RcppThread_2.3.0 htmltools_0.5.9
#> [22] sys_3.4.3 buildtools_1.0.0 class_7.3-23
#> [25] sass_0.4.10 yaml_2.3.12 lazyeval_0.2.3
#> [28] gmp_0.7-5.1 profmem_0.7.0 prodlim_2026.03.11
#> [31] plotly_4.12.0 pillar_1.11.1 jquerylib_0.1.4
#> [34] tidyr_1.3.2 MASS_7.3-65 cachem_1.1.0
#> [37] gower_1.0.2 rpart_4.1.27 parallelly_1.47.0
#> [40] lava_1.9.1 tidyselect_1.2.1 digest_0.6.39
#> [43] future_1.70.0 earthtide_0.1.8 listenv_0.10.1
#> [46] dplyr_1.2.1 purrr_1.2.2 maketools_1.3.2
#> [49] splines_4.6.0 fastmap_1.2.0 grid_4.6.0
#> [52] cli_3.6.6 magrittr_2.0.5 utf8_1.2.6
#> [55] survival_3.8-6 future.apply_1.20.2 withr_3.0.2
#> [58] Rmpfr_1.1-2 scales_1.4.0 timechange_0.4.0
#> [61] lubridate_1.9.5 httr_1.4.8 globals_0.19.1
#> [64] otel_0.2.0 nnet_7.3-20 timeDate_4052.112
#> [67] evaluate_1.0.5 knitr_1.51 hardhat_1.4.3
#> [70] viridisLite_0.4.3 rlang_1.2.0 Rcpp_1.1.1-1.1
#> [73] glue_1.8.1 sparsevctrs_0.3.6 ipred_0.9-15
#> [76] jsonlite_2.0.0 R6_2.6.1