Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Study of subsurface geology in locating arsenic-free groundwater in Bengal delta,West Bengal,India

Over a large area of the Bengal delta in West Bengal, India, arsenic distribution patterns in groundwater were studied. One hundred and ten boreholes at different target locations were made, subsurface sediments were logged and analysed, and arsenic values in sediments vis-à-vis groundwater were compared. The study elucidates the subsurface geology of the western part of Bengal delta and characterises the sediments that were intersected in different boreholes with contrasting values of arsenic in groundwater. It reveals an existence of multiple aquifers stacked over each other. Depending on the color and nature of aquifer-sands and their overlying clay beds six aquifer types (Type-1 to Type-6) are classified and described. Sediment-arsenic for all the varieties of aquifer sands are near similar but the groundwater-arsenic of these six aquifers varies widely. Type-2 and Type-5 aquifers host arsenic-contaminated groundwater whereas the other four aquifers are arsenic-free. Type-2 and Type-5 aquifers are capped by a grey to dark grey soft organic matter-rich clay unit which makes these aquifers semi-confined to leaky-confined. These contribute in releasing arsenic from the sediments. The results of this study are employed in a proposed georemedial measure against this hazardous toxic element.     

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

Simulation of groundwater level recovery in abandoned mines,Fengfeng coalfield,China

Abandoned mines are of high potential risk as they could be a large underground storage of pollutants(heavy metals and organic wastes, etc.). Various physical, chemical and biological reactions would take place when groundwater flows into underground spaces, which makes abandoned mine a huge potential hazard to groundwater environment. The recovery of groundwater level is one of the key elements controlling the reactions and causing such hazards. This paper simulated groundwater level recovery processes in the abandoned mines, Fengfeng coalfield by using the computer program FEFLOW. The paper integrated the pipe flow model, "three zones" model and groundwater inrush(discharge) model in the simulation of groundwater in the complex laneway-aquifer system. Groundwater flow in the laneway systems was considered pipe flow and described in Bernoulli equation. The water-bearing medium(coal seam roof) overlying the laneway systems was divided into "three zones" composed of the caving zone, fissure zone and bending zone based on the disruption degrees of previous mining. Groundwater in the Ordovician limestone aquifer(bottom of coal seam) flowing into laneway systems was considered a major inrush/recharge source, and its flow rate was calculated by an inrush(discharge) model which was newly developed in this study and incorporated into FEFLOW. The results showed that it would take approximately 95 days for groundwater in abandoned mines to recover to regional groundwater level elevation, and the total amount of water filling up would be about 1.41195×10~7 m~3, which is consistent with the actual data. The study could be of theoretical and practical significance to mitigate abandoned mines' hazards and improve mine groundwater utilization.     

Inverse problem in hydrogeology

The state of the groundwater inverse problem is synthesized. Emphasis is placed on aquifer characterization, where modelers have to deal with conceptual model uncertainty (notably spatial and temporal variability), scale dependence, many types of unknown parameters (transmissivity, recharge, boundary conditions, etc.), nonlinearity, and often low sensitivity of state variables (typically heads and concentrations) to aquifer properties. Because of these difficulties, calibration cannot be separated from the modeling process, as it is sometimes done in other fields. Instead, it should be viewed as one step in the process of understanding aquifer behavior. In fact, it is shown that actual parameter estimation methods do not differ from each other in the essence, though they may differ in the computational details. It is argued that there is ample room for improvement in groundwater inversion: development of user-friendly codes, accommodation of variability through geostatistics, incorporation of geological information and different types of data (temperature, occurrence and concentration of isotopes, age, etc.), proper accounting of uncertainty, etc. Despite this, even with existing codes, automatic calibration facilitates enormously the task of modeling. Therefore, it is contended that its use should become standard practice.

---

Resumen Se sintetiza el estado del problema inverso en aguas subterráneas. El énfasis se ubica en la caracterización de acuíferos, donde los modeladores tienen que enfrentar la incertidumbre del modelo conceptual (principalmente variabilidad temporal y espacial), dependencia de escala, muchos tipos de parámetros desconocidos (transmisividad, recarga, condiciones limitantes, etc), no linealidad, y frecuentemente baja sensibilidad de variables de estado (típicamente presiones y concentraciones) a las propiedades del acuífero. Debido a estas dificultades, no puede separarse la calibración de los procesos de modelado, como frecuentemente se hace en otros campos. En su lugar, debe de visualizarse como un paso en el proceso de entendimiento del comportamiento del acuífero. En realidad, se muestra que los métodos reales de estimación de parámetros no difieren uno del otro en lo esencial, aunque sí pueden diferir en los detalles computacionales. Se discute que existe amplio espacio para la mejora del problema inverso en aguas subterráneas: desarrollo de códigos amigables al usuario, acomodamiento de variabilidad a través de geoestadística, incorporación de información geológica y diferentes tipos de datos (temperatura, presencia y concentración de isótopos, edad, etc), explicación apropiada de incertidumbre, etc. A pesar de esto, aún con los códigos existentes, la calibración automática facilita enormemente la tarea de modelado. Por lo tanto, se sostiene que su uso debería de convertirse en práctica standard.

---

Résumé Létat du problème inverse des eaux souterraines est synthétisé. Laccent est placé sur la caractérisation de laquifère, où les modélisateurs doivent jouer avec lincertitude des modèles conceptuels (notamment la variabilité spatiale et temporelle), les facteurs déchelle, plusieurs inconnues sur différents paramètres (transmissivité, recharge, conditions aux limites, etc.), la non linéarité, et souvent la sensibilité de plusieurs variables détat (charges hydrauliques, concentrations) des propriétés de laquifère. A cause de ces difficultés, le calibrage ne peut être séparé du processus de modélisation, comme cest le cas dans dautres cas de figure. Par ailleurs, il peut être vu comme une des étapes dans le processus de détermination du comportement de laquifère. Il est montré que les méthodes dévaluation des paramètres actuels ne diffèrent pas si ce nest dans les détails des calculs informatiques. Il est montré quil existe une large panoplie de techniques d ‹inversion : codes de calcul utilisables par tout-un-chacun, accommodation de la variabilité via la géostatistique, incorporation dinformations géologiques et de différents types de données (température, occurrence, concentration en isotopes, âge, etc.), détermination de lincertitude. Vu ces développements, la calibration automatique facilite énormément la modélisation. Par ailleurs, il est souhaitable que son utilisation devienne une pratique standardisée.

     

Development of groundwater modeling for the Azraq Basin, Jordan

 The three-dimensional groundwater flow model MODFLOW was applied to simulate water level change in the complex multi-aquifer systems (the Upper and Middle Aquifers) of the Azraq basin. The model was calibrated by matching observed and simulated drawdown for steady and transient states over the period 1970–1992. Drawdown data for the period 1993–1997 were used to test the model's ability to predict the response of the aquifers. The model performed well in representing the water level contours of the Upper and Middle Aquifers for steady state calibration. Agreement between the observed and simulated drawdowns was obtained for transient state calibration. To predict the aquifer system responses for the period of 1997–2025, four different pumping schemes (scenarios) have been investigated. The first scenario (present pumping rate) reveals that there will be approximately a 25 m drop in the water level at the well-field area in 2025. However, the worst scenario (pumping rate at 1.5 times the present rate) reveals an approximate 39 m drop in the water level at the well-field area in 2025. The safe yield for the Upper Aquifer System was found to be about 25 million cubic meters (MCM) yearly. Received: 24 June 1999 · Accepted: 30 November 1999     

Quantitative analysis of springs

Growing demand for groundwater resources and stringent environmental concerns has led to large groundwater investigations, including characterization of aquifer systems that are hydraulically connected to springs. A pumping test is one of the most reliable means of quantifying hydraulic characteristics and the response of natural springs to pumping because it yields results that, in general, are representative of a larger area than are results from a single point observation. Recharge to the aquifer sustaining discharge from springs must be evaluated prior to the utilization of springs. The spring hydrograph is analyzed, as the shape of a hydrograph is a reflection of the response of the aquifer to recharge. The form and rate of recession provide significant information on the storage, lithological composition, and structural characteristics of the aquifer system sustaining the spring. Water tracing techniques have been developed and used over a period of centuries to delineate catchment boundaries, estimate groundwater flow velocities, determine areas of recharge, and identify sources of pollution of spring water.     

Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion

Nonlocal moment equations allow one to render deterministically optimum predictions of flow in randomly heterogeneous media and to assess predictive uncertainty conditional on measured values of medium properties. We present a geostatistical inverse algorithm for steady-state flow that makes it possible to further condition such predictions and assessments on measured values of hydraulic head (and/or flux). Our algorithm is based on recursive finite-element approximations of exact first and second conditional moment equations. Hydraulic conductivity is parameterized via universal kriging based on unknown values at pilot points and (optionally) measured values at other discrete locations. Optimum unbiased inverse estimates of natural log hydraulic conductivity, head and flux are obtained by minimizing a residual criterion using the Levenberg-Marquardt algorithm. We illustrate the method for superimposed mean uniform and convergent flows in a bounded two-dimensional domain. Our examples illustrate how conductivity and head data act separately or jointly to reduce parameter estimation errors and model predictive uncertainty.This work is supported in part by NSF/ITR Grant EAR-0110289. The first author was additionally supported by scholarships from CONACYT and Instituto de Investigaciones Electricas of Mexico. Additional support was provided by the European Commission under Contract EVK1-CT-1999-00041 (W-SAHaRA-Stochastic Analysis of Well Head Protection and Risk Assessment).     

Hydrogeology of limited aquifer in a granitic terrain

 Detailed hydrogeological studies in a granitic micro-watershed have been carried out to determine the extent, behavior, and characteristics of the aquifer. The study includes analysis of lithologs, drill time log, pumping tests, and slug tests. Realistic field conditions have been taken into account for characterizing the aquifer system. Slug tests were carried out to estimate aquifer parameters at the wells which could not sustain pumping. Received: 20 November 1997 · Accepted: 23 February 1998