Geology ReportsSearch

Geology topics

W.N. Herkelrath

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

12 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

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