Saltwater Upconing and Decay Beneath a Well Pumping Above an Interface Zone

Saltwater, or brine, underlies fresh water in many aquifers, with a transition zone separating them. Pumping fresh water by wells located above the transition zone produces upconing of the latter, eventually salinizing the pumped water, forcing shut-off. The salinity of the pumped water depends on the pumping rate, on the location of the well's screen, on the fresh water flow regime, and on the difference in density between fresh and salt water, expressed as a dimensionless factor called density difference factor (DDF). Following the well’s shut-off, the upconed saltwater mound undergoes decay, tending to return to the prepumping regime. In this paper, the upconing-decay processes in an axially symmetrical system are investigated to discover how they are affected by the DDF and by the dispersivities. The code FEAS-Brine, developed for the simulation of coupled density-dependent flow and salt transport, is used. In this code, the flow equation is solved by the Galerkin finite element method (FEM), while the advective-dispersive salt transport equation is solved in the EulerianLagrangian framework. This code does not suffer from the instability constraint on the Peclet number in the vicinity of the pumping well, where advection dominates the salt transport. Simulation results show that upconing is very sensitive to the DDF, which, in our work, is in the range from 0 (for ideal tracer) to 0.2 (for brine). It is shown that for the DDF of 0.025 (for seawater), local upconing occurs only for low iso-salinity surfaces, while those of high salt concentration, practically, do not shift toward the pumping well. For an ideal tracer, all isosalinity surfaces rise toward the pumping well. For brine, however, only iso-salinity surfaces of very low salinity upcone towards the pumping well. The decay process is lengthy; it takes a long time for the upconed saltwater to migrate back to the original horizontal transition zone prior to pumping. However, the wider transition zone caused by hydrodynamic dispersion can never return to the initial one. This indicates that once a pumping well is abandoned because of high salinity, it can be reused for groundwater utilization only after a long time. Quanlin Zhou et al. (To TIPM) Page 2 7/23/2004

[1]  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 .

[2]  W. C. Krumbein,et al.  The Flow of Homogeneous Fluids through Porous Media , 1938 .

[3]  C. E. Jacob,et al.  Non-steady radial flow in an infinite leaky aquifer , 1955 .

[4]  Louis A. Schmittroth,et al.  Numerical inversion of Laplace transforms , 1960, Commun. ACM.

[5]  M. S. Hantush.  Modification of the theory of leaky aquifers , 1960 .

[6]  M. S. Hantush.  Hydraulics of Wells , 1964 .

[7]  Gedeon Dagan,et al.  Moving Interface in Coastal Aquifers , 1964 .

[8]  M. J. Mundorff,et al.  Electric-analog studies of Brine coning beneath freshwater wells in the Punjab region, West Pakistan , 1968 .

[9]  S. P. Neuman,et al.  Theory of Flow in a Confined Two Aquifer System , 1969 .

[10]  O. Strack.  Some cases of interface flow towards drains , 1972 .

[11]  E. C. Childs.  Dynamics of fluids in Porous Media , 1973 .

[12]  B. M. Sahni.  Physics of Brine Coning Beneath Skimming Wells , 1973 .

[13]  R. Haubold.  Approximation for Steady Interface Beneath a Well Pumping Fresh Water Overlying Salt Water , 1975 .

[14]  G. Pinder,et al.  Approximate analysis of upconing , 1977 .

[15]  R. C. Earlougher.  Advances in Well Test Analysis , 1977 .

[16]  A. N. Stokes,et al.  An Improved Method for Numerical Inversion of Laplace Transforms , 1982 .

[17]  R. Miller,et al.  Update on the High-Current Injector for the Stanford Linear Collider , 1983, IEEE Transactions on Nuclear Science.

[18]  M. Thiele,et al.  Finite-element analysis of dispersion-affected saltwater upconing below a pumping well , 1984 .

[19]  M. A. Ajiz,et al.  A robust incomplete Choleski‐conjugate gradient algorithm , 1984 .

[20]  A. Moench,et al.  Transient Flow to a Large‐Diameter Well in an Aquifer With Storative Semiconfining Layers , 1985 .

[21]  Randall J. Charbeneau,et al.  Saltwater upconing in unconfined aquifers , 1985 .

[22]  Kwang-Je Kim,et al.  Brightness, coherence, and propagation characteristics of synchrotron radiation , 1986 .

[23]  W. R. Souza,et al.  Variable density flow and solute transport simulation of regional aquifers containing a narrow freshwater‐saltwater transition zone , 1987 .

[24]  A. S. Goodman,et al.  Analysis of saltwater upconing beneath a pumping well , 1987 .

[25]  The Impact of the Dupuit‐Forchheimer Approximation on Salt‐Water Intrusion Simulation , 1987 .

[26]  A. Herbert,et al.  Coupled groundwater flow and solute transport with fluid density strongly dependent upon concentration , 1988 .

[27]  T. D. Streltsova,et al.  Well Testing in Heterogeneous Formations , 1988 .

[28]  Kwang‐Je Kim,et al.  RF AND SPACE-CHARGE INDUCED EMITTANCES IN LASER-DRIVEN RF GUNS , 1988 .

[29]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[30]  S. Majid Hassanizadeh,et al.  On the modeling of brine transport in porous media , 1988 .

[31]  Kwang‐Je Kim,et al.  Rf and space-charge effects in laser-driven rf electron guns , 1989 .

[32]  J. Bear,et al.  Introduction to Modeling of Transport Phenomena in Porous Media , 1990 .

[33]  R. Palmer,et al.  Preliminary study of gun emittance correction , 1990 .

[34]  F. J. Kuchuk,et al.  Pressure transients during injection: Constant rate and convolution solutions , 1993 .

[35]  Alexander H.-D. Cheng,et al.  Multilayered leaky aquifer systems: 1. Pumping well solutions , 1993 .

[36]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[37]  R. Buddemeier,et al.  Modeling saltwater upconing in a freshwater aquifer in south-central Kansas , 1997 .

[38]  G. Moridis.  Semianalytical solutions for parameter estimation in diffusion cell experiments , 1998 .

[39]  Vedat Batu.  Aquifer Hydraulics: A Comprehensive Guide to Hydrogeologic Data Analysis , 1998 .

[40]  Gedeon Dagan,et al.  Free‐surface flow toward a well and interface upconing in stratified aquifers of random conductivity , 1998 .

[41]  L. H. Motz,et al.  Analytical solution for determining the critical condition of saltwater upconing in a leaky artesian aquifer , 1999 .

[42]  Multilayered aquifer systems : fundamentals and applications , 2000 .

[43]  Accurate calculation of specific discharge in heterogeneous porous media , 2001 .

[44]  G. Wittum,et al.  The saltpool benchmark problem – numerical simulation of saltwater upconing in a porous medium , 2002 .

[45]  Time-Dependent Measurement of Strongly Density-Dependent Flow in a Porous Medium via Nuclear Magnetic Resonance Imaging , 2002 .

[46]  James J. Butler,et al.  Pumping‐induced leakage in a bounded aquifer: An example of a scale‐invariant phenomenon , 2002 .

[47]  Hongbin Zhan,et al.  A method of calculating pumping induced leakage , 2006 .