Models and Interpretation Methods for Single-Hole Flowmeter Experiments

: Subsurface and groundwater ﬂow characterization is of great importance for various environmental applications, such as the dispersion of contaminants and their remediation. For single-hole ﬂowmeter measurements, key characteristics, such as wellbore storage, skin factor heterogeneities, and variable pumping and aquifer ﬂow rates, have a strong impact on the system characterization, whereas they are not fully considered in existing models and interpretation methods. In this study, we develop a new semi-analytical solution that considers all these characteristics in a physics-based consistent manner. We also present two new interpretation methods, the Double Flowmeter Test with Transient Flow rate (DFTTF) and the Transient Flow rate Flowmeter Test (TFFT), for interpreting data collected during single and multiple pumping tests, respectively. These solution and methods are used as follows. (i) The impact of wellbore storage, transient pumping rate, and property heterogeneities on the interpretation of data collected during single pumping tests are studied over 49 two-aquifer cases. (ii) The effect of the skin factor heterogeneity on transmissivity and storativity estimates, as well as the variability range of the (non-unique) corresponding solutions, are analyzed for the interpretation of multiple-pumping experiments. The results presented in this work show the importance of the various properties and processes that are considered, and the need for the new models and methods that are provided.


Introduction
Characterizing subsurface and groundwater flow is performed by collecting data from several techniques, including hydraulic tests, thermal experiments and electrical measurements (e.g., [1][2][3][4]), and inverting this data with the most appropriate inversion strategies (e.g., [5][6][7][8]). Among all these methods, the most-used method is the pumping test, which gives global estimates of the hydraulic properties of the system and provide information on its transient behavior when interpreting the data with transient-flow solutions (e.g., [9][10][11][12]). When characterizing the vertical heterogeneities of the system is required, for investigating the dispersion of contaminant and planning remediation strategies for instance, additional information can be acquired by downhole well logging measurements including temperature, vertical flow rates, and direct observations with optical and acoustic imaging tools (e.g., [13][14][15]).
Using heat pulse and electromagnetic tools [13,[16][17][18] results in measuring vertical flow rates that are as small as 0.05 Lpm. These measurements are used either in a qualitative way to identify the most permeable zones of the system, or in a quantitative manner to evaluate the hydraulic properties of these zones [19][20][21]. In the latter case, the properties are estimated from vertical flow rate measurements that are collected above each conductive zone in single-or cross-borehole configurations (e.g., [14,22,23]). Two main methods are used to collect this data. (i) Log measurements are acquired during a single pumping test, which is easy and quick to perform, and results in collecting one or two values of flow rates above each conductive zone while the logging tool is lowered in the well (e.g., [1,24,25]). (ii) Multiple pumping tests are performed with the logging tool being located above a different conductive zone for each test, resulting in monitoring the changes in vertical flow rates over time (e.g., [14,26,27]). Method (ii) corresponds to series of local experiments that provide more information than method (i), while being time consuming, since it requires to conduct more experiments and to wait until that the system goes back to equilibrium between each experiment.
The corresponding data are usually interpreted with standard analytical solutions that represent the identified conductive zones as independent aquifers, whose behavior is described with the Theis Equation [28][29][30]. Alternatively, double-layer numerical models are used to consider more complex configurations, including vertical crossflow between adjacent porous layers or skin effect with a homogeneous skin zone located around the well [31][32][33]. Heterogeneous skin factors, wellbore storage, and transient pumping rate are considered separately in different studies such as the semi-analytical multi-aquifer models developed in [34,35] and the methods formulated in [14,26]. These three features describe phenomena occurring in the pumping well and its vicinity that are important for the interpretation of pumping tests. (i) The skin factor enables one to take into account the skin effects that are related to headloss in the well and in its vicinity in the context of significant differences between the well vicinity and aquifer hydraulic conductivities. These differences are due to, for instance, mud invasion or turbulence in the well and its vicinity [9]. (ii) The wellbore storage impacts the data collected at the beginning of the experiment and the solutions that neglect this characteristic must be applied when the wellbore storage effect is completed. This requires to conduct long-time experiments, in particular for lowpermeability aquifers with positive skin effects and large well diameters [30,36]. (iii) In addition to intentional changes related to the experiment need, changes in pumping flow rate along pumping experiments are very frequent due to the pump functioning, with a decrease in flow rate when the drawdown increases. These changes can also be due to potential clogging of the pump or the command valve, and must be accounted for when the decrease is important [27].
Some of the solutions cited above rely on strong assumptions regarding the vertical heterogeneities of the considered systems, assuming for instance homogeneous hydraulic diffusivities and storativities in the system. Furthermore, none of these solutions considers at the same time the effect of wellbore storage, skin factor heterogeneities, and variable pumping and aquifer flow rates, whereas the importance of each of those characteristics has been demonstrated separately. In this study, we develop a new semi-analytical solution that considers all these key characteristics and that is used to generate synthetic reference data. We also remind the assumptions related to standard interpretation methods and we analyze their impact on the interpretation of single-hole flowmeter measurements that are collected over single pumping tests. For the interpretation of flowmeter data collected during single or double logs along the well, a new interpretation method is presented to improve the estimation of aquifer properties, and evaluate the impact of standard assumptions over a large range of cases. For the interpretation of flowmeter data collected during multiple pumping tests, where the full transient behavior of vertical flow rates above each conductive zone is recorded, we also present a new interpretation method that is easy to implement. This method is used to evaluate the variability range of key properties that drive the transient behavior of the system, which are the storativity and skin factor. The validity and accuracy of each model and interpretation method are discussed through synthetic cases, before providing a discussion and conclusions on this work.

Considered Experiments and Methods
The analysis of flowmeter tests provided in this study is performed by considering the following kinds of experiments (Table 1).

•
Single-pumping single-log experiments. The vertical flow rates and hydraulic heads are measured along the borehole while the flowmeter is lowered into the well during a single pumping test. • Single-pumping double-log experiments. As before, the vertical flow rates and hydraulic heads are measured along the borehole during two log experiments that are conducted under the same pumping test. • Multiple-pumping local-log experiments. A pumping test is performed for each conductive zone that needs to be characterized (except the upper one) with the logging tool localized above this zone. Table 1. Considered flowmeter experiments with the following collected data: h w,i and q i are the hydraulic heads and vertical flow rates, respectively, measured above conductive zone i with N aq the number of conductive zones to characterize, and h j w,i and q j i are their counterpart collected from two logs (j = 1, 2). SFT, DFT, DFTTF, and TFFT are the interpretation methods presented in Appendix A and Section 2.

Experiment Name Collected Data Interp. Methods
Single-pumping single-log Multiple-pumping local-log (h w,i (t), q i (t)) TFFT As described in Table 1, these experiments are interpreted with the methods SFT, DFT, DFTTF, and TFFT. SFT and DFT are standard methods whose basis are reminded in Appendix A, and DFTTF and TFFT correspond to new methods that are described below. For the methods presented below, the vertical flow rates q i (i = 1, . . . , N aq ) measured above aquifer i are expressed in terms of aquifer flow rates Q i (i = 1, . . . , N aq ) with the differentiation method described in expressions (A2).

Double Flowmeter Test with Transient Flow Rate (DFTTF)
The Double Flowmeter Test with Transient Flow rate (DFTTF) method is based on the DFT, which consists of evaluating the transmissivity T i and storativity S i of aquifer i from Theis' solution using two measurements of the vertical flow rates and hydraulic head above each conductive zone (see the description of the DFT method in Appendix A.2). While the DFT method relies on the assumption that the values of the considered vertical flow rates are equal to their average value, the DFTTF method considers these values as distinct without assumption. In the conditions of validity of the logarithmic approximation of the Theis function, the hydraulic properties of aquifer i are expressed as: with the aquifer flow rates Q i (t i1 ) and Q i (t i2 ) and hydraulic heads h wi (t i1 ) and h wi (t i2 ) deduced from measurements performed at times t i1 and t i2 for aquifer i. In expression (1), r wi and σ i are the well radius and skin factor associated with aquifer i, respectively, σ i being considered as homogeneous over the aquifers (i.e., σ i = σ) and defined from pump test interpretation, as in DFT. Note that in expression (1b), S i is expressed as a function of t i1 and T i , but could also be expressed as a function of t i2 and T i by replacing t i1 with t i2 .
Since T i depends on both t i1 and t i2 , S i also depends on both times regardless of the chosen expression.

Transient Flow Rate Flowmeter Test (TFFT)
The TFFT (Transient Flow rate Flowmeter Test) method is developed to interpret multiple-pumping local-log experiments. These experiments consist of performing a pumping test for each aquifer to characterize (except the upper one) with the same pumping rate Q P . The resulting drawdown h wS and the full transient behavior of the vertical flow rates q i above aquifers i are recorded and the corresponding aquifer transient flow rates Q i are deduced as explained before. This leads to the set of data (Q i (t), Q P (t), h wS (t)) associated with each aquifer i (i = 1, . . . , N aq ) and collected during N P = N aq − 1 pumping tests. These data are inverted with the following algorithm.

1.
Model the data Q P (t) with variable pumping flow rate models, as described in Appendix C.1.

2.
Estimate the unknowns T i , S i , σ i by inverting the data (Q i (t), h wS (t)) provided by the multi-aquifer model described in Appendix B using a numerical optimization method. More precisely, the Laplace transform of Q i (t) and h wS (t) are given in expressions (A15) and (A17) and numerically inverted with [37]'s algorithm.
The data inversion proceeds iteratively by a least square method using a gradient algorithm, applied to the transient drawdown and aquifers flow rates, for which different weights can be considered. The objective function is expressed as: with the weights α t and α t i , and the reference (Q data Note that this objective function is used in Section 3.3 with synthetic data that are considered every minute and weights that are set to 1. The inversions are performed from an initial value given by the DFTTF method with null skin factors and convergence is assumed when the relative change of the objective function is less than 10 −12 in the last five iterations.

Considered Configurations
We consider a system of two aquifers, T 1 and S 1 being the transmissivity and storativity of aquifer 1, respectively, and T 2 and S 2 of aquifer 2. We evaluate how these properties are estimated from single-and double-log flowmeter experiments with the interpretation methods SFT, DFT, and DFTTF in various context in terms of wellbore storage and transient pumping rate (Section 3.2). Then, we evaluate these properties by taking into account different skin factors in two-aquifer systems considering wellbore storage and transient pumping rate with an analysis of the range of the storativity and skin factor that can be defined (Section 3.3). This is done by considering series of local flowmeter experiments that are interpreted with the TFFT method. For all the cases presented in this section, the reference synthetic data are provided by the multi-aquifer model that is presented in Appendix B and takes into account the wellbore storage effect, transient pumping flow rate, and heterogeneous skin factor. A homogeneous pumping well of radius 8 cm is considered with aquifers of 1 m thickness, resulting in transmissivity and storativity equal to hydraulic conductivity and specific storage, respectively.

Single-and Double-Log Flowmeter Experiments
The considered two-aquifer system is defined by setting T 1 and S 1 to 10 −4 m 2 /s and 10 −3 , respectively, and having T 2 and S 2 ranging from 10 −7 to 10 −1 m 2 /s and 10 −6 to 1, respectively, with null skin factor in both aquifers. The three following configurations are considered: (i) standard models without accounting for wellbore storage and considering constant pumping flow rate (Config1), (ii) models that account for wellbore storage and consider constant pumping flow rate (Config2), and (iii) no wellbore storage and transient pumping flow rate (Config3). The transmissivities are estimated with the interpretation methods SFT, DFT, and DFTTF, and the storativities with DFT and DFTTF, considering the measurements at time t 1 = 10 min for SFT and at times t 1 and t 2 = 300 min for DFT and DFTTF. Note that the small time t 1 allows to investigate the influence of the wellbore storage and the transient pumping rate. The corresponding results are presented in Figure 1 and the pumping flow rate Q P that is used in Config3 is defined in Figure 2. Note that Q P is set to the constant value 4 Lpm in Config1 and Config2.
The results provided in Figure 1 are described as follows. (i) When the diffusivities of aquifers i (aquifer for which the properties are estimated) and j (the other aquifer) are similar (i.e., ε ij ≈ 1), the ratio of estimated to true properties (T i and S i ) are equal to 1 (or very close to 1) for all the considered configurations and interpretation methods. This shows that for homogeneous hydraulic diffusivity (i.e., ε i ≈ ε j ), the assumptions related to the considered interpretation methods hold, in particular the assumptions related to Theis' model. (ii) For large values of ε ij (i.e., ε ij >> 1), the errors in estimating the hydraulic properties are relatively small in comparison with the errors obtained for small values of ε ij (i.e., ε ij << 1). In most of the cases, we observe that the error in estimating T i and S i increases when ε ij varies from 1 to 10 3 because the assumptions described before that hold when ε ij ≈ 1 are less and less fulfilled. This error then decreases when ε ij varies from 10 3 to 10 6 , which corresponds to configurations where the diffusivity of aquifer j is negligible in comparison with that of aquifer i for which the properties need to be estimated. This implies that aquifer j poorly contributes to the pumping and that the assumption of constant pumping related to aquifer i and used in Theis' model holds. (iii) For small values of ε ij (i.e., ε ij << 1), the transmissivities and storativities are overestimated and underestimated, respectively. ε i < ε j corresponds to configurations where T i < T j or S i > S j , whereas the methods SFT and DFT rely on the assumption of homogeneous diffusivities. Overestimating T i and underestimating S i result in tending to the relations T i ≈ T j and S i ≈ S j , and thus, tending to verify the homogeneous diffusivity assumption.
(iv) For all the estimates of T 1 and T 2 , the results obtained with SFT and DFT are similar or slightly improved by DFT, except for the estimated value of T 1 in Config3. For all the cases, the DFTTF method always provides better estimates than SFT and DFT with the ratio of estimated to true value around 1 for most of the cases and reaching for example a maximum value of 3 for the transmissivity estimates. (v) Only small differences are observed between the three considered configurations, showing that the wellbore storage and transient pumping flow rate do not have a significant impact on the estimation of the considered properties. An exception is observed when estimating the storativities S 1 and S 2 in Config3, for which the transient pumping flow rate results in increasing the ratio of estimated to the true value of 4 to 6 orders of magnitude with the DFT method (blue symbols in Figure 1i). This results in a significant underestimation of S 1 and S 2 .

Series of Local Flowmeter Experiments
We consider now the effect of different skin factors in two-aquifer configurations with T 1 and T 2 set to 5 × 10 −4 and 10 −5 m 2 /s, S 1 and S 2 to 5 × 10 −4 and 10 −3 , and the skin factors σ 1 and σ 2 to 0 and 1, respectively. As before, the pumping well radius is homogeneous and equal to the wellbore storage radius r wS = 8 cm. The applied transient pumping flow rate is the same as in Config3 of the previous example and the resulting transient aquifer flow rates are defined from the multi-aquifer model provided in Appendix B and Figure 2. These data are interpreted with the TFFT method presented in Section 2.3 and the resulting estimated values are presented in Figure 3. with Q t 1 = 4 and Q t 2 = 3.8 Lpm at times t 1 = 0 and t 2 = 300 min with the fitting coefficient b set to 10 5 . Q P is the pumping flow rate considered in Config3 of Section 3.2 and in all the experiments in Section 3.3. Q t 1 and Q t 2 correspond to the aquifer flow rates of the reference experiments considered in Section 3.3. Figure 3a shows the minimum value of the objective function (2) obtained for each value of couple (σ 1 , σ 2 ). These results show that the smallest value is observed for the true values of σ 1 and σ 2 (i.e., σ 1 = 0 and σ 2 = 1), demonstrating that the studied objective function is well defined. The corresponding estimated values of T 1 , T 2 , S 1 , and S 2 provided in   Figure 3d,e show minimum localized values of both S 1 and S 2 for high contrast between σ 1 and σ 2 . This corresponds to configurations where σ 1 is large and σ 2 small, or σ 1 small and σ 2 large, which are represented in the top right and bottom left corners of Figure 3d,e (yellow and orange cells). (iv) We also observe that high values of σ 1 and σ 2 result in overestimating S 1 and S 2 with red cells located at the bottom right corners of Figure 3d,e. High values of σ i correspond to high and low values of the headlosses in the skin and aquifer, respectively, corresponding to a low contribution of the aquifer to the drawdown. Trying to fit reference data that are obtained with smaller values of σ i results in overestimating the storativities to counterbalance the impact of large values of σ i .
The storativity is usually determined from the drawdown recorded in a distant observation well. Determining this property from data collected in the pumping well is a challenge, since an infinity of couples (S, σ) can reproduce the recorded drawdown, as demonstrated in Appendix C.2. The resulting relationship (A23) between S i and σ i is used to estimate an equivalent couple of parameters (S i , σ i ) that are presented in Figure 3f,g. The global behavior of S 1 and S 2 is well reproduced, except for the localized minimum values that are observed at the bottom left corners in Figure 3d,e, and not in Figure 3f,g. These differences are due to relationship (A23), which only depends on (S i , σ i ) when estimating (S i , σ i ), implying that this relationship does not consider the impact of one aquifer on the other. These results show that the localized minimum values observed in Figure 3d,e are due to the impact of the properties of aquifer j when estimating the storativity of aquifer i.

Discussion
We present a new semi-analytical multi-aquifer model relying on an independent aquifers representation. This physics-based solution is formulated for heterogeneous aquifer properties, skin factor, and pumping well radius, taking into account the wellbore storage and various transient pumping rates and transient aquifer flow rates. We also present two new interpretation methods that (i) improve the transmissivity and storativity estimates of multi-aquifer systems, while taking into account wellbore storage and transient pumping flow rates, and (ii) help to conduct sensitivity analysis of these estimates in relation to the skin factor.
Our study shows that the interpretation of single-and double-log flowmeter experiments collected during a single pumping test is improved by using the DFTTF method. By accounting for transient aquifer flow rates, this method gives better estimates of the transmissivity and storativities of aquifers than the standard SFT and DFT methods. We also show that the interpretation of multiple-pumping local-log experiments conducted with the TFFT method leads to consistent results. It demonstrates radically different sensitivity of the transmissivities and storativities to the skin factors considered in the aquifers. A relationship between storativity and skin factor is provided and tested to analyze the (nonunique) couples of solutions related to these parameters in the context of single-hole data, while the storativity is usually estimated from data collected in an observation borehole.
However, the presented solution and methods have been only applied to two-aquifer systems and synthetic data with relatively simple configurations. Additional work is required to demonstrate the efficiency of this solution and methods on complex configurations and their ability to interpret field data. Extending the sensitivity analysis also presented in this work to complex configurations and field data is required as well to fully demonstrate the interest of this work.

Conclusions
For future work and applications, this easy-to-implement solution can be extended to account for inwell headlosses, and for interpreting the recovery phase of flowmeter experiments by applying the superposition principle. We will also consider structural information that helps to reduce the range of acceptable skin factors and associated storativities. For example, small and large storativities are unlikely in semi-confined systems and in weakly fissured confined hard rocks, respectively. Using logs of acoustic or optical borehole images helps to verify the presence of fractures, which gives information on the skin factor. The latter is negative when the well is intersected by local well-open conducts (such as well-open fractures), small with low flow rates, null when considering fractures with homogeneous aperture and rugosity in laminar regime, and positive for fractures filled by drilling mud or sediments. The relevance of the skin factor value, which corresponds to a singular headloss at the infinitesimal interface between the well and the aquifer, can also be checked with equivalent finite skins, such as darcian skins, or laminar-turbulent fracture skins [38].

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Standard Models and Interpretation Methods
Considering a homogeneous and infinite aquifer subject to a constant pumping flow rate Q with negligible wellbore storage, ref. [39]'s solution provides the following expression of the well drawdown: where K, S s , and b are the aquifer hydraulic conductivity, specific storativity, and thickness, respectively, σ and r w are the skin factor and pumping well radius, and W is the Theis function. In standard pumping tests, drawdown data are collected in a distant observation well and expression (A1) is used to estimate global values of the aquifer transmissivity and storativity that are defined as T = Kb and S = S s b.
In a system composed of horizontal and independent aquifers, where T i and S i are the transmissivity and storativity of aquifer i, single-hole flowmeter tests rely on logs that are performed along the pumping well and provide the vertical flow rates q i measured above each conductive zone i. The aquifer flow rate Q i of aquifer i, which is defined as the flow rate provided by the aquifer during the pumping experiment, is deduced from the measurements q i with the following vertical differentiation: where N aq is the number of horizontal aquifers numbered from top to bottom. Estimates of the hydraulic properties of each aquifer i are provided by the SFT and DFT interpretation methods for single and double flowmeter tests, each of those methods relying on different assumptions that are described below.
Appendix A.1. Single Flowmeter Test (SFT) The SFT method consists of interpreting aquifer flow rates Q i deduced from a singlelog of vertical flow rates q i and associated with the hydraulic head h w,i for each conductive zone i. The hydraulic transmissivity T i is estimated at a given time using the following assumptions: (H1) the hydraulic diffusivities ε i = K i /S i are homogeneous (i.e., ε i = ε, ∀i) [40] and (H2) the specific storativities S si are homogeneous (i.e., S si = S s , ∀i) [28]. Homogeneous skin factor is also usually assumed (i.e., σ i = σ) and T i is inverted from (A1) for each aquifer i, analytically under assumption (H1) and numerically under (H2). In our study, we consider the SFT method with assumption (H1). Alternatively, considering in addition negligible inwell headlosses (i.e., h wi (t) = h w (t)) and homogeneous pumping well radius (i.e., r wi = r w ), T i is given by the well-know formula T i = TQ i /Q. This expression has been deduced from numerical simulations in [29] for horizontal layered systems with cross-flow exchanges under pseudo-steady state conditions, and from analytical demonstration in [30], as done here for independent layered systems under transient conditions.

Appendix A.2. Double Flowmeter Test (DFT)
Alternatively, the DFT method relies on two logs of data which provide the aquifer flow rates Q i (t i1 ) and Q i (t i2 ) at two distinct times t i1 and t i2 . Assuming homogeneous skin factors as before, and being in the conditions of validity of the logarithmic approximation of the Theis function [39], the following explicit expressions of T i and S i are obtained: [30] whereQ i is the average between flow rates Q i (t i1 ) and Q i (t i2 ) (according to Theis' assumption of constant flow rate) and r we is the effective radius [41] which integrates the skin factor σ with expression r we = r w e −σ . For this method, and for all the double-log methods presented in this work, if the logarithmic approximation does not hold, T i and S i can be inverted jointly from (A1) applied to aquifer i at times t i1 and t i2 by using a numerical optimization method. However, as stated in [39], the logarithmic approximation is quickly reached in the pumping well, except in rare particular cases.

Appendix B. Multi-Aquifer Model
We present a semi-analytical solution for simulating transient drawdown in radially infinite confined multi-aquifers with homogeneous boundary conditions, wellbore storage, heterogeneous skin effects, and variable pumping flow rate. The presented model also accounts for a variable well radius, which is useful for deep wells in which the well radius decreases with depth. We consider N aq horizontal independent aquifers numbered from top to bottom ( Figure A1) that are intercepted by the pumping well. When applying the pumping flow rate Q P to the well, mass balance in the well is expressed as: where Q P is negative when the flow is extracted from the well, Q i is the (positive) flow rate pumped from aquifer i, and Q S is the (negative) wellbore storage flow rate. Note that Q P here is equivalent to −Q in Appendix A.1. The flow rates used in Equation (A4) are expressed as: where S w is the well capacity defined as S w = πr 2 wS with r wS the wellbore storage radius (i.e., the radius of the zone where the water table fluctuates), h wS is the wellbore storage drawdown (i.e., the water surface level), r wi is the well radius at the depth of aquifer i, b i and K i are the conductive thickness and conductivity of aquifer i, respectively, and h i is the drawdown in aquifer i. In each aquifer i (i = 1, . . . , N aq ), the mass balance equation for a slightly compressible fluid is expressed as: where S si is the specific storage of aquifer i, r is the radial distance from the well center, and t is the time elapsed since the pumping starts. Initial and boundary conditions are defined by: and the well skin effect is represented by a singular head loss as [42]: with the dimensionless skin factor σ i . Applying the Laplace transform to Equation (A6) with the initial condition provided in (A7) leads to: where p is the Laplace variable and ε i is the diffusivity defined as ε i = K i /S si . Considering the boundary condition given in (A7), the solution for Equation (A9) is expressed as: where I and K are the modified Bessel functions of the first and second kind, respectively, γ i = p/ε i , and C 1i = 0 when applying the boundary condition from (A7). The derivative of h i is then expressed in the Laplace domain as: Expression (A8) is formulated in the Laplace domain as: and expressed as follows by using (A10) and (A11): h wi = C 2i A 1i , A 1i = K 0 (γ i r wi ) + r wi σ i γ i K 1 (γ i r wi ).
The following expressions of h i and its derivative in the Laplace domain are deduced from (A13):h Expressing the Laplace transform of the vertical flow rate Q i , defined in (A5), and combining it with the derivative expression in (A14) result in: Combining (A4), (A5) and (A15) result in: and considering that the headloss in the well is negligible (i.e.,h wi =h wS ),h wS is expressed as: All the drawdowns and flow rates expressed in the Laplace domain are inverted in the time domain using [37]'s algorithm.

Appendix C. Models for Transient Parameters and Properties
Appendix C. 1

. Pumping Flow Rate Models
We focus on expressions of flow rate models that have an explicit analytical formulation in the Laplace domain, for example exponential and polynomial formulations, as well as linear combinations of them. As the exponential decrease reproduces well the typical flow rate decrease that is observed due to the pump functioning [27], we consider this formulation with a decrease from flow rates Q t 1 to Q t 2 at times t 1 and t 2 , respectively. This leads to the following pumping flow rate expression: where b is a fitting coefficient that controls the decrease shape. The Laplace transform is expressed as: considering that the pumping starts at time t 1 = 0.
Appendix C.2. Couples of Equivalent Parameters (S si , σ i ) For single-hole tests, the couples of parameters (S si , σ i ) and (S si , σ i ) are equivalent when they produce the same values of h wS , Q S and Q i . From the demonstration of the multi-aquifer model provided in Appendix B, this is equivalent to produce the same values of A 2i and A 2i that are defined as: A 2i = 2πr wi b i K i γ i K 1 (γ i r wi ) K 0 (γ i r wi ) + r wi σ i γ i K 1 (γ i r wi ) , (A20) A 2i = 2πr wi b i K i γ i K 1 (γ i r wi ) K 0 (γ i r wi ) + r wi σ i γ i K 1 (γ i r wi ) where γ i = pS si /K i . Considering A 2i = A 2i results in: Focusing on small arguments of the Bessel functions, the following approximations [43] are used: where ln is the natural logarithm and γ E is the Euler constant. This leads to the following explicit relationship between the couples of properties (S si , σ i ) and (S si , σ i ):