Numerical modelling of the impacts of water abstraction for hydraulic fracturing on groundwater–surface water interaction: a case study from northwestern Alberta,Canada

ABSTRACT

The temporal dynamics of groundwater–surface water interaction under the impacts of various water abstraction scenarios are presented for hydraulic fracturing in a shale gas and oil play area (23 984.9 km²), Alberta, Canada, using the MIKE-SHE and MIKE-11 models. Water-use data for hydraulic fracturing were obtained for 433 wells drilled in the study area in 2013 and 2014. Modelling results indicate that water abstraction for hydraulic fracturing has very small (<0.35%) negative impacts on mean monthly and annual river and groundwater levels and stream and groundwater flows in the study area, and small (1–4.17%) negative impacts on environmental flows near the water abstraction location during low-flow periods. The impacts on environmental flow depend on the amount of water abstraction and the daily flow over time at a specific river cross-section. The results also indicate a very small (<0.35%) positive impact on mean monthly and annual groundwater contributions to streamflow because of the large study area. The results provide useful information for planning long-term seasonal and annual water abstractions from the river and groundwater for hydraulic fracturing in a large study area.     

Influence of Hydraulic Fracturing on Overlying Aquifers in the Presence of Leaky Abandoned Wells

The association between hydrocarbon‐rich reservoirs and organic‐rich source rocks means unconventional oil and gas plays usually occur in mature sedimentary basins—where large‐scale conventional development has already taken place. Abandoned wells in proximity to hydraulic fracturing could be affected by increased fluid pressures and corresponding newly generated fractures that directly connect (frac hit) to an abandoned well or to existing fractures intersecting an abandoned well. If contaminants migrate to a pathway hydraulically connected to an abandoned well, upward leakage may occur. Potential effects of hydraulic fracturing on upward flow through a particular type of leaky abandoned well—abandoned oil and gas wells converted into water wells were investigated using numerical modeling. Several factors that affect flow to leaky wells were considered including proximity of a leaky well to hydraulic fracturing, flowback, production, and leaky well abandonment methods. The numerical model used historical records and available industry data for the Eagle Ford Shale play in south Texas. Numerical simulations indicate that upward contaminant migration could occur through leaky converted wells if certain spatial and hydraulic conditions exist. Upward flow through leaky converted wells increased with proximity to hydraulic fracturing, but decreased when flowback and production occurred. Volumetric flow rates ranged between 0 and 0.086 m³/d for hydraulic‐fracturing scenarios. Potential groundwater impacts should be paired with plausible transport mechanisms, and upward flow through leaky abandoned wells could be unrelated to hydraulic fracturing. The results also underscore the need to evaluate historical activities.     

Groundwater Protection and Unconventional Gas Extraction: The Critical Need for Field‐Based Hydrogeological Research

Unconventional natural gas extraction from tight sandstones, shales, and some coal‐beds is typically accomplished by horizontal drilling and hydraulic fracturing that is necessary for economic development of these new hydrocarbon resources. Concerns have been raised regarding the potential for contamination of shallow groundwater by stray gases, formation waters, and fracturing chemicals associated with unconventional gas exploration. A lack of sound scientific hydrogeological field observations and a scarcity of published peer‐reviewed articles on the effects of both conventional and unconventional oil and gas activities on shallow groundwater make it difficult to address these issues. Here, we discuss several case studies related to both conventional and unconventional oil and gas activities illustrating how under some circumstances stray or fugitive gas from deep gas‐rich formations has migrated from the subsurface into shallow aquifers and how it has affected groundwater quality. Examples include impacts of uncemented well annuli in areas of historic drilling operations, effects related to poor cement bonding in both new and old hydrocarbon wells, and ineffective cementing practices. We also summarize studies describing how structural features influence the role of natural and induced fractures as contaminant fluid migration pathways. On the basis of these studies, we identify two areas where field‐focused research is urgently needed to fill current science gaps related to unconventional gas extraction: (1) baseline geochemical mapping (with time series sampling from a sufficient network of groundwater monitoring wells) and (2) field testing of potential mechanisms and pathways by which hydrocarbon gases, reservoir fluids, and fracturing chemicals might potentially invade and contaminate useable groundwater.     

Assessing Methane in Shallow Groundwater in Unconventional Oil and Gas Play Areas,Eastern Kentucky

The expanding use of horizontal drilling and hydraulic fracturing technology to produce oil and gas from tight rock formations has increased public concern about potential impacts on the environment, especially on shallow drinking water aquifers. In eastern Kentucky, horizontal drilling and hydraulic fracturing have been used to develop the Berea Sandstone and the Rogersville Shale. To assess baseline groundwater chemistry and evaluate methane detected in groundwater overlying the Berea and Rogersville plays, we sampled 51 water wells and analyzed the samples for concentrations of major cations and anions, metals, dissolved methane, and other light hydrocarbon gases. In addition, the stable carbon and hydrogen isotopic composition of methane (δ¹³C‐CH₄ and δ²H‐CH₄) was analyzed for samples with methane concentration exceeding 1 mg/L. Our study indicates that methane is a relatively common constituent in shallow groundwater in eastern Kentucky, where methane was detected in 78% of the sampled wells (40 of 51 wells) with 51% of wells (26 of 51 wells) exhibiting methane concentrations above 1 mg/L. The δ¹³C‐CH₄ and δ²H‐CH₄ ranged from ?84.0‰ to ?58.3‰ and from ?246.5‰ to ?146.0‰, respectively. Isotopic analysis indicated that dissolved methane was primarily microbial in origin formed through CO₂ reduction pathway. Results from this study provide a first assessment of methane in the shallow aquifers in the Berea and Rogersville play areas and can be used as a reference to evaluate potential impacts of future horizontal drilling and hydraulic fracturing activities on groundwater quality in the region.     

Organic Contaminants in Portland Cements Used in Monitoring Well Construction

We report the results of two independent laboratory investigations to evaluate total and leachable concentrations of glycols, glycol ethers, phenol, and other compounds in representative Type I and Type I/II Portland cement products that United States Environmental Protection Agency (EPA), The American Society of Testing and Materials (ASTM) and others recommend as annular sealants in monitoring well completions. Water well drillers also use these cements in their well completions. The EPA has included some of these compounds for analysis in their National Hydraulic Fracturing Study to evaluate the effects of hydraulic fracturing on ground‐ and surface water resources. During any contaminant investigation, materials used in monitoring or water well drilling, completion, development, and sampling must be free of the chemicals being targeted by the regulatory agency. Three of five bulk cement products we tested contained part per million (ppm) (mg/kg) concentrations of diethylene glycol, ethylene glycol, tetraethylene glycol, and triethylene glycol; chemicals added as grinding aids during manufacture. Some cements also contained ppb (µg/kg) concentrations of benzoic acid, phenols, propylene glycol, and 2‐butoxyethanol. Leaching of cured cement samples in water produced ppm (mg/L) concentrations of glycols in the supernatant. These results show that cured cements in monitoring or water wells can contaminate groundwater samples with glycols and phenol. Our findings should help prevent future sample bias and false positives when testing for glycol compounds and phenol in groundwater samples from monitoring or water wells and highlight the need to test materials or products used in monitoring or water well drilling, completions, development, and sampling to avoid false positives when sampling and analyzing for less common analytical constituents.     

Numerical Modeling of Methane Leakage from a Faulty Natural Gas Well into Fractured Tight Formations

Horizontal drilling and hydraulic fracturing have enabled hydrocarbon recovery from unconventional reservoirs, but led to natural gas contamination of shallow groundwaters. We describe and apply numerical models of gas‐phase migration associated with leaking natural gas wells. Three leakage scenarios are simulated: (1) high‐pressure natural gas pulse released into a fractured aquifer; (2) continuous slow leakage into a tilted fractured formation; and (3) continuous slow leakage into an unfractured aquifer with fluvial channels, to facilitate a generalized evaluation of natural gas transport from faulty natural gas wells. High‐pressure pulses of gas leakage into sparsely fractured media are needed to produce the extensive and rapid lateral spreading of free gas previously observed in field studies. Transport in fractures explains how methane can travel vastly different distances and directions laterally away from a leaking well, which leads to variable levels of methane contamination in nearby groundwater wells. Lower rates of methane leakage (≤1 Mcf/day) produce shorter length scales of gas transport than determined by the high‐pressure scenario or field studies, unless aquifers have low vertical permeabilities (≤1 millidarcy) and fractures and bedding planes have sufficient tilt (～10°) to allow a lateral buoyancy component. Similarly, in fractured rock aquifers or where permeability is controlled by channelized fluvial deposits, lateral flow is not sufficiently developed to explain fast‐developing gas contamination (0‐3 months) or large length scales (～1 km) documented in field studies. Thus, current efforts to evaluate the frequency, mechanism, and impacts of natural gas leakage from faulty natural gas wells likely underestimate contributions from small‐volume, low‐pressure leakage events.     

An Analytical Model for Hydraulic Fracturing in Shallow Bedrock Formations

A theoretical method is proposed to estimate post‐fracturing fracture size and transmissivity, and as a test of the methodology, data collected from two wells were used for verification. This method can be employed before hydrofracturing in order to obtain estimates of the potential hydraulic benefits of hydraulic fracturing. Five different pumping test analysis methods were used to evaluate the well hydraulic data. The most effective methods were the Papadopulos‐Cooper model (1967), which includes wellbore storage effects, and the Gringarten‐Ramey model (1974), known as the single horizontal fracture model. The hydraulic parameters resulting from fitting these models to the field data revealed that as a result of hydraulic fracturing, the transmissivity increased more than 46 times in one well and increased 285 times in the other well. The model developed by dos Santos (2008) , which considers horizontal radial fracture propagation from the hydraulically fractured well, was used to estimate potential fracture geometry after hydrofracturing. For the two studied wells, their fractures could have propagated to distances of almost 175 m or more and developed maximum apertures of about 2.20 mm and hydraulic apertures close to 0.30 mm. Fracturing at this site appears to have expanded and propagated existing fractures and not created new fractures. Hydraulic apertures calculated from pumping test analyses closely matched the results obtained from the hydraulic fracturing model. As a result of this model, post‐fracturing geometry and resulting post‐fracturing well yield can be estimated before the actual hydrofracturing.     

Conceptualization of the hydrogeological system of some sedimentary aquifers in Savelugu–Nanton and surrounding areas,Northern Ghana

Groundwater resources play a pivotal role in the rural water delivery system in Ghana. The hydrogeological system of Middle Voltaian terrain was simulated using available data on hydraulic heads and boundary conditions. The objective was to characterize the general groundwater flow pattern and provide local estimates of the distribution of hydraulic conductivity and recharge fields. The results suggest a predominant NE–SW flow direction, which ties in with the general regional structural trend and indicates that the hydrogeological conditions of the rocks are controlled by structural entities created in the wake of fracturing and/or weathering of the rocks whose primary permeabilities are considerably reduced because of high compaction and low‐grade metamorphism. Calibrated hydraulic conductivities range between 1.90 and 10.81 m/d. The spatial distribution appears to reflect the intensity of fracturing and/or weathering of the rock and the proportion of the clay fraction of the weathered zone. Vertical groundwater recharge has been estimated to range between 0.3% and 4.1% of the annual rainfall. This recharge rate is quite low and reflects the imperviousness of the thick overburden because of high clay content in some places and high compaction in others. Despite this apparently low recharge rate, groundwater resources potential in the area appear to be high, and increased abstraction from existing abstraction wells by up to 50% does not appear to register significant effects on groundwater budgets at the simulated recharge rates. This suggests that the well yields are much lower than the potential of the aquifer system. The apparently low yields might be associated with poor well development and the choice of inappropriate well completion materials. This study recommends a monitoring system to be developed for a much more regional groundwater flow simulation under transient conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Trace Element Behavior in Methane‐Rich and Methane‐Free Groundwater in North and East Texas

There is concern about adverse impacts of natural gas (primarily methane) production on groundwater quality; however, data on trace element concentrations are limited. The objective of this study was to compare the distribution of trace elements in groundwater samples with and without dissolved methane in aquifers overlying the Barnett Shale (Hood and Parker counties, 207 samples) and the Haynesville Shale (Panola County, 42 samples). Both shales have been subjected to intensive hydraulic fracturing for gas production. Well clusters with high dissolved methane were previously found in these counties and are thought to be of natural origin. Overall, groundwater in these counties is of excellent quality with typically low elemental concentrations. Several statistical analyses strongly suggest that most trace element concentrations, generally at low background levels, are no higher and even reduced when dissolved methane is present. In addition, trace element concentrations are not correlated with distance to gas wells. The reduction in trace element concentrations is attributed to anaerobic microbial degradation of methane, is associated with a higher pH (>8.5), and, likely, with precipitation of carbonates and pyrite and formation of clays. Trace and other elements are likely incorporated within the precipitating mineral crystalline network or sorbed. High pH values are found throughout these high‐methane clusters (e.g., Parker‐Hood cluster), even in subregions where methane is not present, which is consistent with a pervasive natural origin of dissolved methane rather than a limited gas well source.     

Temporal scaling of groundwater level fluctuations near a stream

Temporal scaling in stream discharge and hydraulic heads in riparian wells was evaluated to determine the feasibility of using spectral analysis to identify potential surface and groundwater interaction. In floodplains where groundwater levels respond rapidly to precipitation recharge, potential interaction is established if the hydraulic head (h) spectrum of riparian groundwater has a power spectral density similar to stream discharge (Q), exhibiting a characteristic breakpoint between high and low frequencies. At a field site in Walnut Creek watershed in central Iowa, spectral analysis of h in wells located 1 m from the channel edge showed a breakpoint in scaling very similar to the spectrum of Q (～20 h), whereas h in wells located 20 and 40 m from the channel showed temporal scaling from 1 to 10,000 h without a well-defined breakpoint. The spectral exponent (β) in the riparian zone decreased systematically from the channel into the floodplain as groundwater levels were increasingly dominated by white noise groundwater recharge. The scaling pattern of hydraulic head was not affected by land cover type, although the number of analyses was limited and site conditions were variable among sites. Spectral analysis would not replace quantitative tracer or modeling studies, but the method may provide a simple means of confirming potential interaction at some sites.     

Field verification of stable perched groundwater in layered bedrock uplands

Data substantiating perched conditions in layered bedrock uplands are rare and have not been widely reported. Field observations in layered sedimentary bedrock in southwestern Wisconsin, USA, provide evidence of a stable, laterally extensive perched aquifer. Data from a densely instrumented field site show a perched aquifer in shallow dolomite, underlain by a shale-and-dolomite aquitard approximately 25 m thick, which is in turn underlain by sandstone containing a 30-m-thick unsaturated zone above a regional aquifer. Heads in water supply wells indicate that perched conditions extend at least several kilometers into hillsides, which is consistent with published modeling studies. Observations of unsaturated conditions in the sandstone over a 4-year period, historical development of the perched aquifer, and perennial flow from upland springs emanating from the shallow dolomite suggest that perched groundwater is a stable hydrogeologic feature under current climate conditions. Water-table hydrographs exhibit apparent differences in the amount and timing of recharge to the perched and regional flow systems; steep hydraulic gradients and tritium and chloride concentrations suggest there is limited hydraulic connection between the two. Recognition and characterization of perched flow systems have practical importance because their groundwater flow and transport pathways may differ significantly from those in underlying flow systems. Construction of multi-aquifer wells and groundwater withdrawal in perched systems can further alter such pathways.     

Tube wave to shear wave conversion at borehole plugs

In hydraulic fracturing experiments, perforation shots excite body and tube waves that sample, and thus can be used to characterize, the surrounding medium. While these waves are routinely employed in borehole operations, their resolving power is limited by the experiment geometry, the signal‐to‐noise ratio, and their frequency content. It is therefore useful to look for additional, complementary signals that could increase this resolving power. Tube‐to‐body‐wave conversions (scattering of tube to compressional or shear waves at borehole discontinuities) are one such signal. These waves are not frequently considered in hydraulic fracture settings, yet they possess geometrical and spectral attributes that greatly complement the resolution afforded by body and tube waves alone. Here, we analyze data from the Jonah gas field (Wyoming, USA) to demonstrate that tube‐to‐shear‐wave conversions can be clearly observed in the context of hydraulic fracturing experiments. These waves are identified primarily on the vertical and radial components of geophones installed in monitoring wells surrounding a treatment well. They exhibit a significantly lower frequency content (10–100 Hz) than the primary compressional waves (100–1000 Hz). Tapping into such lower frequencies could help to better constrain velocity in the formation, thus allowing better estimates of fracture density, porosity and permeability. Moreover, the signals of tube‐to‐shear‐wave conversion observed in this particular study provide independent estimates of the shear wave velocity in the formation and of the tube wave velocity in the treatment well.     

Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Commingled Fluids in Abandoned Boreholes: Proximity Analysis of a Hidden Liability

The interactions between old abandoned wellbores of suspect well integrity with hydraulic fracturing (HF), enhanced oil recovery (EOR), or salt water disposal (SWD) operations can result in upward leakage of deep aqueous liquids into overlying aquifers. This potential for upward fluid migration is largely unquantified as monitoring abandoned wells is rarely done, and leakage may go unnoticed especially when in deeper aquifers. This study performs a proximity analysis between old abandoned wells and HF, EOR, and SWD wells, and identifies commingled old abandoned wellbores, which are those wells where groundwater may flow from one aquifer to one or more other aquifers, to identify the locations with the greatest potential for upward aqueous fluid migration at three study sites in the Western Canadian Sedimentary Basin. Our analysis indicates that at all three study sites there are several locations where HF, EOR, or SWD operations are located in close proximity to a given old abandoned well. Much of this overlap occurs in formations above typically produced hydrocarbon reservoirs but below exploited potable aquifers, otherwise known as the intermediate zone, which is often connected between abandonment plugs in old abandoned wells. Information on the intermediate zone is often lacking, and this study suggests that unanticipated alterations to groundwater flow systems within the intermediate zone may be occurring. Results indicate the need for more field-based research on the intermediate zone.     

Modelling study on the impact of deep building foundations on the groundwater system

Coastal areas are usually the preferred place of habitation for human beings. Anthropogenic activities such as the construction of high‐rise buildings and underground transport systems usually require extensive deep foundations and ground engineering works, which may unintentionally modify the coastal groundwater system because the construction materials of foundations are usually of low hydraulic conductivity. In this paper, the impact of these building foundations on the groundwater regime is studied using hypothetical flow and transport models. Various possible realizations of foundation distributions are generated using stochastic parameters derived from a topographical map of an actual coastal area in Hong Kong. The effective hydraulic conductivity is first calculated for different realizations and the results show that the effective hydraulic conductivity can be reduced significantly. Then a hypothetical numerical model based on FEFLOW is set up to study the change of hydraulic head, groundwater discharge, and saltwater‐fresh water interface. The groundwater level and flow are modified to various degrees, depending on the foundations percentage and the distribution pattern of the buildings. When the foundations percentage is high and the building foundations are aggregated, the hydraulic head is raised significantly and the originally one‐dimensional groundwater flow field becomes complicated. Seaward groundwater discharge will be reduced and some groundwater may become seepage through the ground surface. The transport model shows that, after foundations are added, overall the seawater and fresh groundwater interface moves landward, so extensive foundations may induce seawater intrusion. It is believed that the modification of the coastal groundwater system by building foundations may have engineering and environmental implications, such as submarine groundwater discharge, foundation corrosion, and slope stability. Copyright © 2007 John Wiley & Sons, Ltd.     

Conventional Oil—The Forgotten Part of the Water-Energy Nexus

The impacts of unconventional oil and gas production via high-volume hydraulic fracturing (HVHF) on water resources, such as water use, groundwater and surface water contamination, and disposal of produced waters, have received a great deal of attention over the past decade. Conventional oil and gas production (e.g., enhanced oil recovery [EOR]), which has been occurring for more than a century in some areas of North America, shares the same environmental concerns, but has received comparatively little attention. Here, we compare the amount of produced water versus saltwater disposal (SWD) and injection for EOR in several prolific hydrocarbon producing regions in the United States and Canada. The total volume of saline and fresh to brackish water injected into depleted oil fields and nonproductive formations is greater than the total volume of produced waters in most regions. The addition of fresh to brackish “makeup” water for EOR may account for the net gain of subsurface water. The total amount of water injected and produced for conventional oil and gas production is greater than that associated with HVHF and unconventional oil and gas production by well over a factor of 10. Reservoir pressure increases from EOR and SWD wells are low compared to injection of fluids for HVHF, however, the longer duration of injections could allow for greater solute transport distances and potential for contamination. Attention should be refocused from the subsurface environmental impacts of HVHF to the oil and gas industry as a whole.     

Isotope Approach to Assess Hydrologic Connections During Marcellus Shale Drilling

Water and gas samples were collected from (1) nine shallow groundwater aquifers overlying Marcellus Shale in north‐central West Virginia before active shale gas drilling, (2) wells producing gas from Upper Devonian sands and Middle Devonian Marcellus Shale in southwestern Pennsylvania, (3) coal‐mine water discharges in southwestern Pennsylvania, and (4) streams in southwestern Pennsylvania and north‐central West Virginia. Our preliminary results demonstrate that the oxygen and hydrogen isotope composition of water, carbon isotope composition of dissolved inorganic carbon, and carbon and hydrogen isotope compositions of methane in Upper Devonian sands and Marcellus Shale are very different compared with shallow groundwater aquifers, coal‐mine waters, and stream waters of the region. Therefore, spatiotemporal stable isotope monitoring of the different sources of water before, during, and after hydraulic fracturing can be used to identify migrations of fluids and gas from deep formations that are coincident with shale gas drilling.     

Three-dimensional geostatistical inversion of synthetic tomographic pumping and heat-tracer tests in a nested-cell setup

A main purpose of groundwater inverse modeling lies in estimating the hydraulic conductivity field of an aquifer. Traditionally, hydraulic head measurements, possibly obtained in tomographic setups, are used as data. Because the groundwater flow equation is diffusive, many pumping and observation wells would be necessary to obtain a high resolution of hydraulic conductivity, which is typically not possible. We suggest performing heat tracer tests using the same already installed pumping wells and thermometers in observation planes to amend the hydraulic head data set by the arrival times of the heat signals. For each tomographic combinations of wells, we recommend installing an outer pair of pumping wells, generating artificial ambient flow, and an inner well pair in which the tests are performed. We jointly invert heads and thermal arrival times in 3-D by the quasi-linear geostatistical approach using an efficiently parallelized code running on a mid-range cluster. In the present study, we evaluate the value of heat tracer versus head data in a synthetic test case, where the estimated fields can be compared to the synthetic truth. Because the sensitivity patterns of the thermal arrival times differ from those of head measurements, the resolved variance in the estimated field is 6 to 10 times higher in the joint inversion in comparison to inverting head data only. Also, in contrast to head measurements, reversing the flow field and repeating the heat-tracer test improves the estimate in terms of reducing the estimation variance of the estimate. Based on the synthetic test case, we recommend performing the tests in four principal directions, requiring in total eight pumping wells and four intersecting observation planes for heads and temperature in each direction.     

识别地下水超采区井水位异常性质的理论与方法

长期过量开采地下水,使地下水位持续下降、水质发生变化,动水位观测井断流;地面沉降造成井管上窜,观测管路系统被损坏等,这些现象对地震地下流体观测地震前兆异常的正确判断带来很大困难。应用水文地质理论与方法,分析含水层的水均衡状态、应力-应变状态及其与水位动态的关系,初步探讨了超采区井水位异常性质的理论与方法。结果表明,根据井孔所在区水位下降漏斗的扩散特征,结合以上所提到的理论和方法,依据资料多年变化特征,可以较准确地判断异常的性质。研究结果有助于区分单一集中抽水与长期地下水超采对水位观测的影响,有助于正确识别超采区水位前兆异常,有助于地震分析预报水平的提高