Effect of wellbore storage and finite thickness skin on flow to a partially penetrating well in a phreatic aquifer

An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.     

Semi‐analytical solution for a slug test in partially penetrating wells including the effect of finite‐thickness skin

This paper presents a new semi‐analytical solution for a slug test in a well partially penetrating a confined aquifer, accounting for the skin effect. This solution is developed based on the solution for a constant‐flux pumping test and a formula given by Peres and co‐workers in 1989. The solution agrees with that of Cooper and co‐workers and the KGS model when the well is fully penetrating. The present solution can be applied to simulate the temporal and spatial head distributions in both the skin and formation zones. It can also be used to demonstrate the influences of skin type or skin thickness on the well water level and to estimate the hydraulic parameters of the skin and formation zones using a least‐squares approach. The results of this study indicate that the determination of hydraulic conductivity using a conventional slug‐test data analysis that neglects the presence of a skin zone will give an incorrect result if the aquifer has a skin zone. Copyright © 2008 John Wiley & Sons, Ltd.     

Non-linear laminar flow into well partially penetrating a porous aquifer of finite thickness

In steady state condition, non-linear laminar flow of fluid into a well partially penetrating a porous aquifer of finite thickness is considered. The influence of such a flow on discharge and its dependence on related physical quantities are investigated. It is observed that the discharge into the well decreases as the depth of the well is decreased and the region of non-linear laminar flow is widened, which is quite obvious from physical considerations. As a particular case, result for a fully penetrating well has been deduced.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

Analysis of the water content in an unconfined aquifer with a partially penetrated pumping well

Most methods developed to represent water flow phenomena in an unconfined aquifer with a fully penetrated pumping well are either numerical, such as the well‐known FEMWTER model, or experimental; analytical models of a partially penetrated pumping well are rare. This study employs the linearized Richards equation as the governing equation, with the aid of Fourier Integral Transformation, to obtain an analytical solution of the water content distribution in an unconfined aquifer with a partially penetrated pumping well. The results from this study could serve to substantiate in some sense results from numerical models. In addition, the theory developed herein can be modified to simulate a vacuum‐pressured pumping well since it is derived by considering, among others, the location and length of a well screen with fluxes. Copyright © 2008 John Wiley & Sons, Ltd.     

Two-region non-Darcian flow toward a well in a confined aquifer

We have derived an analytical solution for two-region flow toward a well in a confined aquifer based on a linearization method. The two-region flow includes Izbash non-Darcian flow near the well and Darcian flow in the rest of the aquifer. The wellbore storage is also considered. The type curves in the non-Darcian and Darcian flow domains are obtained by a numerical Laplace inversion method incorporated in MATLAB programs. We have compared our results with the one-region Darcian flow model (Theis). Our solutions agree with those of Sen [Sen Z. Type curves for two-region well flow. J Hydr Eng 1988;114(12):1461–84] which were obtained using the Boltzmann transform at late times for fully turbulent flow, while some difference has been found at early and moderate times. We have defined a dimensionless non-Darcian hydraulic conductivity term which is shown to be a key parameter for analyzing the two-region flow. A smaller dimensionless non-Darcian hydraulic conductivity results in a larger drawdown in the non-Darcian flow region at late times. However, the dimensionless non-Darcian hydraulic conductivity does not affect the slope of the dimensionless drawdown versus the logarithmic dimensionless time in the non-Darcian flow region at late times. The dimensionless non-Darcian hydraulic conductivity does not affect the late time drawdown in the Darcian flow region.     

Radial groundwater flow to a finite diameter well in a leaky confined aquifer with a finite‐thickness skin

A mathematical model that describes the drawdown due to constant pumpage from a finite radius well in a two‐zone leaky confined aquifer system is presented. The aquifer system is overlain by an aquitard and underlain by an impermeable formation. A skin zone of constant thickness exists around the wellbore. A general solution to a two‐zone leaky confined aquifer system in Laplace domain is developed and inverted numerically to the time‐domain solution using the modified Crump (1976) algorithm. The results show that the drawdown distribution is significantly influenced by the properties and thickness of the skin zone and aquitard. The sensitivity analyses of parameters of the aquifer and aquitard are performed to illustrate their effects on drawdowns in a two‐zone leaky confined aquifer system. For the negative‐skin case, the drawdown is very sensitive to the relative change in the formation transmissivity. For the positive‐skin case, the drawdown is also sensitive to the relative changes in the skin thickness, and both the skin and formation transmissivities over the entire pumping period and the well radius and formation storage coefficient at early pumping time. Copyright © 2009 John Wiley & Sons, Ltd.     

Groundwater response to tidal fluctuations in a leaky confined coastal aquifer with a finite length

The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.     

Analytical model for fully three‐dimensional radial dispersion in a finite‐thickness aquifer

Analytical solutions for contaminant transport in a non‐uniform flow filed are very difficult and relatively rare in subsurface hydrology. The difficulty is because of the fact that velocity vector in the non‐uniform flow field is space‐dependent rather than constant. In this study, an analytical model is presented for describing the three‐dimensional contaminant transport from an area source in a radial flow field which is a simplest case of the non‐uniform flow. The development of the analytical model is achieved by coupling the power series technique, the Laplace transform and the two finite Fourier cosine transform. The developed analytical model is examined by comparing with the Laplace transform finite difference (LTFD) solution. Excellent agreements between the developed analytical model and the numerical model certificate the accuracy of the developed model. The developed model can evaluate solution for Peclet number up to 100. Moreover, the mathematical behaviours of the developed solution are also studied. More specifically, a hypothetical convergent flow tracer test is considered as an illustrative example to demonstrate the three‐dimensional concentration distribution in a radial flow field. The developed model can serve as benchmark to check the more comprehensive three‐dimensional numerical solutions describing non‐uniform flow contaminant transport. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of solute transport in a divergent flow tracer test with scale‐dependent dispersion

It has been known for many years that dispersivities increase with solute displacement distance in a subsurface. The increase of dispersivities with solute travel distance results from significant variation in hydraulic properties of porous media and was identified in the literature as scale‐dependent dispersion. In this study, Laplace‐transformed analytical solutions to advection‐dispersion equations in cylindrical coordinates are derived for interpreting a divergent flow tracer test with a constant dispersivity and with a linear scale‐dependent dispersivity. Breakthrough curves obtained using the scale‐dependent dispersivity model are compared to breakthrough curves obtained from the constant dispersivity model to illustrate the salient features of scale‐dependent dispersion in a divergent flow tracer test. The analytical results reveal that the breakthrough curves at the specific location for the constant dispersivity model can produce the same shape as those from the scale‐dependent dispersivity model. This correspondence in curve shape between these two models occurs when the local dispersivity at an observation well in the scale‐dependent dispersivity model is 1·3 times greater than the constant dispersivity in the constant dispersivity model. To confirm this finding, a set of previously reported data is interpreted using both the scale‐dependent dispersivity model and the constant dispersivity model to distinguish the differences in scale dependence of estimated dispersivity from these two models. The analytical result reveals that previously reported dispersivity/distance ratios from the constant dispersivity model should be revised by multiplying these values by a factor of 1·3 for the scale‐dependent dispersion model if the dispersion process is more accurately characterized by scale‐dependent dispersion. Copyright © 2007 John Wiley & Sons, Ltd.     

Wellbore flow‐rate solution for a constant‐head test in two‐zone finite confined aquifers

The solution describing the wellbore flow rate in a constant‐head test integrated with an optimization approach is commonly used to analyze observed wellbore flow‐rate data for estimating the hydrogeological parameters of low‐permeability aquifers. To our knowledge, the wellbore flow‐rate solution for the constant‐head test in a two‐zone finite‐extent confined aquifer has never been reported so far in the literature. This article is first to develop a mathematical model for describing the head distribution in the two‐zone aquifer. The Laplace domain solutions for the head distributions and wellbore flow rate in a two‐zone finite confined aquifer are derived using the Laplace transform, and their corresponding time domain solutions are then obtained using the Bromwich integral method and residue theorem. These new solutions are expressed in terms of an infinite series with Bessel functions and not straightforward to calculate numerically. A large‐time solution for the wellbore flow rate is therefore developed by employing the relationship of small Laplace variable versus large time variable and L'Hospital's rule. The result shows that the large‐time solution is identical to the steady‐state solution obtained after applying the Tauberian theorem into the Laplace domain solution. This large‐time solution can reduce to the Thiem equation in the case of no skin. Finally, the newly developed solution is used to investigate the effects of outer boundary distance and conductivity ratio on the wellbore flow rate. Copyright © 2011 John Wiley & Sons, Ltd.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Non-Darcian flow to a well in an aquifer–aquitard system

In this study, we use a linearization procedure and a finite difference method to solve non-Darcian flow to a well in an aquifer–aquitard system. The leakage effect is considered. Flow in the aquifer is assumed to be non-Darcian and horizontal, whereas flow in the aquitard is assumed to be Darcian and vertical. The Izbash equation [Izbash SV. O filtracii V Kropnozernstom Materiale. USSR: Leningrad; 1931 [in Russian]] is employed to describe the non-Darcian flow. The wellbore storage is also considered in this study. An approximate semi-analytical solution has been obtained by the linearization procedure, and a numerical solution has been obtained by using a finite difference method. The previous solutions for Darcian flow case and non-Darcian flow case without leakage can be described as special cases of the new solutions. The error caused by the linearization procedure has also been analyzed. The relative error caused by the linearization procedure is nearly 100% at early times, and decreases to zero at late times. We have also compared the results in this study with Wen et al. [Wen Z, Huang G, Zhan H. A numerical solution for non-Darcian flow to a well in a confined aquifer using the power law function. J Hydrol, 2008d [in revision]] in which the leakage effect is not considered, and Hantush and Jacob [Hantush MS, Jacob CE. Non-steady radial flow in an infinite leaky aquifer. Trans Am Geophys Union 1955;36(1):95–100] who investigated a similar problem in Darcian flow case. The comparison of this study and Wen et al. (2008d) indicates the dimensionless drawdown in the aquifer with leakage is less than that without leakage, and the leakage has little effect at early times. The comparison between the results of this study and that of Hantush and Jacob (1955) indicates that the dimensionless drawdown in the aquifer for non-Darcian flow is larger at early times and smaller at late times, than their counterparts for Darcian flow. A larger dimensionless non-Darcian conductivity k_D results in a smaller dimensionless drawdown in the aquifer at late times, and leads to a larger dimensionless drawdown in the aquifer at early times. A smaller dimensionless leakage parameter B_D results in a smaller drawdown at late times, and the leakage does not affect the early-time drawdown. The analysis of the dimensionless drawdown inside the well has also been included in this study when the wellbore storage is considered.     

Tide‐induced groundwater head fluctuation in a coastal aquifer system with a submarine outcrop covered by a thin silt layer

We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 John Wiley & Sons, Ltd.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Mixing at the interface between fresh and salt waters in 3D steady flow with application to a pumping well in a coastal aquifer

We consider 3D steady flow of fresh water over a salt water body in a confined aquifer of constant thickness D, with application to a pumping well in a coastal aquifer. With neglect of mixing, a sharp interface separates the two fluid bodies and an existing analytical solution, based on the Dupuit assumption, is adopted. The aim is to solve for the mixing between the fresh and salt waters for α_T/D  1 (α_T transverse dispersivity), as field studies indicate that α_T = O(10⁻³ − 10⁻² m). The mixing zone around the interface is narrow and solutions by existing codes experience numerical difficulties. The problem is solved by the boundary layer (BL) approximation, extending a method, applied previously to two-dimensional flows. The BL equations of variable-density flow are solved by using the Von Karman integral method, to determine the BL thickness and the rate of entrainment of salt water along the interface. Application to the pumping well problem yields the salinity of the pumped water, as function of the parameters of the problem (well discharge, seaward discharge, well distance from the coast and density difference).     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal