Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Wetting mechanism and alteration of nano-sized shale pores: Insights from contrast variation small angle neutron scattering

Wettability of tight shale is crucial for fluid flow and mass transport process in energy geosciences. However, understanding the interfacial chemistry and wetting mechanisms at sub-nano-pore scales remains a formidable challenge. Here, in this study, the Contrast Variation technique of Small Angle Neutron Scattering (CV-SANS) is employed to investigate shale’s interfacial chemistry using reagents that possess a range of different polarities, including water, n-decane, toluene, and dimethyl methanamide. Through five different experimental strategies, we have demonstrated a successful modification of shale wettability, ranging from enhancement, weakening, to reversal. Delving into the mechanisms, we illustrated the crucial role of pre-existing liquid films in these changes, where the uniquely co-existing polar and non-polar functional groups in dimethyl methanamide acted as a conduit for interfacial chemistry adjustments. Furthermore, a solvent immersion led to matrix dilation as well as liberation of residual oil-occupied pores, resulting in altered pore size distributions, with hydrogen bonding playing a significant role in the polar groups. Interestingly, despite shale exhibiting a stronger affinity for oil over water, hydrophilic solvents induced more substantial dilation than lipophilic ones. Collectively, this work elucidates the dynamic change of interfacial chemistry via the configuration of polarity using chemical reagents, and the CV-SANS technique underscores its invaluable utilities in decoding the interfacial wettability traits in nanopore space of shale.

58 GEOSCIENCES↗
DOE OSTI2024

Initial Laboratory Measurements Probing Hydrogen Interactions with Eagle Ford Shale and Pyrite: Potential Implications for Subsurface Hydrogen Storage

Hydrogen (H 2 ) has the potential to be a transformative technology as an enabler to a low-carbon future and promoter of renewable energy. When H 2 is injected and stored in the subsurface, it has the potential to interact with the caprock (usually shale) which overlies and seals the storage reservoir. Here, this study examines geochemical reactions or changes in surface morphology to Eagle Ford Shale, a proxy for caprock, upon exposure to H 2 at 50°C and 10.3 MPa. Reactions were also performed with N 2 to provide an experimental control. Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS), Atomic Force Microscopy (AFM), and Optical Photothermal-Infrared (O-PTIR) Spectroscopy were applied to quantify changes on the microscale and nanoscale level. Fluid chemistry changes were monitored with Ion Chromatography (IC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Exposure of Eagle Ford Shale to H 2 gas alone did not result in any alterations to the shale chemically or any changes in the surface morphology. Exposure of Eagle Ford Shale to both H 2 and water as well as N 2 and water resulted in changes to the surface morphology because of gypsum dissolution and reprecipitation, thus indicating that H 2 is not necessary to promote changes. Pure pyrite was the most reactive with H 2 possibly resulting in a reduction to pyrrhotite. These initial studies suggest that the extent of reactions activated by hydrogen with caprock are minor under the temperature and pressure conditions that would represent underground hydrogen storage.

08 HYDROGEN↗
DOE OSTI2025

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN↗
DOE OSTI2024

Climate forcing controls on carbon terrestrial fluxes during shale weathering

Climate influences near-surface biogeochemical processes and thereby determines the partitioning of carbon dioxide (CO 2 ) in shale, and yet the controls on carbon (C) weathering fluxes remain poorly constrained. Using a dataset that characterizes biogeochemical responses to climate forcing in shale regolith, we implement a numerical model that describes the effects of water infiltration events, gas exchange, and temperature fluctuations on soil respiration and mineral weathering at a seasonal timescale. Our modeling approach allows us to quantitatively disentangle the controls of transient climate forcing and biogeochemical mechanisms on C partitioning. We find that ~3% of soil CO 2 (1.02 mol C/m 2 /y) is exported to the subsurface during large infiltration events. Here, net atmospheric CO 2 drawdown primarily occurs during spring snowmelt, governs the aqueous C exports (61%), and exceeds the CO 2 flux generated by pyrite and petrogenic organic matter oxidation (~0.2 mol C/m 2 /y). We show that shale CO 2 consumption results from the temporal coupling between soil microbial respiration and carbonate weathering. This coupling is driven by the impacts of hydrologic fluctuations on fresh organic matter availability and CO 2 transport to the weathering front. Diffusion-limited transport of gases under transient hydrological conditions exerts an important control on CO 2 (g) egress patterns and thus must be considered when inferring soil CO 2 drawdown from the gas phase composition. Our findings emphasize the importance of seasonal climate forcing in shaping the net contribution of shale weathering to terrestrial C fluxes and suggest that warmer conditions could reduce the potential for shale weathering to act as a CO 2 sink.

58 GEOSCIENCES↗
DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗
DOE OSTI2021

Multiscale Characterization of the Caney Shale — An Emerging Play in Oklahoma

From a hydrocarbon perspective, the Caney Shale has historically been evaluated as a sealing unit, which resulted in limited studies characterizing the rock properties of the Caney Shale and its suitability for hydraulic fracturing. The objective of our research is to help bridge the current knowledge gap through the integration of multiscale laboratory techniques and to characterize the macro- and microscale rock properties of the Caney Shale. We employed an integrated approach for the characterization of the Caney using 200 ft (61 m) of Caney core from a target well in southern Oklahoma. Core observation and petrographic analysis of thin sections were combined to characterize the general rock types and associated fabrics and textures. Mineralogical composition, pore system architecture, and rock fabric were analyzed using x-ray diffraction (XRD), scanning electron microscopy/energy dispersive x-ray spectroscopy (SEM/EDS), and focused ion beam (FIB)-SEM. In addition, rebound hardness and indentation testing were carried out to determine rock hardness (brittleness) and elasticity, respectively. With the integrated multiscale characterization, three mixed carbonate-siliciclastic rock types were identified — mudstone, calcareous siltstone, and silty carbonate — likely representing a spectrum of deposition from low to relatively high energy environments in the distal portions of a ramp system. Silty carbonate contains mostly interparticle pores. The calcareous siltstones and silty mudstones contain a combination of organic matter pores and interparticle pores. Each of the rock types shows unique mineralogical compositions based on XRD. The mudstone lithofacies has the highest clay content and the least carbonate content. Calcareous siltstones show moderate carbonate and clay content. Silty carbonate indicates the highest carbonate content with the least clay content. In an order of mudstone, calcareous siltstone, and silty carbonate, rebound hardness and Young’s modulus show an increasing trend. As a result of rock-fluid interactions, there are potential scaling reactions during completion and production that could ultimately affect permeability and production rates. Overall, the proposed multiscale integration approach is critical for the geologic characterization of most rocks. However, in shale reservoirs dominated by microporosity and microstructure where engineered fractures are expected to provide permeability at a reservoir scale, successful integration is essential. An optimized, integrated geological characterization of the Caney Shale that is well aligned with the engineering designs in drilling, completing, and producing wellbores will ultimately lead to optimal production while providing safe and environmentally responsible operations.

Wang, Yulun↗
DOE OSTI2022

AmeriFlux US-SSH Susquehanna Shale Hills Critical Zone Observatory

This is the AmeriFlux version of the carbon flux data for the site US-SSH Susquehanna Shale Hills Critical Zone Observatory. Site Description - The Susquehanna Shale Hills Critical Zone Observatory is comprised of one first-order catchment in the Susquehanna River basin. This catchment, known as Shale Hills, is about 8 hectares in total area. The stream that defines the Shale Hills catchment flows into Shavers Creek in the Juniata River sub-basin. The vegetation cover at Shale Hills is dominated by deciduous broadleaf forest, with some evergreen needleleaf trees along the stream.

Forsythe, Brandon R.↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Pore Structure and Wettability of Bossier Shale, East Texas, United States: Insights from Integrated Porosimetry, Scattering, and Imbibition Approaches

In this study, ascertaining the pore geometry and wettability characteristics of tight shales is of great significance for revealing the mechanisms of occurrence, migration, and production of shale resources. Taking the Bossier Shale collected from Well A, East Texas, USA, as an example, the pore geometry and wettability behavior in a broad nm-& μm scale pore spectrum were quantified based on integrated techniques, including water immersion porosimetry (WIP), mercury intrusion porosimetry(MIP), (ultra)small-angle X-ray scattering [(U)SAXS], contact angle, and liquid spontaneous imbibition (SI). Mainly owing to differences in the sample sizes used, data interpretation, as well as the detectable pore type (e.g., connected or not) and pore diameter ranges among different methods, sample porosities derived from WIP, MIP, and (U)SAXS were quite different, which range from 7.70 to 13.34%, 5.12 to 9.52%, and 2.48 to 6.44%, respectively. Therefore, a comparison of the porosities derived from different methods must be handled with utmost care. Additionally, although (U)SAXS is an effective technique for detecting both connected and non-connected pores, the fraction of non-connected pores could be underestimated by directly comparing the pore size distribution derived from (U)SAXS and MIP tests. This is because the sample sizes used for and pore information reflected by these two methods are different. Furthermore, the studied Bossier Shales exhibited various contact angle values and imbibition behaviors when using differently polarized liquids, revealing their mixed-wet characteristics. Comprehensively considering the differences in contact angle, imbibition slope, and imbibed liquid volume, three sub-categories of wettability behavior, with respect to more oil-wet, more water-wet, and intermediate mixed-wet, were qualitatively identified.

(U)SAXS↗
DOE OSTI2021

Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery

One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.

03 NATURAL GAS↗
DOE OSTI2023

Experimental flow-through a single fracture with monolayer proppant at reservoir conditions: A case study on Caney Shale, Southwest Oklahoma, USA

Global energy systems are undergoing a crucial transition to clean energy sources, thereby moving away from fossil-fuel based energy with the ultimate aim of reducing carbon dioxide emissions. It is essential to continue to advance in exploring and developing shale gas resources as this will underpin the global supply of natural gas that is needed to support the energy transition. The challenge hereby is that the slow fluid flow of the shale matrix due to its very low permeability requires large-scale volume fracturing to produce a conductive fracture network that enhances the connectivity between the shale formation and the wellbore. This work examines the influence of a thin proppant layer on a single fracture using a flow-through experiment, and also explores whether proppant embedment in the Caney shale is affected by the rock mineralogy, surface roughness, fluids, confining stress, time, temperature and bedding. The experiment was conducted for a duration of nine days (216 h), experimental temperature was varied from ambient temperature to reservoir temperature of 125 °C (257 °F) and confining stress was varied from 367 psi (2.53MPa) to a maximum of 4011.82 psi (27.66 MPa). Here, we find that the conductivity of the fracture is primarily influenced by the layer of proppant used, surface roughness, mineralogy, fluids, temperature and closure stress.

Caney shale↗
DOE OSTI2022

Kinetic and equilibrium reactions on natural and laboratory generation of thermogenic gases from Type II marine shale

The phenomenon that laboratory pyrolysis experiments produce much wetter gases than those in natural reservoirs is a long-recognized and debated problem in the investigation of natural gases in sedimentary basins. In this study, we explore the discrepancy by pyrolyzing a type II kerogen from the Woodford Shale in Oklahoma, compared with the previous results on the produced natural gases from the Arkoma Basin generated from the same source rock (Liu et al., 2019) with the discussion of gas and isotopic compositions at bulk and position-specific (PS) levels. An improved GC-pyrolysis-GC IRMS method is applied for the determination of PS δ 13 C of propane produced in the pyrolysis of the Woodford Shale at Easy %R o from 0.76 to 3.27. Kinetic and thermodynamic considerations of the chemical and isotopic compositions of the natural and laboratory pyrolysis gases suggest that the generation of light hydrocarbons involves uni-directional cracking reactions, exchange reactions with water, and likely reversible reactions among light hydrocarbons and other H-containing volatiles. After the gas generation in the unconventional Woodford Shale reservoirs, the C 1 -C 4 gases might have approached close to chemical equilibrium of C 1 -C 3 and isotope equilibrium of C 2 -C 1 and C 3 -C 1 pairs at their peak temperatures. The capping H for the generation of C 1 -C 4 in the Woodford Shale gases appears to have experienced at least partial exchange with the water, while that in the pyrolysis gases is only originated from organic-bound compounds with large kinetic isotope effects (KIE). Our findings indicate that elevated compound-specific and PS δ 13 C values of propane in the wet-gas cracking stage are significantly influenced by the breakdown of the thermally stable compounds (e.g., remaining kerogen, residues). A first synthesis of PS δ 13 C and δ 2 H isotopic compositions of propane from this study and the literature data suggests relatively similar isotopic structures of propane precursors in kerogens. Finally, this study demonstrates that PS isotope analysis of propane can contribute to identifying various geological (e.g., maturation, wet-gas cracking, H exchange, diffusion) and biodegradation processes.

58 GEOSCIENCES↗
DOE OSTI2023

Microwave-enhanced methane cracking for clean hydrogen production in shale rocks

Here, steam methane reforming (SMR) generates about 95% of hydrogen (H 2 ) in the U.S. using natural gas as a main feedstock. However, this technology also generates a large amount of carbon dioxide (CO 2 ), a major greenhouse gas causing global warming. Carbon capture and storage (CCS) technique is required, but the cost and safety of storing CO 2 underground are a concern. Here we propose a new approach using microwave/electromagnetic irradiation to produce clean hydrogen from unrecovered hydrocarbons within petroleum reservoirs. Solid carbon or CO 2 produced during this process will be simultaneously sequestrated underground without involving CCS. In this paper, we perform a series of experiments to investigate the in-situ hydrogen production from shale gas (methane) conversion by passing a methane stream through a packed shale rock sample heated by microwave. We found that methane conversion was significantly enhanced in the presence of Fe and Fe 3 O 4 particles as catalysts, with a conversion of 40.5% and 100% at reaction temperature of 500 °C and 600 °C, respectively. Methane conversion is promoted at a lower reaction temperature by the catalytic effect of minerals in shale. Additionally, the influences of catalysts, shale rock, and methane flow rate are characterized.

08 HYDROGEN↗
DOE OSTI2022

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Shale gas: Geochemical and physical constraints on genesis, storage, and producibility

Our research project on shale gas included two tasks: (1) Testing of different shales’ ability to geocatalytically generate natural gas at low temperatures; verification of geocatalytic methanogenesis would potentially render shale gas a quasi-renewable resource, and (2) quantification and characterization of porosity and permeability in organic-rich rocks with implication for producibility. Adequate porosity is required for geological storage of shale gas. At the same time, pore connectivity is essential for gas producibility.

03 NATURAL GAS↗
DOE OSTI2024

Electrofracturing of Shale at Elevated Pressure

Electrofracturing deeply buried shale formations could be used to increase reservoir permeability and improve reservoir production without requiring large volumes of freshwater. This paper describes a novel experimental system and initial test results to electrofracture shale under high confining pressures. Core-scale laboratory testing was performed on twelve rock samples recovered from a shale gas reservoir. Each sample was subjected to confining pressures of 20.7 MPa (3000 psi) or 58.6 MPa (8000 psi), representative of overburden pressures at depth. Samples were then subjected to application of high voltage until specimen fracture. The experiments produced deformed samples with multiple fracture types, both parallel and oblique to bedding planes. Electrofracturing increased permeabilities by up to nine orders of magnitude for extended time periods. Rock fracture and throughgoing fractures were demonstrated. Computed tomography images revealed the creation of fractures and tube/tunnel flow channels, which resisted closure under hydrostatic pressures up to 58.6 MPa. The breakdown energy and permeability changes in the sample were independent of applied confining pressure. The cumulative energy input required for fracture depended on applied confining pressure and sample length. The energy required to fracture samples up to 9 cm in length is generally more than 0.5 kJ/cm, but no greater than 1 kJ/cm. Our results show that electrofracture of shales under confining pressure is possible and could be a possible water-free mechanism for reservoir stimulation.

04 OIL SHALES AND TAR SANDS↗