Function to calculate the drawdown of a Hantush-Jacob model
hantush_jacob_solution(ptest, a, t0, lambda, t)
| ptest | A pumping_test object |
|---|---|
| a | Slope of the straight line fitted to the drawdown data using the Cooper-Jacob approach |
| t0 | Intercept of the straight line fitted to the drawdown data using the Cooper-Jacob approach |
| lambda | Leakage factor defined as $$\lambda = \sqrt{\frac{b^{'}T}{k_{h}^{'}}}$$ where \(b^{'}\) and \(k_{h}^{'}\) are the thickness and the vertical hydraulic conductivity of the aquitard, and \(T\) is the transmissivity of the underlying aquifer. |
| t | Numeric vector with the time values |
A numeric vector with the calculated drawdown
Hantush, M. S. Flow to wells in aquifers separated by a semipervious layer. Journal of Geophysical Research, 1967, 72, 1709.
Other hantush_jacob functions: hantush_jacob_WF_LT_dlogt,
hantush_jacob_WF_LT,
hantush_jacob_calculate_parameters,
hantush_jacob_solution_dlogt,
hantush_jacob_solution_initial,
hantush_jacob_well_function
data(hantush_jacob) ptest <- pumping_test('Well1', Q = 6.309e-3, r = 3.048, t=hantush_jacob$t, s = hantush_jacob$s) ptest$additional_parameters$B <- 6.096 sol0 <- hantush_jacob_solution_initial(ptest) sol <- hantush_jacob_solution(ptest,sol0$a, sol0$t0, sol0$lambda, ptest$t) print(sol)#> [1] 3.517643 6.729374 8.286243 9.312756 10.500161 12.693943 13.947944 #> [8] 15.465329 16.371902 17.261669 17.680389 18.114622 18.645372 18.958491 #> [15] 19.189459 19.480924 19.610391 19.696790 19.748772 19.796912 19.820854 #> [22] 19.827123 19.828194 19.828041 19.827638 19.826990 19.826659 19.826531 #> [29] 19.826506 19.826532 19.826556 19.826616 19.826667 19.826719 19.826751 #> [36] 19.826782 19.826812 19.826829 19.826854 19.826863 19.826872 19.826875 #> [43] 19.826876