Combining steam injection with hydraulic fracturing for the in situ remediation of the unsaturated zone of a fractured soil polluted by jet fuel

A steam injection pilot-scale experiment was performed on the unsaturated zone of a strongly heterogeneous fractured soil contaminated by jet fuel. Before the treatment, the soil was stimulated by creating sub-horizontal sand-filled hydraulic fractures at three depths. The steam was injected through one hydraulic fracture and gas/water/non-aqueous phase liquid (NAPL) was extracted from the remaining fractures by applying a vacuum to extraction wells. The injection strategy was designed to maximize the heat delivery over the entire cell (10 m × 10 m × 5 m). The soil temperature profile, the recovered NAPL, the extracted water, and the concentrations of volatile organic compounds (VOCs) in the gas phase were monitored during the field test. GC-MS chemical analyses of pre- and post-treatment soil samples allowed for the quantitative assessment of the remediation efficiency. The growth of the heat front followed the configuration of hydraulic fractures. The average concentration of total hydrocarbons (g/kg of soil) was reduced by ～ 43% in the upper target zone (depth = 1.5-3.9 m) and by ～ 72% over the entire zone (depth = 1.5-5.5 m). The total NAPL mass removal based on gas and liquid stream measurements and the free-NAPL product were almost 30% and 2%, respectively, of those estimated from chemical analyses of pre- and post-treatment soil samples. The dominant mechanisms of soil remediation was the vaporization of jet fuel compounds at temperatures lower than their normal boiling points (steam distillation) enhanced by the ventilation of porous matrix due to the forced convective flow of air. In addition, the significant reduction of the NAPL mass in the less-heated deeper zone may be attributed to the counter-current imbibition of condensed water from natural fractures into the porous matrix and the gravity drainage associated with seasonal fluctuations of the water table.     

WAVE INDUCED FLOW AND TRANSPORT IN SEDIMENT BEDS1

ABSTRACT: The fate of contaminants in large water bodies is highly influenced by the transfer of flow and solutes across the water sediment interface. In this paper, an analytical model is presented where flow in both sediment bed and open channel is coupled at the interface through a boundary layer occupying the upper part of the sediment bed. The presence of this layer allows not only the capture of the inertia effects through a drag term in the generalized Darcy's equation, but also the specification of different soil parameters for the two porous zones. The flow is advective and driven by wave action along the water surface. The resulting system is solved for the pressure and flux in each sediment layer. The generated transport velocity fields are linked to a random walk simulation that is used to examine the trajectories of solute particles. Comparison of these trajectories against experimental tracer tests suggests a pattern very similar to the one attributed to the presence of surface mounds. The results clearly show the significance of the boundary layer and the drag term for soil with high permeability and the impact of both the thickness of the boundary layer and the length of the gravity wave relative to the depth of the water channel on the transport and exchange across the interface. The paper also examines the sensitivity of the mass exchange to the permeability of the two porous zones.     

Nitrate removal in riparian wetlands: interactions between surface flow and soils

Riparian wetlands containing springs are thought to be ineffective at removing nitrate because contact times between the upwelled ground water and the underlying microbially active soils are short. Tracer experiments using lithium bromide (LiBr) and nitrate (NO3-N) injected at the surface were used to quantify residence times and NO3-N removal in a riparian swale characteristic of New Zealand hill-country pasture. An experimental enclosure was used with collecting trays at the downstream end to measure flow and concentration, shallow wells to measure subsurface concentrations, and an array of logging conductivity probes to monitor tracer continuously. The majority of added tracer reached the outlet more slowly than could be explained by surface flow, but more quickly than could be explained by Darcy seepage flow. There was evidence from the wells of tracer diffusing vertically to a depth of at least 5 cm into the surface soil layer, which was permanently saturated and highly porous. During dry weather 24 +/- 9% of added NO3-N was removed over a distance of 1.5 m largely by denitrification. The net uptake length coefficient for this wetland (K = 0.08 +/- 0.03 m(-1)) is slightly higher than the range (K = 0.01-0.07 m(-1)) measured in a small stream channel infested with macrophytes. Nitrate removal is expected to decrease with increasing flow. Seepage flow is estimated to have removed only 7 +/- 4% of the added NO3-N and we hypothesize that vertical diffusion substantially increases NO3-N removal in this type of wetland. Riparian wetlands with springs and surface flows should not be dismissed as having low NO3-N removal potential without checking whether there is significant vertical mixing.     

治理套外返污水的油井、水井

由于固井质量不合格或是套管破裂、错断引起的油井和水井套返污水，造成了吉林油田老油区普遍的环境污染。从１９９０年开始，通过反复的摸索和试验，提出并逐渐完善了两种切实可行的治理方法，即封堵剂法和封隔器法。实践证明，效果良好。     

Flow and diffusion measurements in natural porous media using magnetic resonance imaging

Flow and diffusion of water in natural porous media, quartz sand, and calcareous gravel were measured using a 1.5-T clinical magnetic resonance tomograph. The spatial resolution of the dynamic measurements was 1.32 x 1.32 x 5 mm3, and the time between two cross-sectional measurements was approximately 10 s. The measured coefficients of molecular diffusion for water were in good agreement with theoretical data. Flow was measured without any tracer at velocities between 0.15 and 6.67 mm/s. The results, based on a calibration within one part of the column, were in good agreement with data obtained from a tracer experiment and from a numerical model. It was possible to measure the flow velocity in larger pores and preferential flow paths directly. The results of the flow measurements in smaller pores reflected the mean velocity within that volume element. In that case the obtained values were close to the average linear velocity. Since the time resolution is high a monitoring of flow processes is possible. The pore space was imaged with a spatial resolution of 0.5 x 0.5 x 0.5 mm3. Here, the porosity of pores that are larger than 0.2 mm can be measured directly; for smaller pores a calibration is necessary.     

SPATIAL MODELLING WITH GEOGRAPHIC INFORMATION SYSTEMS FOR DETERMINATION OF WATER RESOURCES VULNERABILITY APPLICATION TO AN AREA IN MASSIF CENTRAL (FRANCE)1

ABSTRACT: A method for water resources protection based on spatial variability of vulnerability is proposed. The vulnerability of a water resource is defined as the risk that the resource will become contaminated if a pollutant is placed on the surface at one point as compared to another. A spatial modelling method is defined in this paper to estimate a travel time between any point of a catchment and a resource (river or well). This method is based on spatial analysis tools integrated in Geographical Information Systems (GIS). The method is illustrated by an application to an area of Massif Central (France) where three different types of flow appear: surface flow, shallow subsurface flow, and permanent ground water flow (baseflow). The proposed method gives results similar to classical methods of estimation of travel time. The contribution of GIS is to improve the mapping of vulnerability by taking the spatial variability of physical phenomena into account.     

Three-dimensional nickel ion transport through porous media using magnetic resonance imaging

The transport of Ni2+ ions in a column, filled with porous media, was observed in three dimensions and time by magnetic resonance imaging (MRI) in a clinical scanner. For porous media we used glass beads or quartz sand in a saturated continuous flow mode. The magnetic moment of Ni2+ decreased the T1 relaxation time of 1H in aqueous solution. This concentration-dependent effect was used by a fast low angle shot (FLASH) MRI sequence for imaging the concentration of the dissolved ions. Since Ni2+ behaves as a conservative tracer under the chosen conditions, the tracer motion was representative for the water flow in the porous medium. Currently, we can achieve an isotropic spatial resolution of 1.5 mm and a temporal resolution of 170 s. The transport observation gives direct access to hydraulic flow properties of the porous media. The fluid flow velocity field was calculated by a fronttracking method and the statistical properties of the velocities were investigated. We also compared the experimental data with the three-dimensional particle tracking model PARTRACE, which uses the experimental flow field as input.     

PERCOLATE WATER AND BROMIDE MOVEMENT IN THE ROOT ZONE OF EFFLUENT IRRIGATION SITES1

ABSTRACT: A bromide tracer was used to evaluate percolate water and ion movement in the upper 1.2 m of soil at a proposed sewage effluent irrigation site located in the Missouri Ozarks. Two plots representing Doniphan silt loam and Crider silt loam soils were sprinkler irrigated with local ground water at a rate of 7.62 cm/week from June through August 1976. Soil water potential, percent soil moisture by volume, and background levels of bromide in soil water, ground water, and precipitation were measured at the study plots. Bromide exchange properties and saturated hydraulic conductivity of the soils were determined in the laboratory. During two selected time periods, irrigation water, was spiked with NaBr (5.0 mg/l Br). Bromide movement through the upper profile was quantified by soil water samples and post-sampling neutron activation analysis. Soil moisture was near saturatin in both soils when the Br tracer was applied. Bromide concentrations above background levels (0.023 mg/l Br, Doniphan silt loam and 0.016 mg/l Br, Crider silt loam) were detected within 2.60 hours at 0.9 m in the Doniphan soil and within 3.75 hours at that depth in the Crider soil. The rate of Br movement in the profile was greater in both soils than the measured saturated hydraulic conductivity, Bromide concentrations above background levels were present in soil water from the study plots for a minimum of 21 days after irrigation with the Br tracer.     

HISTORICAL CHANGES IN UPPER MISSISSIPPI RIVER WATER AREAS AND ISLANDS1

ABSTRACT: Periodic surveys of the upper Mississippi River since 1866 and a discharge record of nearly equal length provided an opportunity to learn more about the magnitudes and rates of geomorphic processes at work in large stream systems. Furthermore, geomorphic and hydrologic adjustments could be evaluated in relation to watershed land use changes, small‐scale climate fluctuations, and considerable modifications to the channel and floodplain during the period of record. The present study uses GIS mapping to quantitatively compare historical changes in mapped land and water phenomena in the upper Mississippi River Pool 10, located along southwest Wisconsin's border. Modest channel widening and decreases in island area throughout the study reach during the last century are detectable. Flood magnitudes and frequencies also have varied during this time, and stages and low flow discharges have increased since the 1940s. The latter hydrologic change appears to be closely associated with the reach's geomorphic adjustments. Results are representative of a valley reach where a major tributary contributes a large sand bedload, forming an alluvial fan of considerable size in the floodplain.     

Experimental investigation of wellbore integrity and CO2–brine flow along the casing–cement microannulus

Wellbore integrity is one of the key performance criteria in the geological storage of CO₂. It is significant in any proposed storage site but may be critical to the suitability of depleted oil and gas reservoirs that may have 10’s to 1000’s of abandoned wells. Much previous work has focused on Portland cement which is the primary material used to seal wellbore systems. This work has emphasized the potential dissolution of Portland cement. However, an increasing number of field studies (e.g., Carey et al., 2007), experimental studies (e.g., Kutchko et al., 2006) and theoretical considerations indicate that the most significant leakage mechanism is likely to be flow of CO₂ along the casing–cement microannulus, cement–cement fractures, or the cement–caprock interface.In this study, we investigate the casing–cement microannulus through core-flood experiments. The experiments were conducted on a synthetic wellbore system consisting of a 5-cm diameter sample of cement that was cured with an embedded rectangular length of steel casing that had grooves to accommodate fluid flow. The experiments were conducted at 40 $°$ C and 14 MPa pore pressure for 394 h. During the experiment, 6.2 l of a 50:50 mixture of supercritical CO₂ and 30,000 ppm NaCl-rich brine flowed through 10-cm of limestone before flowing through the 6-cm length cement–casing wellbore system. Approximately 59,000 pore volumes of fluid moved through the casing–cement grooves. Scanning electron microscopy revealed that the CO₂–brine mixture impacted both the casing and the cement. The Portland cement was carbonated to depths of 50–250 $\mu$m by a diffusion-dominated process. There was very little evidence for mass loss or erosion of the Portland cement. By contrast, the steel casing reacted to form abundant precipitates of mixed calcium and iron carbonate that lined the channels and in one case almost completely filled a channel. The depth of steel corroded was estimated at 25– $30\mu$m and was similar in value to results obtained with a simplified corrosion model.The experimental results were applied to field observations of carbonated wellbore cement by Carey et al. (2007) and Crow et al. (2009) to show that carbonation of the field samples was not accompanied by significant CO₂–brine flow at the casing–cement interface. The sensitivity of standard-grade steel casing to corrosion suggests that relatively straight-forward wireline logging of external casing corrosion could be used as a useful indicator of flow behind casing. These experiments also reinforce other studies that indicate rates of Portland cement deterioration are slow, even in the high-flux CO₂–brine experiments reported here.     

Phosphorus transport in overland flow in response to position of manure application

Phosphorus (P) loss in overland flow varies with spatial distribution of soil P, management, and hydrological pathways. The effect of flow time, flowpath length, and manure position on P loss in overland flow from two central Pennsylvania soils packed in boxes of varying length (0.5, 1.0, 1.5, 2.75, and 4.0 m long x 15 cm wide x 5 cm deep) were examined by collecting flow samples at 5-min intervals for 30 min (50 mm h(-1) rainfall) without and with 75 kg P ha(-1) applied as swine (Sus scrofa) manure over 0.5 m of the box slope length at distances of 0 to 3.5 m from the downslope collection point. Dissolved reactive P concentration was more closely related to the proportion of clay in sediment of overland flow before (r = 0.98) than after (r = 0.56) manure application. This was attributed to the transport of larger, low-density particles after applying manure. The concentration of dissolved and particulate P fractions decreased with increasing flowpath length, due to dilution rather than sorption of P by surface soil during overland flow. Total P loss (mainly as particulate P) from the Watson channery silt loam (fine-loamy, mixed, active, mesic Typic Fragiudult) was more than from Berks channery silt loam (loamy-skeletal, mixed, active, mesic Typic Dystrudept), even with manure applied. Thus, while P loss in overland flow is affected by where manure is applied relative to flowpath length, initial soil P concentration should not be discounted when looking at areas of potential P loss within a watershed.     

HYDROLOGIC RESPONSE TO PUMPING AND CONTAMINANT ADVECTION IN A FRACTURED ROCK ENVIRONMENT1

ABSTRACT: Ground water flow and supply at the Whiteshell Research Area (WRA) in southeastern Manitoba and the advective movement of contaminants from a hypothetical nuclear fuel waste disposal vault to the adjacent biosphere and a nearby ground water supply well are simulated using finite-element modeling and numerical particle-tracking technique. The hypothetical vault is located at a depth of 500 m, below the water table, in low-permeability plutonic rock of the Canadian Shield. The rock mass is intersected by high-permeability fracture zones (aquifers), which also act as conduits for vault contaminants to migrate to the ground surface. The ground water resource is, therefore, limited in quantity and quality and should be explored with care. A 30 m deep well, which pumps water at a rate of 120 m³/yr from a low-dipping fracture zone, LD1, reduces natural discharge from the system to augment natural recharge. At this pumping rate, a 100 m or 200 m deep well neither reduces discharge nor induces recharge into the system. Thus, at the WRA, a 30 m deep domestic water supply well best meets the water requirements of a one-person household at the rate of 120 m³/yr. The 100 m and 200 m wells best meet the requirements of a family of six and a family of six with light irrigation, respectively, without capturing contaminants’pathways from the vault to the ground surface. By virtue of the proximity of the 200 m well intake to the hypothetical vault, this well performs best as a purge well at pumping rates of 0,000 m³/yr and greater. Finite-element modeling is useful in evaluating the water supply potential of a fractured rock environment in which a nuclear waste disposal vault is proposed to be sited.     

GROUND WATER HYDRAULICS CONSIDERATIONS REGARDING LANDFILLS1

ABSTRACT: Landfill siting and design guidelines or regulations differ from state to state. Most include hydrogeological criteria, referring to hydraulic conductivities, aquifers, ground water flow patterns, contaminant travel times, and distance between landfill and sensitive targets for contaminants, etc. However, almost all of the existing hydrogeological guidelines are incomplete, inconsistent, or both. The aquitard between landfill and regional aquifer frequently offers less resistance to leachate migration than compliance with regulations may suggest. Residence times of leachate, that makes it through the landfill liner, is often overestimated. Monitoring wells in the regional aquifer are unreliable detectors of local leaks in a landfill. If a landfill does leak, costly aquifer restoration is called for. For traditional landfill designs, ground water monitoring considerations suggest the siting over homogeneous sand and gravel aquifers, rather than over complex till environments. An alternative landfill design criterion is suggested, which is based on a negative hydraulic gradient underneath the landfill. This design guarantees ground water protection, simplifies landfill monitoring, and generally enhances the landfill economy.     

ANALYSIS AND EXTRACTION OF LOW FLOW RECESSION CHARACTERISTICS1

ABSTRACT: It is important to extract and assess low flow recession characteristics for water resources planning in the upper reaches of streams. However, it is very difficult to express synthetically the low flow recession characteristics for a stream flow. In this paper, first a new method of constructing the master recession curve based on the exponential expression is proposed and applied with the restriction that there are no regulation or diversion structures in the upper reaches above the measurement station. Daily precipitation and stream flow were used for the analysis. Second, analysis for a recession constant was conducted and the relationship between the recession constant and low flow and/or geology was qualitatively examined. In conclusion, the application of the proposed method indicated that it is objective and useful for constructing the master recession curve. It became apparent that the recession constant of a master recession curve may be defined as the total index of low flow characteristics. In addition, it was found that baseflow value increases in the order of Paleozoic, Mesozoic, Tertiary, and Quaternary.     

NEAR‐STREAM LANDUSE EFFECTS ON STREAMWATER NUTRIENT DISTRIBUTION IN AN URBANIZING WATERSHED1

ABSTRACT: We investigated spatial and temporal relationships among surface and subsurface watershed attributes and stream nutrient concentrations in urbanizing Johnson Creek watershed in northern Oregon. We sampled stream water at eight urban and five nonurban locations from March 1998 through December 1999. We sampled eight wells distributed over the two primary aquifers in the watershed. Using a Geographic Information System (GIS), percentages of landuse attributes within a radius of 30, 91, and 152 m from each sample site were quantified. We analyzed relationships between (1) nutrient concentrations and percentage cover of different landuse attributes, and (2) nutrient concentrations and underlying hydrologic units. We did not find a significant relationship between ground water chemistry and stream water chemistry. We found elevated levels of phosphorus (P) concentrations correlated with urban landuse, while higher nitrogen (N) concentrations were correlated with nonurban (primarily agricultural) landuse. We concluded that elevated levels of N in nonurban areas of Johnson Creek watershed were associated with agricultural practices. We further concluded that urban development factors such as increases in storm drains, dry wells, and impermeable surfaces may be responsible for higher input of P to the stream in urbanizing areas of the Johnson Creek watershed.     

Effect of reservoir flushing on downstream river water quality

The effect of short-term reservoir flushing on downstream water quality in the Geum River, Korea was studied using field experiments and computer simulations. The reservoir release was increased from 30 to 200 m(3)/s within 6 h for the purpose of this experiment. The flushing discharge decreased the concentrations of soluble nitrogen and phosphorus species considerably, but the experimental results revealed a negative impact on organic forms of nutrients and biochemical oxygen demand (BOD). A dynamic river water quality model was applied to simulate the river hydraulics and water quality variations during the event. The model showed very good performance in predicting the travel time of flushing flow and the variations of dissolved forms of nitrogen and phosphorus constituents. However, it revealed a limited capability in simulating organic forms of nutrients and BOD because it does not consider the re-suspension mechanism of these constituents from sediment during the wave front passage.     

PERMEABILITY OF AIR AND IMMISCIBLE ORGANIC LIQUIDS IN POROUS MEDIA1

ABSTRACT: The conductivity of air, mineral oil (relative viscosity 77), and a light nontoxic oil (relative viscosity 4.7) was measured in three porous media: a sand, loamy sand, and a silt loam. The measurements were made over a range of water contents for each porous medium. Small volumes of air were present as well as significant amounts of water during most of the oil conductivity measurements. The results were compared to two methods for calculating conductivities of immiscible fluids in water-wet porous media. A new equation that accounts for swelling and for the gas slippage effect in very small pores was formulated for use with these methods The observed conductivities, spanning seven orders of magnitude, agree reasonably well with calculated values. Only three soil parameters are required to calculate the conductivities: (1) the saturated conductivity of water, (2) the saturated conductivity of the immiscible fluid of interest, and (3) a pore size index value that is obtained from an estimate from the water release curve of the porous material. Remediation of organic liquid spills is briefly discussed to illustrate the practical applications of gas phase conductivities, as well as those for immiscible organic liquid phases. It is concluded that, in light of spatial variation under field conditions, the method presented for calculating values of three-phase conductivities will be useful in the management of immiscible organic liquid spills and leaks.     

THERMAL AND DISSOLVED OXYGEN CHARACTERISTICS OF A SOUTH CAROLINA COOLING RESERVOIR1

ABSTRACT: Temperature and dissolved oxygen concentrations were measured monthly from January 1971 to December 1982 at 1-m depth intervals at 13 stations in Keowee Reservoir in order to characterize spatial and temporal changes associated with operation of the Oconee Nuclear Station. The reservoir water column was i to 4°C warmer in operational than in non-operational years. The thermo-dine was at depths of 5 to 15 m before the operation of Oconee Nuclear Station, but was always below the upper level of the intake (20 m) after the station was in full operation; this suggests that pumping by the Oconee Nuclear Station had depleted all available cool hypolimnetic water to this depth. As a result summer water temperatures at depths greater than 10 m were usually 10°C higher after plant operation began than before. By fall the reservoir was nearly homothemious to a depth of 27 m, where a thermocine developed. Seasonal temperature profiles varied with distance from the plant; a cool water plume was evident in spring and a warm water plume was present in the summer, fall, and winter. A cold water plume also developed in the northern section of the reservoir due to the operation of Jocassee Pumped Storage Station. Increases in the mean water temperature of the reservoir during operational periods were correlated with the generating output of the power plant. The annual heat load to the reservoir increased by one-third after plant operations began. The alteration of the thermal stratification of the receiving water during the summer also caused the dissolved oxygen to mix to greater depths.     

RECIRCULATING WELLS: GROUND WATER REMEDIATION AND PROTECTION OF SURFACE WATER RESOURCES1

ABSTRACT: Several chlorinated solvent plumes threaten the sole‐source aquifer underlying the Massachusetts Military Reservation at the western end of Cape Cod. Sensitive surface water features including ponds, cranberry bogs, and coastal wetlands are hydraulically connected to the aquifer. For one of the plumes (CS‐10 the original remedy of 120 extraction and reinjection wells has the potential for significant disruption of surface water hydrology, through the localized drawdown and mounding of the water table. Recirculating wells with in‐well air stripping offer a cost‐effective alternative to conventional pump‐and‐treat technology that does not adversely affect the configuration of the water table. Pilot testing of a two well system, pumping 300 gpm, showed a capture radius of > 200 feet per well, in‐well trichloroethylene (TCE) removal efficiencies of 92 to 98 percent per recirculation cycle, an average of three recirculation cycles within the capture zone, and no measurable effect on water table elevations at any point within the recirculation/treatment zone. During 120 days of operation, the mean concentration of TCE in the treatment zone was reduced by 83 percent, from 1,111 μg/l to 184 μg/l. Full‐scale design projections indicate that 60 wells at an average spacing of 160 feet, having an aggregate recirculation 11 MGD, can contain the CS‐b plume without ground water extraction or adverse hydraulic effects on surface water resources. The estimated capital costs for such a system are about $7 million, and annual operations‐and‐maintenance costs should be about $1.4 million, 40 percent of those associated with a pump and treat system over a 20‐year period.     

A Stochastic Modeling Approach to Address Hydrogeologic Uncertainties in Modeling Wellhead Protection Boundaries in Karst Aquifers1

Abstract: Development of any numerical ground‐water model is dependent on hydrogeologic data describing the subsurface. These data are obtained from geologic core analyses, stratigraphic analyses, aquifer performance tests, and geophysical studies. But typically in remote areas, these types of data are very sparse and site‐specific in terms of the aerial extent of the resource to be modeled. Uncertainties exist as to how well the available data from a few locations defines a heterogeneous surficial aquifer such as the Biscayne Aquifer in Miami‐Dade County, Florida. This is particularly the case when an exceptionally conductive horizontal flow zone is detected at one site due to specialized testing that was not historically conducted at the other at sites that provided data for the model. Not adequately accounting for the potential effect of the high flow zone in the aquifer within a ground‐water numerical model, even though the zone may be of very limited thickness, might underpredict the well field protection capture boundaries. Applied Stochastic ground‐water modeling in determining well field protection zones is steadily becoming important in addressing the uncertainty of the hydrogeologic subsurface parameters, specifically in karstic heterogeneous aquifers. This is particularly important in addressing the uncertainty of a 60‐day travel time capture zone in the Northwest Well Field, Miami‐Dade County, where a predominantly high flow zone controls much of the flow in the production wells. A stochastic ground‐water modeling application along with combination of pilot points and regularization technique is presented to further consolidate the uncertainty of the subsurface.