Function to calculate the drawdown using the Theis solution
theis_solution(ptest, a, t0, t)
| ptest | A pumpint 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 |
| t | Numeric vector with the time values |
A numeric vector with the drawdown
Theis, C. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage. Transactions of the American Geophysical Union, 1935, 2, 519-524
Other theis functions: exponential_integral,
theis_WF_LT_dlogt,
theis_WF_LT,
theis_calculate_parameters,
theis_solution_dlogt,
theis_solution_initial,
theis_well_function_large,
theis_well_function_small,
theis_well_function
# Data from a confined aquifer data(theis) ptest <- pumping_test('Well1', Q = 1.3e-3, r = 200, t = theis$t, s = theis$s) res <- theis_solution_initial(ptest) s_model <- theis_solution(ptest,res$a, res$t0, theis$t) #Compare drawdown from model with measurements mse <- sqrt(mean((s_model-theis$s)^2)) print(mse)#> [1] 0.1137923