Automated Estimation of Aquifer Parameters from Arbitrary‐Rate Pumping Tests in Python and MATLAB

Abstract

Inspired by the analysis by Mishra et al. (2012) of variable pumping rate tests using piecewise‐linear reconstructions of the pumping history, this article contains a derivation of the convolutional form of pumping tests in which the pumping history may take any possible form. The solution is very similar to the classical Theis (1935) equation but uses the Green's function for a pumped aquifer given by taking the time derivative of the well function . This eliminates one integration inside another and renders the convolution including the pumping history about as computationally demanding as calculating the well function alone, so that the convolution can be completed using handy mathematical software. It also allows nonlinear well losses, and because an easily‐computed deterministic model exists for all data points and pumping history, an objective function may include all data, so that errors are reduced in calculating any nonlinear‐well losses. In addition, data from multiple observation wells may be used simultaneously in the inversion. We provide codes in MATLAB and Python to solve for drawdown resulting from an arbitrary pumping history and compute the optimal aquifer parameters to fit the data. We find that the subtleties in parameter dependencies and constructing an appropriate objective function have a substantial effect on the interpreted parameters. Furthermore, the optimization from step‐drawdown tests is typically nonunique and strongly suggests that a Bayesian inversion should be used to fully estimate the joint probability density of the parameter vector.

Analysis of Groundwater Time Series With Limited Pumping Information in Unconfined Aquifer: Response Function Based on Lagging Theory

“…For any time‐variant pumping rate, it can be evaluated using convolution: $$${h}_{p}=\int \nolimits_{0}^{t}Q(\tau ){\theta }_{p}\left(t-\tau \right)d\tau $$$ with $$${\theta }_{p}(t)=\frac{d{{\Theta }}_{p}(t)}{dt}$$$ where θ p ( t ) is the so‐called impulse response function, which be readily obtained using numerical differentiation techniques. To evaluate the convolution integrals Equation based on finite discretization, the fastest convolution method uses Fast Fourier Transform, which evenly discretizes the time over which the convolution is performed (Benson, 2023).…”

“…where θ p (t) is the so-called impulse response function, which be readily obtained using numerical differentiation techniques. To evaluate the convolution integrals Equation 8a based on finite discretization, the fastest convolution method uses Fast Fourier Transform, which evenly discretizes the time over which the convolution is performed (Benson, 2023).…”

Weighted joint curve inversion of aquifer parameters in pumping tests with multiple observation wells

“…Interdisciplinary research has developed other methods that have been used to determine aquifer parameters, including hydrochemical dynamics (Cao et al 2009), inversion methods based on groundwater microdynamics (Hsieh et al 1987;Zhang et al 2019), hydraulic tomography (Illman et al 2007), and intelligent optimization algorithms (Lu et al 2011). However, these interdisciplinary methods either require data that are difficult to obtain or substantial computational effort, thereby limiting their practical application despite their theoretical foundations (Benson 2024). Currently, field pumping tests are the main method used to acquire aquifer parameters (Avci et al 2013;Houben et al 2022).…”