Search

In the United States, hundreds of thousands of undocumented orphan wells have been abandoned, leaving the burden of managing environmental hazards to governmental agencies or the public. These wells, a result of over a century of fossil fuel extraction without adequate regulation, lack basic information like location and depth, emit greenhouse gases, and leak toxic substances into groundwater. For most of these wells, basic information such as well location and depth is unknown or unverified. Addressing this issue necessitates innovative and interdisciplinary approaches for locating, characterizing, and mitigating their environmental impacts. Our survey of the United States revealed the need for tools to identify well locations and assess conditions, prompting the development of technologies including machine learning to automatically extract information from old records (95%+ accuracy), remote sensing technologies like aero-magnetometers to find buried wells, and cost-effective methods for estimating methane emissions. Notably, fixed-wing drones equipped with magnetometers have emerged as cost-effective and efficient for discovering unknown wells, offering advantages over helicopters and quadcopters. Efforts also involved leveraging local knowledge through outreach to state and tribal governments as well as citizen science initiatives. These initiatives aim to significantly contribute to environmental sustainability by reducing greenhouse gases and improving air and water quality.

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. We use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

Continuous microseismic monitoring using three-component (3C) accelerometers deployed in multiple boreholes allows for tracking the detailed evaluation of mesoscale (∼10 m scale) fracture growth during the fracture stimulation experiments at the first Enhanced Geothermal Systems (EGS) Collab testbed. Building on a well-constrained microseismic event catalog, we invert for moment tensor of the events to better understand the fracture geometry and stress orientations. However, it is challenging because of the unknown orientation of 3C accelerometers and low signal-to-noise-ratio nature of high-frequency (several kHz) monitoring. To address these challenges, we first perform the hodogram analysis on the continuous active-source seismic monitoring (CASSM) data to determine the orientations of the 18 3C accelerometers. We then apply the principal component analysis (PCA) to the observed microseismic waveforms to improve the signal-to-noise ratios. We perform a grid search for the full moment tensor by fitting the PCA-denoised waveforms at a frequency range of 5 to 8 kHz. The moment tensor results show both the creation of hydraulic fractures and the reactivation of natural fractures during the hydraulic stimulations. Our stress inversion based on the inverted moment tensors reveals the alteration of stress regime caused by hydraulic fracture stimulations.

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

The water-gas relative permeability is an important parameter to characterize multiphase flow in sediments. To study the water-gas relative permeability of hydrate-bearing porous media, multiphase flow simulations were carried out at the pore scale using the lattice Boltzmann method. The effects of hydrate saturation and hydrate-growth habits on the water-gas relative permeability, which is scaled by the relative permeability considering the hydrate only, were evaluated in a two-dimensional porous medium. Results show that the increase of hydrate saturation causes the decrease of water-gas effective permeability as expected. However, the effect of hydrate saturation on the water-gas relative permeability is different from that of hydrate saturation on the water-gas effective permeability. The water-gas relative permeability increases with the increase of hydrate saturation in the pore-filling case. The water-gas relative permeability decreases with the increase of hydrate saturation in the grain-coating case. The wettability of solid phase has a different effect on the relative permeability of wetting phase and nonwetting phase. The Jamin effect (phase blocking) was observed and may exist in the production of gas from natural gas hydrate reservoirs. This seriously affects the multiphase flow characteristics. The changes of microscale fluid distribution effect the changes of water-gas relative permeability. The relationship between the water-gas relative permeability and the characterization parameters of microscale fluid distribution was analyzed.

Heat recovery from an enhanced geothermal system (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS commercialization but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150-250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: This increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: A thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. These experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, we show the classic fingering and transitions in flow regimes at (Formula presented.) (capillary number) from −7.2 to −3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO2 (scCO2) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

Estimation of CO2 storage capacity in fractured porous reservoirs requires a better understanding of the CO2-water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO2-water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. The core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO2-water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO2 saturations in the matrix increaze with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO2 storage efficiency in fractured porous media.

The final version of the above article was posted prematurely on 16 July 2021, owing to a technical error. The final, corrected version of record will be made fully available at a later date.

In cryogenic fracturing, a rock surface exposed to cryogenic fluids undergoes a large thermal gradient, and the resultant local tensile stress overcomes rock strength and initiates fractures. This study investigates the development of cracks generated from the cryogenic treatment of a borehole under no external confining stress on specimens. The experiments were performed on transparent PMMA specimens to observe fracture proliferation around boreholes. Liquid nitrogen was flowed through the boreholes to cool the borehole surface. The results show that initial fracture growth is characterized by abrupt starts and stops, and as the fracture propagates outward, the growth appears more continuous. In an early stage, horizontal/radial fractures and vertical fractures are the defining patterns. Horizontal fractures tend to be separated by a specific exclusion distance (i.e., spacing between cracks). While distinct horizontal/vertical fractures and exclusion distance manifest themselves at an early stage, fractures resulting from fracture interactions and curvatures can develop into complex shapes at later stages. Cryogenic thermal loading induces distinctively curved fractures. The tendency of curvature may prevent greater penetration. An increase in the borehole pressure during liquid nitrogen flow, however, can lessen fracture tortu-osity and facilitate radial propagation. A high flow pressure and rate are also advantageous in that they accelerate cooling and fracture propagation.

There is an interest to ensure sub-saturated water content in lines containing carbon dioxide in applications such as enhanced oil recovery and carbon sequestration, to reduce risks of hydrate blockage and corrosion. The water content of carbon dioxide at various temperatures and pressures has been measured in the past, but there is no consistent set of measurements that could be used for carbon dioxide storage and transportation design work. The solubility of water in a carbon dioxide rich gas phase at hydrate forming conditions was measured in this work. Pressures ranged from 12.06 to 29.30 bar along two isotherms, 1 °C and −7 °C, all within the gaseous carbon dioxide and hydrate stability zone. For the first time in these types of measurements, the solid phase was also characterized and confirmed to be carbon dioxide hydrate via X-ray computed tomography, simultaneous with water content measurements of the gas phase. Once carbon dioxide hydrate conversion had reached a maximum value (65% estimated by X-ray computed tomography), the equilibrium water content was measured. Prior to reaching this maximum carbon dioxide hydrate conversion, the water content in carbon dioxide was observed to decrease as liquid water converted to carbon dioxide hydrate. This slow conversion to hydrate, metastability of the hydrate phase, or unexpected phases may be responsible for the large discrepancy between prior data sets for similar carbon dioxide water content measurements.

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO2 (scCO2) is used in enhanced oil recovery. We performed scCO2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from −8.1 to −4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO2 saturation increase at both capillary fingering and crossover flow regimes (−8.1 ≤ logCa ≤ − 6.1), the one heterogeneous wetting to scCO2 results in ~0.09 saturation increase only at the crossover flow regime (−7.1 ≤ logCa ≤ − 6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > − 6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO2 saturation and distribution suggest the importance of wettability on CO2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiagvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

Flow in natural media with nanoporous regions is very complicated, with many governing processes. Welldeveloped numerical codes to integrate and model flow through these media are available. A great deal of very creative work has been done to understand individual processes governing flow in natural media with nanoporous regions, and both simple and elaborate tools have been used to gain process understanding. The complications in understanding and predicting flow in these rocks, particularly with multiscale heterogeneities, anisotropies, and the presence of multiple phases and large gradients are enormous. Here, we examine factors governing flow through natural porous media containing nanoporous regions. We present a conceptual model of the media, touch on the flow physics, and describe the techniques used to examine pore space in these rocks. In addition, we briefly describe some modeling of flow through these media. A number of processes which need better description are identified.

Geological carbon storage in deep aquifers involves displacement of resident brine by supercritical CO2 (scCO2), which is an unstable drainage process caused by the invasion of less viscous scCO2. The unstable drainage is greatly complicated by aquifer heterogeneity and anisotropy and regarded as one of the key factors accounting for the uncertainty in storage capacity estimates. The impacts of pore-scale characteristics on the unstable drainage remain poorly understood. In this study, scCO2 drainage experiments were conducted at 40 °C and 9 MPa using a homogeneous elliptical micromodel with low or high anisotropy, a homogeneous/isotropic hexagonal micromodel, and a heterogeneous sandstone-analog micromodel. Each initially water-saturated micromodel was invaded by scCO2 at different rates with logCa (the capillary number)ranging from −7.6 to −4.4, and scCO2/water images were obtained. The measured CO2 saturations in these centimeter-scale micromodels vary considerably from 0.08 to 0.93 depending on the pore-scale characteristics and capillary number. It was also observed that scCO2 drainage follows the classic flow-regime transition from capillary fingering through crossover to viscous fingering for either of the low-anisotropy elliptical and heterogeneous micromodels, but with disparate crossover zones. The crossover zones of scCO2 saturation were then unified with the minimum scCO2 saturation occurring at logCa*=-4.0 using the complete capillary number (Ca*)that considers pore characteristics. For the hexagonal and the high-anisotropy elliptical micromodels, a monotonic increase in scCO2 saturation with increasing Ca* (without crossover)was observed. It appears that the complete capillary number is more appropriate than the classic capillary number when characterizing flow regimes and CO2 saturation in different pore networks.

Describing the pore habit of methane hydrate in sediment matrices is essential for understanding natural distribution of methane hydrate, methane trace (transport and solidification) in the hydrate stability zone, physical properties of hydrate-bearing sediments, and the associated influence on potential gas production. Pore habit visualization in natural media at pore scale even with laboratory synthesized cores has been challenging due to the similar densities of methane hydrate from pore liquid. In this work, we used phase-contrast assisted micro-CT with potassium iodine-doped brine to visualize four phases: sand particles, pore fluid, methane hydrate and methane gas. This study visualizes the pore habit of methane hydrate at various stages including during hydrate formation in excess-gas systems, its evolution after brine injection to replace pore fluid, and hydrate formation in excess-water systems. Hydrate tends to adopt round and smooth surfaces when in contact with water while exhibits relatively angular interfaces when in contact with methane gas. Hydrate formation in excess-gas systems results in a partial cementing and partial mineral-coating pore habit, while hydrate in excess-water systems develops mainly as pore-filling, and locally cementing or mineral-coating where big gas pockets exist at the initial state. Pore liquid replacement from methane gas to brine triggers a shift of hydrate pore habit towards pore-filling. Methane hydrate evolution over time produces bigger hydrate particles but with less contact area with sand particles. The effects of hydrate pore habit become less important as hydrate particle size exceeds the pore size. Additionally, hydrate formation could trap residual methane gas and brine as inclusions.

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrology models. Field-scale investigations were conducted at an intensively-instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Samplers arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table, then was transported horizontally. Horizontal flow occurred in more locations and at higher velocities of fluxes in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

Dissolution trapping is one of the most important mechanisms for geological carbon storage (GCS). Recent laboratory and field experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and coupled scCO2 dissolution and water flow, i.e., scCO2 dissolution at local pores/pore throats creating new water-flow paths, which in turn enhance dissolution by increased advection and interfacial area. However, the impacts of pore-scale characteristics on these coupled processes have not been investigated. In this study, imbibition and dissolution experiments were conducted under 40 °C and 9 MPa using a homogeneous/isotropic hexagonal micromodel, two homogeneous elliptical micromodels with low or high anisotropy, and a heterogeneous sandstone-analog micromodel. The four micromodels, initially saturated with deionized (DI)-water, were drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at different rates with logCa (the capillary number) ranging from −6.56 to −4.34. Results show that bypass of scCO2 by displacing water is the dominant mechanism contributing to the residual CO2 trapping, triggered by heterogeneity in pore characteristics or pore-scale scCO2-water distribution. Bypass can be enhanced by pore heterogeneity or reduced by increasing transverse permeability, resulting in relatively low (

A laboratory study of cryogenic fracturing was performed to test its ability to improve oil/gas recovery from low-permeability reservoirs. Our objective is to develop well-stimulation technologies using cryogenic fluids [e.g., liquid nitrogen (LN)] to increase permeability in a large reservoir volume surrounding wells. The new technology has the potential to reduce formation damage caused by current stimulation methods and minimize or eliminate water usage. The concept of cryogenic fracturing is that a sharp thermal gradient (thermal shock) created at the surfaces of formation rocks by applying cryogenic fluid can cause strong local tensile stress and start fractures. We developed a laboratory system for cryogenic fracturing under true-triaxial loading, with LN-delivery/control and -measurement systems. The loading system simulates confining stresses by independently loading each axis up to approximately 5,000 psi on 888-in. cubes. Temperature in boreholes and at block surfaces and fluid pressure in boreholes were continuously monitored. Acoustic and pressure-decay measurements were obtained before and at various stages of stimulations. Cubic blocks (8 88-in.) of Niobrara shale, concrete, and sandstones were tested, and stress levels and anisotropies varied. Three schemes were considered: Gas fracturing without cryo-stimulation, gas fracturing after low-pressure cryogen flow-through, and gas fracturing after high-pressure cryogen flow-through. Results from pressure-decay tests show that LN stimulation clearly increases permeability, and repeated stimulations further increase the permeability. Acoustic velocities and amplitudes decreased significantly after cryo-stimulation, indicating fracture creation. In the gas fracturing without the stimulation, breakdown (complete fracturing) occurs suddenly without any initial leaking, and major fracture planes form along the plane containing principal-stress and intermediate-stress directions, as expected theoretically. However, in the gas fracturing after cryogenic stimulations, breakdown occurred gradually and with massive leaking because of thermal fractures created during stimulation. In addition, the major fracture direction does not necessarily follow the plane containing the principal-stress direction, especially at low confining-stress levels. In tests, we observed that cryogenic stimulation seems to disrupt the internal stress field. The increase in borehole temperature after stimulation affects the permeability of the specimen. When a stimulated specimen is still cold, it maintains high permeability because fractures remain open and local thermal tension is maintained near the borehole. When the rock warms back, fractures close and permeability decreases. In these tests, we have not used proppants. Overall, fractures are clearly generated by low- and high-pressure thermal shocks. The added pressure of the high-pressure thermal shocks helps to further propagate cryogenic fractures generated by thermal shock. Breakdown pressure is significantly lowered by LN stimulation, with observed breakdown-pressure reductions up to approximately 40%.

We form methane hydrate by injecting methane gas into a brine-saturated, coarse-grained sample under hydrate-stable thermodynamic conditions. Hydrate forms to a saturation of 11%, which is much lower than that predicted assuming three-phase (gas-hydrate-brine) thermodynamic equilibrium (67%). During hydrate formation, there are temporary flow blockages. We interpret that a hydrate skin forms a physical barrier at the gas-brine interface. The skin fails periodically when the pressure differential exceeds the skin strength. Once the skin is present, further hydrate formation is limited by the rate that methane can diffuse through the solid skin. This process produces distinct thermodynamic states on either side of the skin that allows gas to flow through the sample. This study illuminates how gas can be transported through the hydrate stability zone and thus provides a mechanism for the formation of concentrated hydrate deposits in sand reservoirs. It also illustrates that models that assume local equilibrium at the core-scale and larger may not capture the fundamental behaviors of these gas flow and hydrate formation processes.

For better understanding of frequency dependence (dispersion) of seismic wave velocities caused by stress-induced fluid flow, broadband laboratory measurements were performed on a suite of synthetic glass media containing both equant pores and thermal cracks. Complementary forced oscillation, resonant bar, and ultrasonic techniques provided access to millihertz-hertz frequencies, ~1 kHz frequency, and ~1 MHz frequency, respectively. The wave speeds or effective elastic moduli and associated dissipation were measured on samples under dry, argon- or nitrogen-saturated, and water-saturated conditions in sequence. The elastic moduli, in situ permeability, and crack porosity inferred from in situ X-ray computed tomography all attest to strong pressure-induced crack closure for differential (confining-minus-pore) pressures

Depressurization of hydrate-bearing sediments (HBS) can cause the movement of fine particles, and in turn, such fines migration affects fluid flow and mechanical behavior of sediments, ultimately affecting long-term hydrocarbon production and wellbore stability. This study investigated how and to what extent depressurization of HBS causes fines migration using X-ray computed tomography (CT) imaging. Methane hydrate was synthesized in sediments with 10% fines content (FC), composed of sands with silt and/or clay, and the hydrate-bearing samples were stepwisely depressurized while acquiring CT images. The CT images were analyzed to quantify the spatial changes in FC in the host sediment and thus to capture the fines migration during depressurization. It was found that the FC changes began occurring from the hydrate dissociation regions. This confirms that the multiphase flow caused by depressurization accompanies fines migration. Depressurization of HBS with a hydrate saturation of ~20–40% caused FC reduction from ~10% to ~6–9%, and the extent of fines migration differed with the particle sizes of the host sands and the types of fines. It was found that fines migration was more pronounced with coarse sands and with silty fines. Such observed level of FC reduction is estimated to increase sediment permeability by several factors based on the Kozeny-type permeability model. Our results support the notion that the extent of fines migration and its effect on fluid flow behavior need to be assessed in consideration of physical properties of host sediment and fine particles to identify optimum depressurization strategies.

Niobrara shale cubes of 20 cm from Colorado were employed to investigate gas and supercritical CO2 injection-induced fracturing in naturally fractured caprocks of deep aquifers/depleted reservoirs and fractured shale reservoirs. Under tri-axial stresses, gas or supercritical CO2 was injected into the center of the cubes to induce fracturing. Real-time pressure and temperature, acoustic wave, pressure decay, fracture coloring, and gas fracturing were used to characterize the fracturing process and fracture morphology. Without pore pressure, CO2 injection-induced fracturing occurred and completed instantly, accompanied by an evident temperature drop. Strongly bonded fractures barely affected transverse fracture propagation, whereas weakly bonded or open fractures arrested the injected fluid first and then allowed it to generate new fractures perpendicular to the minimum horizontal stress. Breakdown pressures for cubes with preexisting fractures using gas and supercritical CO2 are much lower than both poroelastic predictions and slick-water fracturing pressure, and some are even lower than the minimum horizontal stress. This is attributed to unconformable preexisting fractures and the low viscosity of CO2. Moreover, decreasing tri-axial stress levels and increasing stress differences tend to lower the breakdown pressure. This study is instructive for understanding and tackling geomechanical issues related to CO2 geological storage and fracturing of shale reservoirs.

Leakage risk assessment is an inevitable procedure in permanent sequestration and storage of CO2 in deep saline aquifers and depleted oil and gas reservoirs, where the integrity of caprock is most critical. Low porosity and low permeability concrete cubes were employed as caprock analogs to investigate the supercritical CO2 injection-induced fracturing processes under true tri-axial stress conditions. A systematic experimental procedure, consisting of active acoustic emission measurement, pressure decay, injection pressure and temperature monitoring, fracture coloring, and gas fracturing, is formulated to qualitatively and quantitatively characterize the injection-induced fracturing processes and fracture morphology. Occurrence of injection-induced fracturing can be directly identified from peaks of borehole pressure profiles as well as sharp drops on temperature profiles due to CO2 expansion, but generally there was no fracture propagation plateau appearing for the 20 cm rock cubes with zero pore pressure. Acoustic wave signatures, including both waveform change and arrival time delay, can effectively capture the extension of induced fractures inside the opaque rock. Initiation and propagation of supercritical CO2 injection-induced fractures are highly dominated by the tri-axial stresses, following the general trend of continuum mechanics at high stress levels with large stress difference. As the stress difference decreases, induced fractures branch off in relatively arbitrary directions. For these supercritical CO2 injection-induced fracturing experiments, poroelastic mechanics model makes a decent fit with the measured breakdown pressure of the caprock analogs. Findings in this study are valuable for risk analysis and operation optimization of geological CO2 sequestration and storage as well as for CO2 fracturing design and implementation in shale and tight reservoirs.

The concept of cryogenic fracturing is that a sharp thermal gradient developed by applying a cryogenic fluid on a rock surface causes a strong local tensile stress that initiates fractures. Prior field tests suggest that field application with special equipment rated for cryogenic temperatures may deliver potential benefits. The tests did not, however, identify the fracture mechanisms at work in downhole conditions. In this study, we present our laboratory designs and procedures developed for studying cryogenic fracturing mechanisms in a well environment, and examine typical data indicative of the performance of the system. The experimental apparatus and procedures were specifically designed to conduct cryogenic fracturing tests in specimens under confining stress, with integrated cryogen transport, measurements, and fracture characterization. A true-triaxial loading system was built to simulate reservoir stress levels and anisotropic stress application, and was designed to avoid thermal stresses in tests involving cryogen by arranging discrete components in an open chamber. To maximize thermal shock on the wellbore surface, tubing and wellhead are configured so that liquid nitrogen enters the borehole and flow out after contributing to thermal shock. The temperature at boreholes and along flow lines, borehole pressure, and liquid nitrogen consumption are monitored throughout treatments. Acoustic transmission and pressure-decay measurements are used to characterize fractures before and after the experiments. Breakdown tests performed on un-stimulated specimens and stimulated specimens compare breakdown pressures of the two groups and thus evaluate the performance of thermal fracturing. The laboratory design was able to effectively apply cryogenic stimulations to laboratory rock specimens. The characterization methods were able to capture fracture creation and rock property changes due to cryogenic fracturing.

The rates of mineral dissolution reactions in porous media are difficult to predict, in part because of a lack of understanding of mineral reactive surface area in natural porous media. Common estimates of mineral reactive surface area used in reactive transport models for porous media are typically ad hoc and often based on average grain size, increased to account for surface roughness or decreased by several orders of magnitude to account for reduced surface reactivity of field as opposed to laboratory samples. In this study, accessible mineral surface areas are determined for a sample from the reservoir formation at the Nagaoka pilot CO2 injection site (Japan) using a multi-scale image analysis based on synchrotron X-ray microCT, SEM QEMSCAN, XRD, SANS, and FIB-SEM. This analysis not only accounts for accessibility of mineral surfaces to macro-pores, but also accessibility through connected micro-pores in smectite, the most abundant clay mineral in this sample. While the imaging analysis reveals that most of the micro- and macro-pores are well connected, some pore regions are unconnected and thus inaccessible to fluid flow and diffusion. To evaluate whether mineral accessible surface area accurately reflects reactive surface area a flow-through core experiment is performed and modeled at the continuum scale. The core experiment is performed under conditions replicating the pilot site and the evolution of effluent solutes in the aqueous phase is tracked. Various reactive surface area models are evaluated for their ability to capture the observed effluent chemistry, beginning with parameter values determined as a best fit to a disaggregated sediment experiment (Beckingham et al., 2016) described previously. Simulations that assume that all mineral surfaces are accessible (as in the disaggregated sediment experiment) over-predict the observed mineral reaction rates, suggesting that a reduction of RSA by a factor of 10–20 is required to match the core flood experimental data. While the fit of the effluent chemistry (and inferred mineral dissolution rates) greatly improve when the pore-accessible mineral surface areas are used, it was also necessary to include highly reactive glass phases to match the experimental observations, in agreement with conclusions from the disaggregated sediment experiment. It is hypothesized here that the 10–20 reduction in reactive surface areas based on the limited pore accessibility of reactive phases in core flood experiment may be reasonable for poorly sorted and cemented sediments like those at the Nagaoka site, although this reflects pore rather than larger scale heterogeneity.

Fine-grained sedimentary rocks – namely mudrocks, including their laminated fissile variety — shales – make up about two thirds of all sedimentary rocks in the Earth's crust and a quarter of the continental land mass. Organic-rich shales and mudstones are the source rocks and reservoirs for conventional and unconventional hydrocarbon resources. Mudrocks are relied upon as natural barriers for geological carbon storage and nuclear waste disposal. Consideration of mudrock multi-scale physics and multi-scale spatial and temporal behavior is vital to address emergent phenomena in shale formations perturbed by engineering activities. Unique physical characteristics of shales arise as a result of their layered and highly heterogeneous and anisotropic nature, low permeability fabric, compositional complexity, and nano-scale confined chemical environments. Barriers of lexicon among geoscientists and engineers impede the development and use of conceptual models for the coupled thermal-hydraulic-mechanical-chemical-biological (THMCB) processes in mudrock formations. This manuscript reviews the THMCB process couplings, resulting emergent behavior, and key modeling approaches. We identify future research priorities, in particular fundamental knowledge gaps in understanding the phase behavior under nano-scale confinement, coupled chemo-mechanical effects on fractures, the interplay between physical and chemical processes and their rates, and issues of non-linearity and heterogeneity. We develop recommendations for future research and integrating multi-disciplinary conceptual models for the coupled multi-scale multi-physics behavior of mudrocks. Consistent conceptual models across disciplines are essential for predicting emergent processes in the subsurface, such as self-focusing of flow, time-dependent deformation (creep), fracture network development, and wellbore stability.

Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and a prolonged depletion of residual scCO2. In this study, pore-scale scCO2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO2 into the sandstone-analogue pore network initially saturated by water without dissolved CO2 (dsCO2). During the experiments, time-lapse images of dye intensity, reflecting water pH, were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO2 dissolution and phase equilibrium occurs when scCO2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO2 dissolution at phase interfaces and diffusion of dsCO2 at the pore scale (10–100 µm) observed after scCO2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. This finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO2 at narrow pore throats.

The few prethaw observations of tundra carbon fluxes suggest that there may be large spring releases, but little is known about the scale and underlying mechanisms of this phenomenon. To address these questions, we combined ecosystem eddy flux measurements from two towers near Barrow, Alaska, with mechanistic soil-core thawing experiment. During a 2 week period prior to snowmelt in 2014, large fluxes were measured, reducing net summer uptake of CO2 by 46% and adding 6% to cumulative CH4 emissions. Emission pulses were linked to unique rain-on-snow events enhancing soil cracking. Controlled laboratory experiment revealed that as surface ice thaws, an immediate, large pulse of trapped gases is emitted. These results suggest that the Arctic CO2 and CH4 spring pulse is a delayed release of biogenic gas production from the previous fall and that the pulse can be large enough to offset a significant fraction of the moderate Arctic tundra carbon sink.

We conducted laboratory studies on the geophysical signals from Arctic saline permafrost soils to help understand the physical and mechanical processes during freeze-thaw cycles. Our results revealed low electrical resistivity (

During the past two decades, hydraulic fracturing has significantly improved oil and gas production from shale and tight sandstone reservoirs in the United States and elsewhere. Considering formation damage, water consumption, and environmental impacts associated with water-based fracturing fluids, efforts have been devoted to developing waterless fracturing technologies because of their potential to alleviate these issues. Herein, key theories and features of waterless fracturing technologies, including Oil-based and CO2 energized oil fracturing, explosive and propellant fracturing, gelled LPG and alcohol fracturing, gas fracturing, CO2 fracturing, and cryogenic fracturing, are reviewed. We then experimentally elaborate on the efficacy of liquid nitrogen in enhancing fracture initiation and propagation in concrete samples, and shale and sandstone reservoir rocks. In our laboratory study, cryogenic fractures generated were qualitatively and quantitatively characterized by pressure decay tests, acoustic measurements, gas fracturing, and CT scans. The capacity and applicability of cryogenic fracturing using liquid nitrogen are demonstrated and examined. By properly formulating the technical procedures for field implementation, cryogenic fracturing using liquid nitrogen could be an advantageous option for fracturing unconventional reservoirs.

In modeling of geological carbon storage, dissolution of supercritical CO2 (scCO2) is often assumed to be instantaneous with equilibrium phase partitioning. In contrast, recent core-scale imbibition experiments have shown a prolonged depletion of residual scCO2 by dissolution, implying a non-equilibrium mechanism. In this study, eight pore-scale scCO2 dissolution experiments in a 2D heterogeneous, sandstone-analog micromodel were conducted at supercritical conditions (9MPa and 40°C). The micromodel was first saturated with deionized (DI) water and drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at constant flow rates after scCO2 drainage was completed. High resolution time-lapse images of scCO2 and water distributions were obtained during imbibition and dissolution, aided by a scCO2-soluble fluorescent dye introduced with scCO2 during drainage. These images were used to estimate scCO2 saturations and scCO2 depletion rates. Experimental results show that (1) a time-independent, varying number of water-flow channels are created during imbibition and later dominant dissolution by the random nature of water flow at the micromodel inlet, and (2) a time-dependent number of water-flow channels are created by coupled imbibition and dissolution following completion of dominant imbibition. The number of water-flow paths, constant or transient in nature, greatly affects the overall depletion rate of scCO2 by dissolution. The average mass fraction of dissolved CO2 (dsCO2) in water effluent varies from 0.38% to 2.72% of CO2 solubility, indicating non-equilibrium scCO2 dissolution in the millimeter-scale pore network. In general, the transient depletion rate decreases as trapped, discontinuous scCO2 bubbles and clusters within water-flow paths dissolve, then remains low with dissolution of large bypassed scCO2 clusters at their interfaces with longitudinal water flow, and finally increases with coupled transverse water flow and enhanced dissolution of large scCO2 clusters. The three stages of scCO2 depletion, common to experiments with time-independent water-flow paths, are revealed by zoom-in image analysis of individual scCO2 bubbles and clusters. The measured relative permeability of water, affected by scCO2 dissolution and bi-modal permeability, shows a non-monotonic dependence on saturation. The results for experiments with different injection rates imply that the non-equilibrium nature of scCO2 dissolution becomes less important when water flow is relatively low and the time scale for dissolution is large, and more pronounced when heterogeneity is strong.

CO2 flooding is a proven enhanced oil recovery (EOR) technique and is also considered as a potential method for CO2 sequestration. Despite having successful field trials on CO2 EOR, the effects of reservoir wettability and permeability heterogeneity on the efficiency of miscible CO2 flooding are not well understood. In this work, laboratory investigations have been carried out to evaluate the influence of these properties on the miscible CO2 EOR performance. The wettability of hydrophilic Berea core samples was altered to be oil-wet by vacuum saturation of the clean and dry core samples with n-hexadecane. The permeability heterogeneity was obtained by combining two half pieces of axially split water-wet core samples of different permeabilities. Core flooding experiments were conducted for n-hexadecane - synthetic brine - CO2 systems at 1400 psig backpressure to achieve minimum miscibility pressure (MMP) of CO2 in n-hexadecane at the test temperature (24 ± 1 °C). It was found that wettability strongly influences CO2 EOR. For the alternate cases of previously brine flooded (to remaining oil saturation) oil-wet and water-wet core samples, five pore volumes (PVs) of CO2 recovered 100% and only 43% of remaining oil in place (ROIP) respectively. Three PVs of CO2 could recover only about 0-5% ROIP from the split core samples. The mechanisms underlying these results are discussed. This study sheds light on the significant influence of reservoir wettability and permeability heterogeneity on the performance of miscible CO2 EOR.

Shallow permafrost distribution and characteristics are important for predicting ecosystem feedbacks to a changing climate over decadal to century timescales because they can drive active layer deepening and land surface deformation, which in turn can significantly affect hydrologic and biogeochemical responses, including greenhouse gas dynamics. As part of the U.S. Department of Energy Next-Generation Ecosystem Experiments- Arctic, we have investigated shallow Arctic permafrost characteristics at a site in Barrow, Alaska, with the objective of improving our understanding of the spatial distribution of shallow permafrost, its associated properties, and its links with landscape microtopography. To meet this objective, we have acquired and integrated a variety of information, including electric resistance tomography data, frequency-domain electromagnetic induction data, laboratory core analysis, petrophysical studies, high-resolution digital surface models, and color mosaics inferred from kite-based landscape imaging. The results of our study provide a comprehensive and high-resolution examination of the distribution and nature of shallow permafrost in the Arctic tundra, including the estimation of ice content, porosity, and salinity. Among other results, porosity in the top 2 m varied between 85% (besides ice wedges) and 40%, and was negatively correlated with fluid salinity. Salinity directly influenced ice and unfrozen water content and indirectly influenced the soil organic matter content. A relatively continuous but depth-variable increase in salinity led to a partially unfrozen saline layer (cryopeg) located below the top of the permafrost. The cryopeg environment could lead to year-round microbial production of greenhouse gases. Results also indicated a covariability between topography and permafrost characteristics including icewedge and salinity distribution. In addition to providing insight about the Arctic ecosystem, through integration of lab-based petrophysical results with field data, this study also quantified the key controls on electric resistivity at this Arctic permafrost site, including salinity, porosity, water content, ice content, soil organic matter content, and lithologic properties.

Injection of carbon dioxide (CO2) into saline aquifers confined by low- permeability cap rock will result in a layer of CO2 overlying the brine. Dissolution of CO2 into the brine increases the brine density, resulting in an unstable situation in which more-dense brine overlies less-dense brine. This gravitational instability could give rise to density-driven convection of the fluid, which is a favorable process of practical interest for CO2 storage security because it accelerates the transfer of buoyant CO2 into the aqueous phase, where it is no longer subject to an upward buoyant drive. Laboratory flow visualization tests in transparent Hele-Shaw cells have been performed to elucidate the processes and rates of this CO2 solute-driven convection (CSC). Upon introduction of CO2 into the system, a layer of CO2-laden brine forms at the CO2-water interface. Subsequently, small convective fingers form, which coalesce, broaden, and penetrate into the test cell. Images and time-series data of finger lengths and wavelengths are presented. Observed CO2 uptake of the convection system indicates that the CO2 dissolution rate is approximately constant for each test and is far greater than expected for a diffusion-only scenario. Numerical simulations of our system show good agreement with the experiments for onset time of convection and advancement of convective fingers. There are differences as well, the most prominent being the absence of cell-scale convection in the numerical simulations. This cell-scale convection observed in the experiments may be an artifact of a small temperature gradient induced by the cell illumination.

Measurements on hydrate-bearing laboratory and field samples are necessary in order to provide realistic bounds on parameters used in numerically modeling the production of natural gas from hydrate-bearing reservoirs. The needed parameters include thermal conductivity, permeability, relative permeability-saturation(s) relationships, and capillary pressure-saturation(s) relationships. We have developed a technique to make hydrate-bearing samples ranging in scale from core plug-size to core-size in the laboratory to facilitate making these measurements. In addition to pressure and temperature measurements, we use x-ray computed tomography scanning to provide high-resolution data providing insights on processes occurring in our samples. Several methods are available to make gas hydrates in the laboratory, and we expect that the method used to make the hydrate will impact the behavior of the hydrate sample, and the parameters measured.

Large-scale water injection at The Geysers, California, has generated substantial benefits in terms of sustaining reservoir pressures and production rates, as well as improving steam composition by reducing the content of non-condensible gases (NCGs). Two effects have been recognized and discussed in the literature as contributing to improved steam composition, (1) boiling of injectate provides a source of "clean" steam to production wells, and (2) pressurization effects induced by boiling of injected water reduce upflow of native steam with large NCG concentrations from depth. In this paper we focus on a possible additional effect that could reduce NCGs in produced steam by dissolution in a condensed aqueous phase.Boiling of injectate causes pressurization effects that will fairly rapidly migrate outward, away from the injection point. Pressure increases will cause an increase in the saturation of condensed phase due to vapor adsorption on mineral surfaces, and capillary condensation in small pores. NCGs will dissolve in the additional condensed phase which, depending upon their solubility, may reduce NCG concentrations in residual steam.We have analyzed the partitioning of HCl between vapor and aqueous phases, and have performed numerical simulations of injection into superheated vapor zones. Our simulations provide evidence that dissolution in the condensed phase can indeed reduce NCG concentrations in produced steam.

The Gardner and van Genuchten models of relative permeability are widely used in analytical and numerical solutions to flow problems. However, the applicab ility of the Gardner model to real problems is usually limited, because empirical relative permeability data to calibrate the model are not routinely available. In contrast, van Genuchten parameters can be estimated using more routinely available matric potential and saturation data. However, the van Genuchten model is not amenable to analytical solutions. In this paper, we introduce generalized conversion formulae that reconcile these two models. In general, we find that the Gardner parameter alpha G is related to the van Genuchten parameters alpha vG and n by alpha G=alpha vG ~; 1:3 n. This conversion rule will allow direct recasting of Gardner-based analytical solutions in the van Genuchten parameter space. The validity of the proposed formulae was tested by comparing the predicted relative permeability of various porous media with measured values.

We discuss the range of activities at Lawrence Berkeley National Laboratory in support of gas production from natural hydrates. Investigations of production from the various classes of hydrate deposits by numerical simulation indicate their significant promise as potential energy sources. Laboratory studies are coordinated with the numerical studies and are designed to address knowledge gaps that are important to the prediction of gas production. Our involvement in field tests is also briefly discussed.