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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Probing Pore Size and Connectivity in Porous Silicas Using 13C MAS NMR Spectroscopy of Supercritical Methane

Measuring the pore size and pore-size distributions and exploring the fluid-exchange dynamics between different types of pores in porous materials remains a significant experimental challenge but is critical to understanding catalysis, chromatography, nutrient cycling, and a whole range of geochemical phenomena, including shale gas and tight gas extraction. Here, we present the results of 1D 13 C NMR and 2D exchange spectroscopy (EXSY) NMR investigations of a porous silica using supercritical methane (scCH 4 ) as a direct probe of pore size and fluid exchange between pore types. The results show that the 13 C chemical shift of scCH 4 adsorbed in nanometer-scale silica pores becomes more negative with increasing pore diameter, in agreement with trends reported for gas hydrates, zeolites, MOFs, and clays and other microporous (<10 nm pores) geochemical materials. These 13 C chemical shifts follow a natural log–linear trend with pore size in the vacuum-dried porous silicas studied here, allowing one to predict pore size based on the 13 C chemical shift in dry silica nanopores. The EXSY and 1D 13 C NMR results both show that CH 4 exchanges between pore and bulk fluid environments over rate scales from 0.01 to 1 kHz in vacuum-dried samples and that exchange dynamics must be considered when interpreting 13 C NMR of pore-adsorbed CH 4 . Introducing H 2 O to the system causes the 13 C chemical shift to become more negative with increasing H 2 O content as a result of H 2 O preferentially filling small pores, forcing CH 4 to occupy larger pores that are better connected to the bulk environment. Likewise, the 1D 13 C data show a decreased distribution of exchange rates between pore-adsorbed and bulk CH 4 in the presence of H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Effect of Confinement on the Density and Diffusivity of Organic Molecules in Single-Digit Nanopores Relative to Bulk Fluid Behavior

Prediction of fluid properties in microporous and mesoporous confinement is non-trivial due to distinct influences from pore surface functionality and shape and fluid molecular characteristics. The vast majority of literature focuses nearly exclusively on a small number of components that are generally low molecular weight gases and vapors. Missing are studies that utilize a relatively large number of compounds that are liquid at ambient conditions to probe the impact of molecular shape and size on confined fluid properties. Such materials are of great interest in the extraction of oil and gas from shale reservoirs or for the design of selective adsorbents for chemical separation. In this work, we reveal non-intuitive relationships that emerge when comparing bulk fluid properties to those when confined in amorphous silica single-digit nanopores (approx. 5 nm SBA-15), specifically the confined density and self-diffusivity. A counterintuitive result was obtained where larger molecules had relatively faster self-diffusivities; however, this behavior was revealed to be caused by these larger compounds exhibiting a relatively low density when confined. As the prediction of confined fluid properties impacts a diversity of disciplines, the results highlight the need for additional investigations where consistent and large-scale experimental procedures are employed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Electrochemically Engineered, Highly Energy-Efficient Conversion of Ethane to Ethylene and Hydrogen below 550 °C in a Protonic Ceramic Electrochemical Cell

Ethylene is one of the largest building blocks in the petrochemical industry, mainly produced by steam cracking of ethane derived from naphtha or shale gas at high temperatures (>800 °C). Despite its technical maturity and economic competitiveness, the thermal steam cracking of ethane is highly energy-intensive. In this work, an electrochemically engineered direct conversion process of ethane to produce hydrogen and ethylene using a planar protonic ceramic membrane reactor with a bi-functional three-dimensional catalytic electrode is reported, with a single-pass ethane conversion of 40% and ethylene yield of 26.7% at 550 °C. Compared with the industrial ethane steam cracking, this method saves process energy input by 45.1% and improves process energy efficiency by 50.6%, based on comprehensive process simulation using Aspen Plus software. Further, steam electrolysis treatment under the solid oxide electrolysis cell mode can regenerate the system’s catalytic performance and significantly alleviate catalytic degradation by 74%, demonstrating high techno-economic viability.

bi-functional electrode↗
DOE OSTI2022

Membrane Bioreactor Pretreatment of High-Salinity O&G Produced Water

Produced water (PW) from oil and gas production contains variable constituents that are difficult to remove with conventional treatment processes. The focus of this study was to explore the long-term performance of a membrane bioreactor (MBR) for the removal of organic constituents from PW and how performance and microbial community composition are affected by progressively increasing salinity and introduction of PW from different shale basins around the United States. Dissolved organic carbon removal from the PW remained consistent throughout the study, averaging 86% from the Denver-Julesburg basin PW and 66% removal from the Permian basin PW. Surfactant removal was less consistent, showing 87% removal of poly(ethylene glycols) (PEGs) at a total dissolved solid (TDS) concentration of 40 g/L but only 58% removal at a TDS concentration of 100 g/L. Diversity in the microbial community decreased during reactor establishment but increased at TDS concentrations above 80 g/L. Finally, the results of this study suggest that MBRs can be effective PW pretreatment processes even at high salinities.

42 ENGINEERING↗
DOE OSTI2023

Direct Conversion of Ethane to Oxygenates, Ethylene, and Hydrogen in a Noncatalytic Biphasic Plasma Microreactor

We selectively upgrade ethane (C 2 H 6 ) to ethanol (C 2 H 5 OH), methanol (CH 3 OH), and acetic acid (CH 3 COOH) in a catalyst-free, continuous, argon/water biphasic plasma microreactor. The water (H 2 O) evaporates and electron- dissociates into OH· radicals. OH· recombines with alkyl radicals, produced via electron dissociation of ethane, to generate the oxygenates that absorb into H 2 O. A plasma-assisted path, reminiscent of the low-temperature thermocatalytic ethane steam reforming, leads to significant H 2 coproduction. The gaseous stream also comprises CO 2 and C 2 H 4 . Up to 1.3 and 1 μmol min –1 of liquid C 2 H 5 OH and CH 3 OH are attained, respectively. Compared to CO 2 -assisted ethane plasma conversion, which produces many oxygenates with low selectivity, the carbon selectivity can range from >70% C 2 H 5 OH, CH 3 OH, and CH 3 COOH to 60% C 2 H 4 . Furthermore, the low carbon footprint, electrified, modular, intensified process using a reactive evaporation and separation plasma could pave the way for the valorization of underutilized shale gas resources in remote areas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗
DOE OSTI2026

Reacting CO 2 with Light Alkanes to Value-Added Products

Catalytic conversion of anthropogenic carbon dioxide (CO 2 ) into value-added products is a promising strategy to mitigate global carbon emissions. Concurrently, the shale gas revolution has provided an abundant supply of light alkanes (methane, ethane, propane, and butane), presenting a unique opportunity to employ these underutilized hydrocarbons as an effective, low-cost hydrogen source for CO 2 reduction. In this Perspective, we summarize past efforts, current state, and future opportunities for reacting CO 2 with light alkanes to generate a diverse range of value-added products. Compared with direct alkane conversion, the introduction of CO 2 fundamentally alters reaction thermodynamics and kinetics, enabling selective C–H and C–C bond activation while suppressing catalyst deactivation from coke formation. Building on decades of research in dry reforming and CO 2 -assisted dehydrogenation, recent advances in catalyst design have enabled CO 2 -assisted dehydrogenation processes that approach chemical equilibrium for the selective production of olefins and syngas. Importantly, advances in catalyst design and reactor engineering have further expanded the product scope beyond gas-phase (syngas and olefins) to include liquid-phase (oxygenates and aromatics), and solid-phase products (carbon nanomaterials). We highlight key catalyst design principles for controlling reaction pathways and discuss major challenges and opportunities in developing selective and versatile platforms for the simultaneous upgrading of CO 2 and light alkanes.

CO2↗
DOE OSTI2020

Viscous Dissipation and Apparent Permeability of Gas Flow in Nanoporous Media

The nanopore confinement effect can increase the mobility for gas flow in shale rocks, and the apparent permeability is usually adopted to account for this effect. However, most of the apparent permeability calculation models are derived based on straight channels (capillaries) and neglect the influence of end effect. The end effect is caused by the bending of streamlines which yields more viscous dissipation and additional flow resistance for gas flow in nanoporous media. Here, for the first time the viscous dissipation for gas flow in nanoporous media is analyzed. In addition, an improved apparent permeability calculation model is proposed based on the Beskok-Karniadakis (B-K) model by introducing the influence of end effect, which is characterized by the ratio of the length to the width of a short channel (L/H). The accuracy of this model is validated by lattice Boltzmann simulations of gas flow in three different geometries: an infinite-length straight channel, a finite-length straight channel, and nanoporous media. For the infinite-length straight channel, there is no end effect therefore both the B-K model and the improved model can well predict the apparent permeability. However, for the finite-length straight channel and nanoporous media, the original B-K model overestimates the gas apparent permeability as it neglects the influence of end effect, while the improved model can still well predict the gas apparent permeability. In summary, the improved apparent permeability calculation model is more reasonable in predicting gas mobility in nanoporous media.

58 GEOSCIENCES↗
DOE OSTI2021

Injection Parameters That Promote Branching of Hydraulic Cracks

Fluid injection into rock formations can either produce complex branched hydraulic fractures, create simple planar fractures, or be dominated by porous diffusion. Currently, the optimum injection parameters to create branched fractures are unknown. We conducted repeatable hydraulic fracturing experiments using analog-rock samples with controlled heterogeneity to quantify the fluid parameters that promote fracture branching. A large range of injection rates and fluid viscosities were used to investigate their effects on induced fracture patterns. Paired with a simple analytical model, our results identify the threshold at which fracture transitions from an isolated planar crack to branched cracks when closed natural fractures exist. These results demonstrate that this transition can be controlled by injection rate and fluid viscosity. In relation to the field practices, the present model predicts slickwater and lower viscosity fluid injections promote fracture branching, with the Marcellus shale used as an example.

58 GEOSCIENCES↗
DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗
DOE OSTI2024

A Tale of Two Catchments: Causality Analysis and Isotope Systematics Reveal Mountainous Watershed Traits That Regulate the Retention and Release of Nitrogen

Abstract Mountainous watersheds are characterized by variability in functional traits, including vegetation, topography, geology, and geomorphology, which determine nitrogen (N) retention, and release. Coal Creek and East River are two contrasting catchments within the Upper Colorado River Basin that differ markedly in total nitrate (NO 3 − ) export. The East River has a diverse vegetation cover, and sinuous floodplains, and is underlain by N‐rich marine shale. At 0.21 ± 0.14 kg ha −1 yr −1 , the East River exports ∼3.5 times more NO 3 − relative to the conifer‐dominated Coal Creek (0.06 ± 0.02 kg ha −1 yr −1 ). While this can partly be explained by the larger size of the East River, the distinct watershed traits of these two catchments imply different mechanisms controlling the aggregate N‐export signal. A causality analysis shows physical and biogenic processes were critical in determining NO 3 − export from the East River catchment. Stable isotope ratios of NO 3 − (δ 15 N NO3 and δ 18 O NO3 ) show the East River catchment is a strong hotspot for biogeochemical processing of NO 3 − at the hillslope soil‐saprolite. By contrast, the conifer‐dominated Coal Creek retained nearly all atmospherically deposited NO 3 − , and its export was controlled by catchment hydrological traits (i.e., snowmelt periods and water table depth). The conservative N‐cycle within Coal Creek is likely due to the abundance of conifer trees, and smaller riparian regions, retaining more NO 3 − overall and reduced processing prior to export. This study highlights the value of integrating isotope systematics to link watershed functional traits to mechanisms of watershed element retention and release.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Caprock Remains Water Wet Under Geologic CO 2 Storage Conditions

Abstract Carbon storage technology is primarily targeted in saline formations, which is a porous rock matrix filled with brine, sealed with a low permeability caprock. There are significant variations of CO 2 wetting properties, typically reported in the literature as contact angle of CO 2 and brine interacting with a rock material, suggesting that CO 2 could become wetting under geostorage conditions and negatively impact containment effectiveness. Here, we performed the first controlled laboratory measurements of CO 2 ‐brine contact angles on shale rocks from low permeability sealing formations with distinctive mineralogic properties—calcite‐rich, quartz‐rich, and dolomite‐rich. We targeted temperatures at 40° and 100°C, pressures at 8.3, 34.5, and 62.1 MPa, and salinity at 35,000 and 260,000 ppm. Results show no significant change in contact angle with mineralogy, temperature, pressure, salinity, and CO 2 bubble size. We conclude that caprocks will remain water‐wet at geologic CO 2 storage conditions and will not lose their capillary sealing capacity.

Tapriyal, Deepak↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2025

Deep Roots Supply Reactivity and Enhance Silicate Weathering in the Bedrock Vadose Zone

In upland environments, roots commonly extend deep below soil into partially saturated bedrock. This Bedrock Vadose Zone (BVZ) has been shown to store and circulate water, host organic carbon respiration and serve as a critical source of rock-derived nutrients. However, the extent to which deep roots influence chemical weathering rates remains poorly understood. Here, we report 4 years of depth-resolved major ion chemistry over a 16-m thick BVZ hosting a deep rhizosphere in a catchment subject to a Mediterranean climate. These data allow development and validation of a reactive transport model (RTM), revealing that the timescales of water storage and drainage in the BVZ are sufficient to facilitate substantial chemical weathering of the shale bedrock. However, observed solute concentrations are only reproduced by the RTM when we explicitly include measured rates of CO 2(g) production meters below soil driven by the deeply rooted forest. By combining direct observations and a process-based RTM we conclude that the carbon respiration promoted by deep roots significantly enhances chemical weathering rates in the BVZ, constituting 43% ± 3% of total solute flux from the base of the BVZ to the water table.

Osorio‐Leon, Ivan D. [University of Illinois Urban↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2020

Carbon footprint of global natural gas supplies to China

As natural gas demand surges in China, driven by the coal-to-gas switching policy, widespread attention is focused on its impacts on global gas supply-demand rebalance and greenhouse gas (GHG) emissions. Here, for the first time, we estimate well-to-city-gate GHG emissions of gas supplies for China, based on analyses of field-specific characteristics of 104 fields in 15 countries. Results show GHG intensities of supplies from 104 fields vary from 6.2 to 43.3 g CO 2 eq MJ -1 . Due to the increase of GHG-intensive gas supplies from Russia, Central Asia, and domestic shale gas fields, the supply-energy-weighted average GHG intensity is projected to increase from 21.7 in 2016 to 23.3 CO 2 eq MJ -1 in 2030, and total well-to-city-gate emissions of gas supplies are estimated to grow by ~3 times. While securing gas supply is a top priority for the Chinese government, decreasing GHG intensity should be considered in meeting its commitment to emission reductions.

03 NATURAL GAS↗