Modelling boreholes in complex heterogeneous aquifers

Reliable estimates of the sustainable yield of supply boreholes are critical to ensure that groundwater resources are managed sustainably. Sustainable yields are dependent on the pumped groundwater level in a borehole, how this relates to vertical aquifer heterogeneity, and features of the borehole itself. This paper presents a 3D radial flow model (SPIDERR), based on the Darcy-Forchheimer equation, for simulating the groundwater level response in supply boreholes in unconfined, heterogeneous aquifers. The model provides a tool for investigating the causes of non-linear behaviour in abstraction boreholes, which can have a significant impact on sustainable yields. This is demonstrated by simulating a variable-rate pumping test in a Chalk abstraction borehole. The application suggests the non-linear response to pumping is due to a combination of factors: a reduction in well storage with depth due to changes in the borehole diameter, a reduction in hydraulic conductivity with depth, and non-Darcian flow.

[1]  L. V. Beek,et al.  Water balance of global aquifers revealed by groundwater footprint , 2012, Nature.

[2]  Guanhua Huang,et al.  Non-Darcian flow to a well in a leaky aquifer using the Forchheimer equation , 2011 .

[3]  Z. Şen.  Analytical solution incorporating nonlinear radial flow in confined aquifers , 1988 .

[4]  L. Lebbe,et al.  MODFLOW procedure to simulate axisymmetric flow in radially heterogeneous and layered aquifer systems , 2014, Hydrogeology Journal.

[5]  B. Connorton,et al.  A numerical model for the prediction of long term well yield in an unconfined chalk aquifer , 1978, Quarterly Journal of Engineering Geology.

[6]  R. Sakthivadivel,et al.  Numerical analysis of large-diameter wells with a seepage face — Reply , 1990 .

[7]  A. Butler,et al.  Characterising the vertical variations in hydraulic conductivity within the Chalk aquifer , 2006 .

[8]  Lewis Clark,et al.  The analysis and planning of step drawdown tests , 1977, Quarterly Journal of Engineering Geology.

[9]  Z. Şen.  Non-Darcian flow in fractured rocks with a linear flow pattern , 1987 .

[10]  Hongbin Zhan,et al.  Approximate Solutions for Forchheimer Flow to a Well , 2008 .

[11]  M. Mansour,et al.  Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer , 2018, Hydrogeology Journal.

[12]  J. Barker.  A generalized radial flow model for hydraulic tests in fractured rock , 1988 .

[13]  K. Rushton.  Seepage and groundwater flow , 1979 .

[14]  Guanhua Huang,et al.  Non-Darcian flow in a single confined vertical fracture toward a well , 2006 .

[15]  James J. Butler,et al.  Pumping tests in nonuniform aquifers — The radially symmetric case , 1988 .

[16]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[17]  S. Booth,et al.  Pumping-test analysis using a discrete time-discrete space numerical method , 1976 .

[18]  M. Owen,et al.  2. Characteristics and yield in fissured chalk , 1978 .

[19]  Hongbin Zhan,et al.  Experimental evidence of scale-dependent hydraulic conductivity for fully developed turbulent flow in a single fracture , 2007 .

[20]  S. Mathias,et al.  Linearized Richards' equation approach to pumping test analysis in compressible aquifers , 2006 .

[21]  Yi‐Feng Chen,et al.  A generalized Forchheimer radial flow model for constant-rate tests , 2017 .

[22]  Andrew McKenzie,et al.  The physical properties of major aquifers in England and Wales , 1997 .

[23]  J. Barker,et al.  Pumping Tests in Patchy Aquifers , 1982 .

[24]  Z. Şen.  Nonlinear radial flow in confined aquifers toward large‐diameter wells , 1990 .

[25]  Mark Bakker,et al.  Simulation of axi-symmetric flow towards wells: A finite-difference approach , 2012, Comput. Geosci..

[26]  K. Howard,et al.  The Unreliability of Open Observation Boreholes in Unconfined Aquifer Pumping Tests , 1982 .

[27]  H. H. Cooper,et al.  Drawdown in a well of large diameter , 1967 .

[28]  J. L. Amprako.  The United Nations World Water Development Report 2015: Water for a Sustainable World , 2016 .

[29]  C. V. Theis.  The relation between the lowering of the Piezometric surface and the rate and duration of discharge of a well using ground‐water storage , 1935 .

[30]  David Andrew Barry,et al.  MODFLOW equipped with a new method for the accurate simulation of axisymmetric flow. , 2003 .

[31]  K. Rushton,et al.  Pumping Test Analysis When Parameters Vary with Depth , 1976 .

[32]  Christian D Langevin,et al.  Modeling Axisymmetric Flow and Transport , 2008, Ground water.

[33]  Z. Şen.  Non-Darcian groundwater flow in leaky aquifers , 2000 .

[34]  P. Döll.  Vulnerability to the impact of climate change on renewable groundwater resources: a global-scale assessment , 2009 .

[35]  C. E. Jacob,et al.  A generalized graphical method for evaluating formation constants and summarizing well‐field history , 1946 .

[36]  K. Rushton.  Significance of a seepage face on flows to wells in unconfined aquifers , 2006, Quarterly Journal of Engineering Geology and Hydrogeology.

[37]  R. J. Hopkirk,et al.  Observation and simulation of non‐Darcian flow transients in fractured rock , 1997 .

[38]  S. Mathias,et al.  Step‐drawdown tests and the Forchheimer equation , 2010 .

[39]  I. Sobol.  On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[40]  S. Mathias,et al.  Numerical simulation of Forchheimer flow to a partially penetrating well with a mixed-type boundary condition , 2015 .

[41]  N. Brozović,et al.  Effects of initial aquifer conditions on economic benefits from groundwater conservation , 2017 .

[42]  K. R. Rushton,et al.  NUMERICAL METHOD OF PUMPING TEST ANALYSIS USING MICROCOMPUTERS , 1984 .

[43]  A. Gallagher,et al.  Analysis of flow processes in fractured chalk under pumped and ambient conditions (UK) , 2009 .

[44]  James J. Butler,et al.  Pumping tests in nonuniform aquifers: The radially asymmetric case , 1993 .

[45]  Hongbin Zhan,et al.  Numerical modeling of Forchheimer flow to a pumping well in a confined aquifer using the strong-form mesh-free method , 2014, Hydrogeology Journal.

[46]  S. Mathias,et al.  Recovering tracer test input functions from fluid electrical conductivity logging in fractured porous rocks , 2005 .

[47]  M. Mansour,et al.  Pumping test analysis using a layered cylindrical grid numerical model in a complex, heterogeneous Chalk Aquifer , 2011 .

[48]  S. P. Neuman.  Theory of flow in unconfined aquifers considering delayed response of the water table , 1972 .

[49]  Z. Şen.  Nonlinear Flow Toward Wells , 1989 .

[50]  C. E. Jacob.  Drawdown Test to Determine Effective Radius of Artesian Well , 1947 .

[51]  Rory N. Mortimore,et al.  Lithostratigraphy for mapping the Chalk of southern England , 1997 .