USGS Science⌕ Search

Geology topics

W.N. Herkelrath

Publications and source records attributed to W.N. Herkelrath.

16 recordsLinked to original sources

Simulation of fluid distributions observed at a crude oil spill site incorporating hysteresis, oil entrapment, and spatial variability of hydraulic properties

Subsurface oil, water, and air saturation distributions were determined using 146 samples collected from seven boreholes along a 120-m transect at a crude oil spill site near Bemidji, Minnesota. The field data, collected 10 years after the spill, show a clearly defined oil body that has an oil saturation distribution that appears to be influenced by sediment heterogeneities and water table fluctuations. The center of the oil body has depressed the water-saturated zone boundary and the oil appears to have migrated laterally within the capillary fringe. A multiphase cross-sectional flow model was developed and used to simulate the movement of oil and water at the spill site. Comparisons between observed and simulated oil saturation distributions serve as an indicator of the appropriateness of using such models to predict the actual spread of organic immiscible liquids at spill sites. Sediment hydraulic properties used in the model were estimated from particle size data. The general large-scale features of the observed oil body were reproduced only when hysteresis with oil entrapment and representations of observed spatial variability of hydraulic properties were incorporated into the model. The small-scale details of the observed subsurface oil distribution were not reproduced in the simulations. The discrepancy between observed and simulated oil distributions reflects the considerable uncertainty in model parameter estimates and boundary conditions, three-phase capillary pressure-saturation-relative permeability functions, representations of spatial variability of hydraulic properties, and hydrodynamics of the groundwater flow system at the study site.

Water Resources Research↗

Laboratory investigations of steam flow in a porous medium

Experiments were carried out in the laboratory to test a theory of transient flow of pure steam in a uniform porous medium. This theory is used in modeling pressure transient behavior in vapor dominated geothermal systems. Transient, superheated steam flow experiments were run by bringing a cylinder of porous material to a uniform initial pressure and then making a step increase in pressure at one end of the sample while monitoring the pressure transient breakthrough at the other end. It was found in experiments run at 100°, 125°, and 146°C that the time required for steam pressure transients to propagate through an unconsolidated material containing sand, silt, and clay was 10–25 times longer than predicted by conventional superheated steam flow theory. It is hypothesized that the delay in the steam pressure transient was caused by adsorption of steam in the porous sample. In order to account for steam adsorption, a sink term was included in the conservation of mass equation. In addition, energy transfer in the system has to be considered because latent heat is released when steam adsorption occurs, increasing the sample temperature by as much as 10°C. Finally, it was recognized that the steam pressure was a function of both the temperature and the amount of adsorption in the sample. This function was assumed to be an equilibrium adsorption isotherm, which was determined by experiment. By solving the modified mass and energy equations numerically, subject to the empirical adsorption isotherm relationship, excellent theoretical simulation of the experiments was achieved.

Water Resources Research↗

A sample-freezing drive shoe for a wire line piston core sampler

Loss of fluids and samples during retrieval of cores of saturated, noncohesive sediments results in incorrect measures of fluid distributions and an inaccurate measure of the stratigraphic position of the sample. To reduce these errors, we developed a hollow drive shoe that freezes in place the lowest 3 inches (75 mm) of a 1.88‐inch‐diameter (48 mm), 5‐foot‐long (1.5 m) sediment sample taken using a commercial wire line piston core smapler. The end of the core is frozen by piping liquid carbon dioxide at ambient temperature through a steel tube from a bottle at the land surface to the drive shoe where it evaporates and expands, cooling the interior surface of the shoe to about ‐ 109°F (‐ 78°C). Freezing a core end takes about 10 minutes. The device was used to collect samples for a study of oil‐water‐air distributions, and for studies of water chemistry and microbial activity in unconsolidated sediments at the site of an oil spill near Bemidji, Minnesota. Before freezing was employed, samples of sandy sediments from near the water table sometimes flowed out of the core barrel as the sampler was withdrawn. Freezing the bottom of the core allowed for the retention of all material that entered the core barrel and lessened the redistribution of fluids within the core. The device is useful in the unsaturated and shallow saturated zones, but does not freeze cores well at depths greater than about 20 feet (6 m) below water, possibly because the feed tube plugs with dry ice with increased exhaust back‐pressure, or because sediment enters the annulus between the core barrel and the core barrel liner and blocks the exhaust.

Ground Water Monitoring and Remediation↗

Use of soil moisture probes to estimate ground water recharge at an oil spill site

Soil moisture data collected using an automated data logging system were used to estimate ground water recharge at a crude oil spill research site near Bemidji, Minnesota. Three different soil moisture probes were tested in the laboratory as well as the field conditions of limited power supply and extreme weather typical of northern Minnesota: a self‐contained reflectometer probe, and two time domain reflectometry (TDR) probes, 30 and 50 cm long. Recharge was estimated using an unsaturated zone water balance method. Recharge estimates for 1999 using the laboratory calibrations were 13 to 30 percent greater than estimates based on the factory calibrations. Recharge indicated by the self‐contained probes was 170 percent to 210 percent greater than the estimates for the TDR probes regardless of calibration method. Results indicate that the anomalously large recharge estimates for the self‐contained probes are not the result of inaccurate measurements of volumetric moisture content, but result from the presence of crude oil, or borehole leakage. Of the probes tested, the 50 cm long TDR probe yielded recharge estimates that compared most favorably to estimates based on a method utilizing water table fluctuations. Recharge rates for this probe represented 24 to 27 percent of 1999 precipitation. Recharge based on the 30 cm long horizontal TDR probes was 29 to 37 percent of 1999 precipitation. By comparison, recharge based on the water table fluctuation method represented about 29 percent of precipitation.

Journal of the American Water Resources Associatio↗

Progression of methanogenic degradation of crude oil in the subsurface

Our results show that subsurface crude-oil degradation rates at a long-term research site were strongly influenced by small-scale variations in hydrologic conditions. The site is a shallow glacial outwash aquifer located near Bemidji in northern Minnesota that became contaminated when oil spilled from a broken pipeline in August 1979. In the study area, separate-phase oil forms a subsurface oil body extending from land surface to about 1 m (3.3 ft) below the 6–8-m (20–26 ft)-deep water table. Oil saturation in the sediments ranges from 10–20% in the vadose zone to 30–70% near the water table. At depths below 2 m (6.6 ft), degradation of the separate-phase crude oil occurs under methanogenic conditions. The sequence of methanogenic alkane degradation depletes the longer chain n-alkanes before the shorter chain n-alkanes, which is opposite to the better known aerobic sequence. The rates of degradation vary significantly with location in the subsurface. Oil-coated soils within 1.5 m (5 ft) of land surface have experienced little degradation where soil water saturation is less than 20%. Oil located 2–8 m (6.6–26 ft) below land surface in areas of higher recharge has been substantially degraded. The best explanation for the association between recharge and enhanced degradation seems to be increased downward transport of microbial growth nutrients to the oil body. This is supported by observations of greater microbial numbers at higher elevations in the oil body and significant decreases with depth in nutrient concentrations, especially phosphorus. Our results suggest that environmental effects may cause widely diverging degradation rates in the same spill, calling into question dating methods based on degradation state.

Environmental Geosciences↗

Impacts of petroleum production on ground and surface waters: Results from the Osage-Skiatook Petroleum Environmental Research A site, Osage County Oklahoma

Kharaka has been a research hydrogeochemist with the U.S. Geological Survey since1975. His current research covers the fields of water-rock-gas interactions, CO 2 sequestration, contamination from agricultural drainage and petroleum produced water, and fluid-fault interactions. Kharaka has authored more than 100 scientific papers and book chapters and has delivered about 200 presentations. He received his Ph.D. from the University of California, Berkeley, in 1971. Thordsen is a geologist, employed with the U.S. Geological Survey since 1981 and presently working in the project team Hydrogeochemistry of Fine Grained Sediments, in the Water Resources Division, National Research Program. His recent research covers water-rock interactions and geochemical modeling in geothermal systems, fault-fluid systems, and petroleum produced water brine impacts. He received his M.S. in geology from the Ohio State University in 1988. Kakouros has been a hydrogeologist in the U.S. Geological Survey Water Resources Division in Menlo Park, California, since 1999. His current research covers the fields of Se contamination from agricultural drainage, CO 2 sequestration, and contamination from petroleum and petroleum-produced water. Kakouros received his M.Sc. from the San Jose State University in 2001. Herkelrath is a physicist in the U.S. Geological Survey's Water Resources Division National Research Program in Menlo Park, California. Herkelrath received his Ph.D. in physics from the University of Wisconsin at Madison. His recent research interests include multiphase flow and solute transport modeling in shallow ground-water systems where oil or brine is present.

Oklahoma↗

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m , p -xylene, and o -xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o -xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m , p -xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o -xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.

Minnesota↗

Ground-water contamination by crude oil: Section B in U.S. Geological Survey Toxic Substances Hydrology Program: Proceedings of the technical meeting, Charleston, South Carolina, March 8-12, 1999: Volume 3 (Part C)

Ground water contamination by crude oil, and other petroleum-based liquids, is a widespread problem. An average of 83 crude-oil spills occurred per year during 1994-96 in the United States, each spilling an average of about 50,000 barrels of crude oil (U.S. Office of Pipeline Safety, electronic commun., 1997). An understanding of the fate of organic contaminants (such as oil and gasoline) in the subsurface is needed to design innovative and cost-effective remedial solutions at contaminated sites.

Water-Resources Investigations Report↗

Determination of subsurface fluid contents at a crude-oil spill site

Measurement of the fluid-content distribution at sites contaminated by immiscible fluids, including crude oil, is needed to better understand the movement of these fluids in the subsurface and to provide data to calibrate and verify numerical models and geophysical methods. A laboratory method was used to quantify the fluid contents of 146 core sections retrieved from boreholes aligned along a 120-m longitudinal transect at a crude-oil spill site near Bemidji, Minnesota, U.S.A. The 47-mm-diameter, minimally disturbed cores spanned a 4-m vertical interval contaminated by oil. Cores were frozen on site in a dry ice-alcohol bath to prevent redistribution and loss of fluids while sectioning the cores. We gravimetrically determined oil and water contents using a two-step method: (1) samples were slurried and the oil was removed by absorption onto strips of hydrophobic porous polyethylene (PPE); and (2) the samples were oven-dried to remove the water. The resulting data show sharp vertical gradients in the water and oil contents and a clearly defined oil body. The subsurface distribution is complex and appears to be influenced by sediment heterogeneities and water-table fluctuations. The center of the oil body has depressed the water-saturated zone boundary, and the oil is migrating laterally within the capillary fringe. The oil contents are as high as 0.3cm 3 cm −3 , which indicates that oil is probably still mobile 10 years after the spill occurred. The thickness of oil measured in wells suggests that accumulated thickness in wells is a poor indicator of the actual distribution of oil in the subsurface. Several possible sources of error are identified with the field and laboratory methods. An error analysis indicates that adsorption of water and sediment into the PPE adds as much as 4% to the measured oil masses and that uncertainties in the calculated sample volume and the assumed oil density introduce an additional ±3% error when the masses are converted to fluid contents.

Journal of Contaminant Hydrology↗

Air encapsulation during infiltration

A series of field and laboratory experiments were performed to measure the effects of air encapsulation within the soil's transmission zone upon several infiltration properties. In the field, infiltration rates were measured using a double-cap infiltrometer (DCI), and soil-water contents were measured using time-domain reflectometry (TDR). Before half of the infiltration experiments, CO 2 was injected through the DCI into the soil to reduce the amount of air encapsulation in the soil's transmission zone. For a gravelly loam as steady infiltration rates were approached, the average volumetric water content was 0.38 cm 3 cm −3 for control experiments and 0.43 cm 3 cm −3 for CO 2 experiments. The average steady infiltration rate was 0.42 cm min −1 for the control experiments compared to 4.40 cm min −1 for the CO 2 experiments. For a sandy loam as steady infiltration rates were approached, the average volumetric water content was 0.43 cm 3 cm −3 for control experiments compared with 0.45 cm 3 cm −3 for CO 2 experiments. The average final infiltration rate was 0.09 cm min −1 for the control experiments compared with 0.42 cm min −1 for the CO 2 experiments. In the laboratory, infiltration experiments were performed using repacked soil columns (15-cm i.d. by 140 cm long), again using TDR and CO 2 flooding. For a medium sand as steady infiltration rates were approached, the average volumetric water content was 0.29 cm 3 cm −3 for the control experiments and 0.36 cm 3 cm −3 for the CO 2 experiments. The average steady infiltration rate was 0.25 cm min −1 for the control experiments and 1.23 cm min −1 for the CO 2 experiments. For a loam as steady infiltration rates were approached, the average volumetric water content was 0.45 cm 3 cm −3 for the control experiments and 0.50 cm 3 cm −3 for the CO 2 experiments. The average steady infiltration rate was 0.02 cm min −1 for the control experiments and 0.10 cm min −1 for the CO 2 experiments. These results suggest that a significant portion of the total encapsulated air resided in interconnected pores within the soil's transmission zone. For the time scale considered, this residual air caused the effective hydraulic conductivity of the transmission zone to remain at a level no greater than 20% of the saturated hydraulic conductivity of the soil.

Soil Science Society of America Journal↗

Submersible pressure outflow cell for measurement of soil water retention and diffusivity from 5 to 95 degrees C

A technique was developed to measure water content in soil as a function of capillary pressure from 5 to 95°C. To overcome problems encountered at high temperature, a modified Tempe pressure cell containing a soil sample is suspended in a constant-temperature water bath. The cell's porous plate is in direct contact with circulating bath water, thus eliminating the problem of entrapped air bubbles. A balance located above the water bath measures water content changes in the soil by weighing the entire pressure cell under water. The technique is designed to measure soil water retention characteristics and to make transient outflow estimates of the soil water diffusivity at temperatures from 5 to 95°C. We also used the technique to determine the isobaric temperature dependence of water retention in soil. Results indicate that at constant capillary pressure, the relationship between moisture content and temperature is hysteretic.

Soil Science Society of America Journal↗

Laboratory investigations of the physics of steam flow in a porous medium

Experiments were carried out in the laboratory to test a theory of transient flow of pure steam in a uniform porous medium. This theory is used extensively in modeling pressure-transient behavior in vapor-dominated geothermal systems. Transient, superheated steam-flow experiments were run by bringing a cylinder of porous material to a uniform initial pressure, and then making a step increase in pressure at one end of the sample, while monitoring the pressure-transient breakthrough at the other end. It was found in experiments run at 100?, 125?, and 146?C that the time required for steam-pressure transients to propagate through an unconsolidated material containing sand, silt, and clay was 10 to 25 times longer than predicted by theory. It is hypothesized that the delay in the steam-pressure transient was caused by adsorption of steam in the porous sample. In order to account for steam adsorption, a sink term was included in the conservation of mass equation. In addition, energy transfer in the system has to be considered because latent heat is released when steam adsorption occurs, increasing the sample temperature by as much as 10?C. Finally, it was recognized that the steam pressure was a function of both the temperature and the amount of adsorption in the sample. For simplicity, this function was assumed to be in equilibrium adsorption isotherm, which was determined by experiment. By solving the modified mass and energy equations numerically, subject to the empirical adsorption isotherm relationship, excellent theoretical simulation of the experiments was achieved. The experiments support the hypothesis that adsorption of steam can strongly influence steam pressure-transient behavior in porous media; the results suggest that the modified steam-flow theory, which includes steam adsorption terms, should be used in modeling steam flow in vapor-dominated geothermal systems.

Open-File Report↗