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 361 records · Page 20

Completion design improvement using a deep convolutional network

Maximizing stimulated natural and hydraulic fracture network is one of the primary hydraulic fracturing concerns for economic production from a horizontal shale gas well. Geomechanical facies and preexisting fractures in each stage are identified based on similarities in formation characteristics to optimize the locations of perforation clusters. This often requires analyzing large volumes of drilling, Logging While Drilling (LWD) and Measurement While Drilling (MWD) data. In this paper, we develop a methodology that calculates the mechanical specific energy (MSE) using real-time drill string acceleration signals directly from its definition. High resolution vibration signals have been collected using a tri-axial accerlometer, which was an auxiliary tool included in acoustic borehole imager. This technique provides a cost-efficient solution for engineered completion design. Furthermore, we adopt deep Convolutional Neural Network (CNN) with signal processing to build a data pipeline that effectively extracts patterns from dynamic acceleration signals for rock lateral MSE classification. First, we apply discrete wavelet transform and Short-Time Fourier Transform (STFT) for signal denoising and pattern recognition. Then we construct an image dataset using multi-scale image fusion at pixel level from 3 sensor channels, including axial, lateral acceleration spectrograms and zero-padded revolutions per minute (RPM). The resulted RGB image dataset includes 4,000 images of 5 MSE ranges with various rock strength conditions. Our results demonstrate that the proposed deep learning model can achieve more than 90% classification accuracy. The deep learning results, as a reference source, were applied in selected Marcellus Shale Energy and Environmental Lab (MSEEL) wells engineered completion located in the Marcellus shale gas site.

03 NATURAL GAS↗

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗

Computed Tomography Scanning and Geophysical Measurements of UW Enterprises LP 1-250512-129 Well in Southwestern Indiana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Illinois Basin core of the Upper Devonian-Early Mississippian New Albany Shale Formation from Posey County, Indiana. The primary impetus of this work is a collaboration between Indiana Geological and Water Survey at Indiana University Bloomington, NETL, and Woolsey Operating Company LLC to characterize and make publicly available core information from the New Albany Shale of the Illinois Basin. Core characterization of this unconventional oil/gas well will aid in understanding the lithology changes and the fracture complexity of the New Albany Shale. There is a potential for this formation to be developed in the future for critical mineral and rare earth element (CM/REE) extraction. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/uw-enterprises. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on this core. None of the equipment used was suitable for direct visualization of the pore space in the fine-grained structures studied; however, fractures, discontinuities, and millimeter-scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Targeted higher resolution CT scanning of select sections was performed with NETL’s industrial and micro-CT scanner. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗

Source characterization of volatile organic compounds at Carlsbad Caverns National Park

Carlsbad Caverns National Park (CAVE), located in southeastern New Mexico, experiences elevated ground-level ozone (O 3 ) exceeding the National Ambient Air Quality Standard (NAAQS) of 70 ppbv. It is situated adjacent to the Permian Basin, one of the largest oil and gas (O&G) producing regions in the US. In 2019, the Carlsbad Caverns Air Quality Study (CarCavAQS) was conducted to examine impacts of different sources on ozone precursors, including nitrogen oxides (NO x ) and volatile organic compounds (VOCs). Here, we use positive matrix factorization (PMF) analysis of speciated VOCs to characterize VOC sources at CAVE during the study. Seven factors were identified. Three factors composed largely of alkanes and aromatics with different lifetimes were attributed to O&G development and production activities. VOCs in these factors were typical of those emitted by O&G operations. Associated residence time analyses (RTA) indicated their contributions increased in the park during periods of transport from the Permian Basin. These O&G factors were the largest contributor to VOC reactivity with hydroxyl radicals (62%). Two PMF factors were rich in photochemically generated secondary VOCs; one factor contained species with shorter atmospheric lifetimes and one with species with longer lifetimes. RTA of the secondary factors suggested impacts of O&G emissions from regions farther upwind, such as Eagle Ford Shale and Barnett Shale formations. The last two factors were attributed to alkenes likely emitted from vehicles or other combustion sources in the Permian Basin and regional background VOCs, respectively. Implications: Carlsbad Caverns National Park experiences ground-level ozone exceeding the National Ambient Air Quality Standard. Volatile organic compounds are critical precursors to ozone formation. Measurements in the Park identify oil and gas production and development activities as the major contributors to volatile organic compounds. Emissions from the adjacent Permian Basin contributed to increases in primary species that enhanced local ozone formation. Observations of photochemically generated compounds indicate that ozone was also transported from shale formations and basins farther upwind. Therefore, emission reductions of volatile organic compounds from oil and gas activities are important for mitigating elevated O 3 in the region.

54 ENVIRONMENTAL SCIENCES↗

Characterizing Radium Distributions, Source Attribution and Impacts from Marcellus Produced Water

This presentation characterizes spatial trends in radium (Ra) concentrations, sources, and potential impacts generated from Marcellus Shale produced waters. Here, we combine previously published produced water data with additional Ra activity measurements extracted from residual waste reports submitted to the Pennsylvania Department of Environmental Protection from 2012-2023. This is the largest, spatially referenced compilation of Ra activity measurements from produced water sourced from the Marcellus Shale to date. The median Marcellus Shale produced water Ra activity presented here is 6,906 pCi/L (interquartile range = 12,515 pCi/L). Radium activities in the positively skewed dataset can exceed 30,000 pCi/L in some parts of the state. The regional trends in Ra concentrations and isotopic ratios will be characterized and potential environmental and health implications discussed.

Mackey, Justin↗

Three-Dimensional Imaging and Quantification of Gas Storativity in Nanoporous Media via X-rays Computed Tomography

This study provides the engineering science underpinnings for improved characterization and quantification of the interplay of gases with kerogen and minerals in shale. Natural nanoporous media such as shale (i.e., mudstone) often present with low permeability and dual porosity, making them difficult to characterize given the complex structural and chemical features across multiple scales. These structures give nanoporous solids a large surface area for gas to sorb. In oil and gas applications, full understanding of these media and their sorption characteristics are critical for evaluating gas reserves, flow, and storage for enhanced recovery and CO 2 sequestration potential. Other applications include CO 2 capture from industrial plants, hydrogen storage on sorbent surfaces, and heterogeneous catalysis in ammonia synthesis. Therefore, high-resolution experimental procedures are demanded to better understand the gas–solid behavior. In this study, CT imaging was applied on the sub-millimeter scale to shale samples (Eagle Ford and Wolfcamp) to improve quantitative agreement between CT-derived and pulse decay (mass balance) derived results. Improved CT imaging formulations are presented that better match mass balance results, highlighting the significance of gas sorption in complex nanoporous media. The proposed CT routine implemented on the Eagle Ford sample demonstrated a 17% error reduction (22% to 5%) when compared to the conventional CT procedure. These observations are consistent in the Wolfcamp sample, emphasizing the reliability of this technique for broader implementation of digital adsorption studies in nanoporous geomaterials.

42 ENGINEERING↗

Multi-Size Proppant Pumping Schedule of Hydraulic Fracturing: Application to a MP-PIC Model of Unconventional Reservoir for Enhanced Gas Production

Slickwater hydraulic fracturing is becoming a prevalent approach to economically recovering shale hydrocarbon. It is very important to understand the proppant’s transport behavior during slickwater hydraulic fracturing treatment for effective creation of a desired propped fracture geometry. The currently available models are either oversimplified or have been performed at limited length scales to avoid high computational requirements. Another limitation is that the currently available hydraulic fracturing simulators are developed using only single-sized proppant particles. Motivated by this, in this work, a computationally efficient, three-dimensional, multiphase particle-in-cell (MP-PIC) model was employed to simulate the multi-size proppant transport in a field-scale geometry using the Eulerian–Lagrangian framework. Instead of tracking each particle, groups of particles (called parcels) are tracked, which allows one to simulate the proppant transport in field-scale geometries at an affordable computational cost. Then, we found from our sensitivity study that pumping schedules significantly affect propped fracture surface area and average fracture conductivity, thereby influencing shale gas production. Motivated by these results, we propose an optimization framework using the MP-PIC model to design the multi-size proppant pumping schedule that maximizes shale gas production from unconventional reservoirs for given fracturing resources.

42 ENGINEERING↗

Organic acid solution having at least one ionic salt and at least one organic acid which is used for rare earth extraction

One or more embodiments relates to a process for extracting Rare Earth Elements (REEs) from REE-bearing underclays, claystones, shales, coal-mining waste, and waste coal. In at least one embodiment the process includes contacting the REE-bearing underclays, claystones, shales, coal-mining waste, and waste coal with an Organic Acid Solution (OAS) comprising at least one organic acid and at least one ionic salt at a predetermined ambient temperature and predetermined pH; and separating the REE from the REE-bearing underclays, claystones, shales, coal-mining waste, and waste coal, forming REE+Yttrium (REY) concentrate.

Verba, Circe↗

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

Hydrogen adsorption and transport in clay-rich geomaterials: Implications for large-volume underground hydrogen storage

Depleted oil and gas reservoirs, characterized by impermeable clay-rich caprocks, are promising sites for large-scale underground hydrogen storage (UHS), which is a key strategy to support hydrogen-based energy systems. However, experimental data on hydrogen storage in clay-rich geomaterials remain scarce. In this work, we experimentally investigated hydrogen adsorption and migration in clay-rich geomaterials in the presence of nitrogen and water under controlled temperatures. Experimental observations showed that hydrogen was adsorbed in dry illite. A dual-porosity transport model was developed to interpret hydrogen transport between large-pore and small-pore domains in illite. The large-pore domain is the space between clay particles (i.e., inter-particle space), whereas the small-pore domain is the nanoscale pore space between clay mineral layers (i.e., inter-layer or intra-particle space). In contrast, nitrogen showed no evidence of adsorption in dry illite because it cannot move into the small-pore domain due to the relatively large kinetic diameter, referred to as the molecular sieving effect. Here, we found that 0.7–1.3 nm is the length scale regulating this molecular sieving effect, matching the interlayer spacing in illite, suggesting that nitrogen is a promising cushion gas in UHS, which aims to maintain adequate pressure in the reservoir for economic operations. In wetted illite, hydrogen was not adsorbed into the interlayer space due to the occupation of adsorption sites by interlayer water, which highlights the critical role of the clay hydration state in controlling hydrogen-clay interactions. Additionally, hydrogen adsorption experiments on crushed shale indicated that the shale surface possessed adsorption sites more favorable for hydrogen than for nitrogen. Through these experiments, we provide new insights into hydrogen storage mechanisms in clay-rich geomaterials and offer valuable laboratory data for evaluating the performance of large-scale UHS systems.

Adsorption↗

Surfactants are Ineffective for Reducing Imbibition of Water-Based Fracturing Fluids in Deep Gas Reservoirs

Minimizing loss of injected hydraulic fracturing fluids into shale along fracture-matrix boundaries is desired because imbibed water restricts gas production and wastes valuable water resources. This problem has motivated the addition of surfactants into water-based hydraulic fracturing fluids in order to reduce the capillary driving force for imbibition. Here, we show that reduction in interfacial tension and wettability alteration has negligible ability to reduce imbibition in deep gas reservoirs. The effectiveness of altering capillary forces acting at the wetting front also depends on the injection pressure acting at the fracture-matrix boundary. The pressure at the interface between the fracture and the shale matrix is constrained between the reservoir pore pressure and formation pressure (rock fracture pressure, also known as breakdown pressure of the rock) and increases with depth to magnitudes that greatly exceed that of capillary pressures. Furthermore, the analyses presented here show that even maximum alteration of interfacial properties that result in strongly hydrophobic interactions between the fracturing fluid and reservoir rock is incapable of significantly reducing imbibition in deep reservoirs. Instead of using surfactants, this analysis points to decreases in wellbore shut-in pressures and shut-in times as practical options for reducing imbibition losses of water-based fluids.

03 NATURAL GAS↗

Effect of Nanoconfinement on NMR Relaxation of Heptane in Kerogen from Molecular Simulations and Measurements

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nanoconfinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. Here, we use atomistic MD simulations to investigate the effects of nanoconfinement on the 1H NMR relaxation times T 1 and T 2 of heptane in kerogen. In the case of T 1 , we discover the important role of confinement in reducing T 1 by ~3 orders of magnitude from that of bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of T 2 , we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces T 2 by ~5 orders of magnitude from the value for bulk heptane. We use the simulated T 2 to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nanopores in kerogen, without additional experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Mineralogy in Controlling Fracture Formation

The presence of fractures in caprocks can pose increased risks in subsurface energy systems and processes like CO 2 sequestration by introducing high-permeability leakage paths. Fracture apertures and permeability can be altered through mineral dissolution and precipitation reactions, but the reactive evolution of fractures is not well understood. In fractures, minerals that are otherwise inaccessible to reactive fluids can become exposed, resulting in mineral reactions unpredicted by bulk formation data. This work seeks to understand the relationship between mineralogy and fracture formation to enhance our understanding of reactive fracture evolution and CO 2 leakage potential. Here, the mineral compositions of mechanically induced fracture surfaces in samples of the Mancos and Marcellus shales have been quantified and compared to those of the near-fracture matrices using imaging and bulk X-ray diffraction (XRD) data. In the Mancos shale, the concentrations of clay minerals are enhanced along fracture surfaces with respect to the bulk, and the fracture is most likely to form at kaolinite–kaolinite interfaces. Further, evaluation of the mineralogical spatial variability through cross-correlation analysis of the surrounding matrix in images of samples cut perpendicular to the fracture shows that clay is 16.7 times more likely to be present than carbonate minerals near the fracture surface. The high correlation persists roughly 200 μm into the surrounding matrix for the Mancos sample and implies that the fracture formed within a defined clay-rich lithofacies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗