Predicting Water Resource Impacts of Unconventional Gas Using Simple Analytical Equations.

The rapid expansion in unconventional gas development over the past two decades has led to concerns over the potential impacts on groundwater resources. Although numerical models are invaluable for assessing likelihood of impacts at particular sites, simpler analytical models are also useful because they help develop hydrological understanding. Analytical approaches are also valuable for preliminary assessments and to determine where more complex models are warranted. In this article, we present simple analytical solutions that can be used to predict: (1) the spatial extent of drawdown from horizontal wells drilled into the gas-bearing formation, and rate of recovery after gas production ceases; (2) the potential for upward transport of contaminants from the gas-bearing formation to shallow aquifers during hydraulic fracturing operations when pressures in the gas-bearing formation are greatly increased; and (3) the potential downward leakage of water from shallow aquifers during depressurization of gas-bearing formations. In particular, we show that the recovery of pressure after production ceases from gas-bearing shale formations may take several hundred years, and we present critical hydraulic conductivity values for intervening aquitards, below which the impact on shallow aquifers will be negligible. The simplifying assumptions inherent in these solutions will limit their predictive accuracy for site-specific assessments, compared to numerical models that incorporate knowledge of spatial variations in formation properties and which may include processes not considered in the simpler solutions.

Variable-density groundwater flow and solute transport in heterogeneous porous media: approaches, resolutions and future challenges.

In certain hydrogeological situations, fluid density variations occur because of changes in the solute or colloidal concentration, temperature, and pressure of the groundwater. These include seawater intrusion, high-level radioactive waste disposal, groundwater contamination, and geothermal energy production. When the density of the invading fluid is greater than that of the ambient one, density-driven free convection can lead to transport of heat and solutes over larger spatial scales and significantly shorter time scales than compared with diffusion alone. Beginning with the work of Lord Rayleigh in 1916, thermal and solute instabilities in homogeneous media have been studied in detail for almost a century. Recently, these theoretical and experimental studies have been applied in the study of groundwater phenomena, where the assumptions of homogeneity and isotropy rarely, if ever, apply. The critical role that heterogeneity plays in the onset as well as the growth and/or decay of convective motion is discussed by way of a review of pertinent literature and numerical simulations performed using a variable-density flow and solute transport numerical code. Different styles of heterogeneity are considered and range from continuously "trending" heterogeneity (sinusoidal and stochastic permeability distributions) to discretely fractured geologic media. Results indicate that both the onset of instabilities and their subsequent growth and decay are intimately related to the structure and variance of the permeability field. While disordered heterogeneity tends to dissipate convection through dispersive mixing, an ordered heterogeneity (e.g., sets of vertical fractures) allows instabilities to propagate at modest combinations of fracture aperture and separation distances. Despite a clearer understanding of the processes that control the onset and propagation of instabilities, resultant plume patterns and their migration rates and pathways do not appear amenable to prediction at present. The classical Rayleigh number used to predict the occurrence of instabilities fails, in most cases, when heterogeneous conditions prevail. The incorporation of key characteristics of the heterogeneous permeability field into relevant stability criteria and numerical models remains a challenge for future research.