Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale Formations”

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

Sensitivity and history match analysis of a carbon dioxide “huff-and-puff” injection test in a horizontal shale gas well in Tennessee

Due to improvements in well development, shale gas production has gained importance recently, especially in the United States. To improve gas production and develop a better understanding of shale reservoirs, researchers are conducting field tests to monitor CO 2 storage and enhanced gas recovery. Reservoir simulations can then utilize these field results for sensitivity studies, uncertainty analysis, history matching, gas production forecasts, and CO 2 storage capacity estimations. One such field test was performed for the Chattanooga Shale formation in Morgan County, Tennessee. Approximately 463 tonnes (510 tons) of CO 2 was successfully injected into a hydraulically fractured horizontal well over a thirteen-day period in March 2014. The injection test was achieved in four steps: pre-injection, injection, soaking, and flowback. In this paper, those steps were modeled with a reservoir simulator to match the historic production and injection rates with gas component compositions and forecast the production for five years. The reservoir fluid was modeled as a multi-component gas, including CH 4 , C 2 H 6 , C 3 H 8 , N 2 , and CO 2 . Langmuir constants, reservoir pressure, and fracture network volume were adjusted to match simulation results with observed production rates and gas composition. Results showed that, under the same reservoir conditions, each gas component behaves differently by way of compositional changes in production. CH4 behaves like N2, while CO 2 behaves like heavier hydrocarbons such as C 2 H 6 and C 3 H 8 . CO 2 plume results showed that, after injection, the produced CO 2 is mostly derived from fractured limestone (i.e., Fort Payne Formation) because the injected CO 2 is not adsorbed in limestone but rather by shale formations.

03 NATURAL GAS↗
DOE OSTI2023

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↗
DOE OSTI2020

Characterizing Pore-Scale Geochemical Alterations in Eagle Ford and Barnett Shale from Exposure to Hydraulic Fracturing Fluid and CO 2 /H 2 O

As demand increases for an affordable energy source that is tied to an environmental obligation to reduce greenhouse gas emissions and water usage, there is a growing consideration in shale production utilizing processes such as 1) enhancing hydrocarbon recovery via carbon dioxide (CO 2 ) flooding, 2) using CO 2 as a fracturing agent to minimize water use, and 3) storing CO 2 in depleted shale formations to mitigate emissions to the atmosphere. Understanding the geochemical reactions and alterations that occur as shale is exposed to fluids and CO 2 is necessary to develop and optimize each of these processes for field applications. While the majority of shale formations are stimulated using traditional fracturing fluid, some may be fractured using CO 2 or other non-traditional means. We examine the effect fracturing fluid has on shale and how it behaves with secondary exposure to dry CO 2 or CO 2 -saturated water using in situ Fourier Transform infrared spectroscopy (FTIR), feature relocation scanning electron microscopy (SEM), and surface area and pore size analysis using volumetric gas sorption. These techniques were performed on Eagle Ford and Barnett shale samples that were exposed to fracturing fluid and unexposed (as received). Shales that have been exposed to traditional fracturing fluid experienced two reaction fronts. The first reaction front was formed during exposure to the fracturing fluid (pH of ~1.4). A secondary reaction front was formed as a result of CO 2 -saturated fluid exposure in the form of carbonic acid (pH ~5.6). These two different reaction mechanisms drove multiple dissolution and precipitation cycles which altered petrophysical properties of the shale and could lead to a significant impact on flow pathways. FTIR showed that equilibration of carbonate dissolution and precipitation cycles could take as long as 35 days. Samples exposed to fracturing fluid showed significantly less carbonate reactivity compared to those exposed to water. Pore size analysis results indicate exposure to fracturing fluid blocked small nanopores (10-15 nm) reducing BET surface area and total pore volume. SEM results show barite precipitated heavily during exposure to fracturing fluid. It appeared that carbonic acid was able to extract sulfur from organic matter to form gypsum evaporites. The mineralogical (barite precipitation and calcite dissolution/precipitation) and pore-scale alterations observed in these samples may lead to enhancement of flow pathways for injected CO 2 or produced hydrocarbons.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Accessibility of Pores to Methane in New Albany Shale Samples of Varying Maturity Determined Using SANS and USANS

The accessibility of pores to methane has been investigated in Devonian New Albany Shale Formation early-mature (Ro = 0.50%) to post-mature (Ro = 1.40%) samples. A Marcellus Shale Formation sample was included to expand the maturation range to Ro 2.50%. These are organic matter-rich rocks with total organic carbon (TOC) values of 3.4 to 14.4% and porosity values of 2.19 to 6.88%. Contrast matching small-angle neutron scattering (SANS) and ultra-small angle neutron scattering (USANS) techniques were used to generate porosity-related data before and after pressure cycling under hydrostatic (in a vacuum and at 500 bar of deuterated methane) and uniaxial stress (0 to ca. 350 bar) conditions. Our results showed that the accessible porosity was small for the samples studied, ranging from zero to 2.9%. No correlation between the accessible porosity and TOC or mineralogical composition was revealed, and the most likely explanation for porosity variation was related to the thermal transformation of organic matter and hydrocarbon generation. Pressure caused improvements in accessible porosity for most samples, except the oil window sample (Ro = 0.84%). Our data show that densification of methane occurs in nanopores, generally starting at diameters smaller than 20 nm, and that the distribution of methane density is affected by pressure cycling.

03 NATURAL GAS↗
DOE OSTI2023

Induced microseismicity and tremor signatures illuminate different slip behaviours in a natural shale fault reactivated by a fluid pressure stimulation (Mont Terri)

SUMMARY Fault slip induced by fluid perturbation in shale formations may only lead to a sparse seismicity. However, fault slip may strongly impact the integrity of shale formations that serve as caprocks for geological reservoirs holding buoyant fluids such as CO2, natural gas or hydrogen. A better understanding of the fluid reactivation processes of fault and the seismic triggering process is therefore critical for reservoir monitoring and fault stability. Here we analyse the seismic responses of a shale fault exposed to fluid pressurization during an in situ field-scale injection experiment at ∼340 m depth in the Mont Terri underground research laboratory (Switzerland). Two main types of seismic signals are observed as the fault was activated and started to slowly slip. After an aseismic phase, we observed tremor signatures and an increase in noise amplitude, which were directly associated with the slowly propagating fault slip in response to fluid injection. These signatures were later followed by micro-earthquakes that seem to occur further away from the fluid-pressurized area. We interpret these micro-earthquakes to be triggered by stress perturbations from the main slip growth. These two classes of seismic responses therefore highlight two different processes. Tremors seem to be a more direct observation for the fluid-induced slip propagation than micro-earthquakes. Even hidden in the noise, they precede earthquake failures, thus providing a useful tool for monitoring fluid leakage activated by slow deformation on low permeable shale faults, with applications for sealing integrity of caprocks.

De Barros, Louis (ORCID:0000000255419162)↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2023

Small-Angle Neutron Scattering Investigation of Oil Recovery in Mineralogically Distinct Wolfcamp Shale Strata

Understanding and improving hydrocarbon yields during enhanced oil recovery (EOR) in unconventional reservoirs is complicated by the intrinsic mineralogical and geochemical heterogeneity of shale formations. Here, in this study, we utilized small-angle neutron scattering (SANS) and ultra-small-angle neutron scattering (USANS) to investigate the degree of oil retention and its location in the nanoporous shale matrix for two mineralogically distinct shale samples. The two samples, dubbed “dark” and “light” based on their color, were taken from adjacent strata in a Wolfcamp shale core. While both samples contained kerogen, the dark sample contained more kerogen and clay (43.7 wt %) while the light sample contained more calcite (54.9 wt %). Samples were presaturated with decane, a model hydrocarbon, prior to pressure cycling with methane. Results showed significantly more retention of decane in 1.5–10 nm radius pores of both, likely indicating that oil is retained within kerogen nanopores. Although the dark sample had a higher porosity of 8.7%, versus 3.3% for the light sample, more pores were accessible to decane and a higher percentage of the imbibed decane was removable from the light sample compared to the dark sample. The majority of decane was not recoverable for the dark sample, indicating that EOR with methane can be challenging. These new findings can help to model expected recoveries of in-place oil from heterogeneous shale formations, as well as inform improved EOR strategies.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Comparison of Geochemical Reactivity of Marcellus and Caney Shale Based on Effluent Analysis

ABSTRACT: In this comparative study, we analyzed the changes in elemental concentrations of hydraulic fracturing fluids after interaction with the Marcellus and Caney Shale formations. The focus was on assessing the inherent risks and environmental implications associated with flowback waters, including their impact on soil, ground, and drinking water quality, and human health safety. The chemical compositions of effluents were determined through Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES) and Mass spectrometry (ICP-MS). The Marcellus Shale, showed concentrations of Cd, averaging 0.380 ppm far exceeding safe water thresholds. Significant levels of As, Se, B, and Pb were detected in both shales, raising concerns about soil and water contamination. Analyzing the effluents from representative samples of sections of the Marcellus (S2 and S7) and Caney (R1 and R2) indicates different geochemical responses over 4 weeks of experiments. This comparison underscores the chemical changes and environmental considerations linked to hydraulic fracturing across shale formations, suggesting the value of tailored monitoring and regulatory measures for each type. 1. INTRODUCTION The Caney and Marcellus shales, differing in geology and geochemistry, represent distinct unconventional reservoirs. The Caney shale, is more ductile, with higher produced water volumes, (Smith et al., 2022) contrasts with the brittle Marcellus shale known for lower brine production but significant data availability. This study aims to elucidate the possible environmental health and safety issues that may arise from the hydraulic fracturing processes. The interaction between fracturing fluids and clays presents a significant challenge. The primary base of these fluids is water, which, when introduced to clay, can induce swelling and constrict flow pathways. This phenomenon is attributed to water molecules infiltrating the layers of clay, particularly in 2:1-type clays, leading to an increased distance between layers. To counteract this, clay stabilizers are employed (Awejori et al., 2021). Despite their effectiveness, these stabilizers are considered temporary solutions. Upon completion of the fracturing process, a concomitant amount of contaminated waters (flow back), with varied content is collected at the surface. From these, we can infer the geochemical reactions and the environmental challenges associated with these waters. Flowback waters can be reinjected or used for other purposes such as irrigation. This requires adequate screening and treatment for safe use.

Dje, L. B.↗
DOE OSTI2023

Organic matter sulfurization and organic carbon burial in the Mesoproterozoic

Throughout the Proterozoic Era, sedimentary organic carbon burial helped set the pace of global oxygenation and acted as a major modulator of atmospheric CO 2 and climate. Although Proterozoic rocks generally contain low concentrations of organic matter (OM), there are key exceptions to this rule, including the relatively OM-rich Arctic Bay shales from Baffin Island, Canada (Bylot Supergroup, Borden Basin, ~1.05 Ga). The mechanisms driving elevated OM concentrations in these and other Proterozoic shales remain poorly understood. In the Mesozoic and Cenozoic, organic matter sulfurization can be a major driver of enhanced OM burial across a range of redox conditions comparable to those inferred for many Proterozoic environments. Therefore, in this study, we evaluate the role of sulfurization in driving OM preservation in the Mesoproterozoic Borden Basin and discuss its relevance to Proterozoic systems in general. We present the first evidence for syngenetic-to-early-diagenetic OM sulfurization in a Proterozoic basin, which begins to fill a several-billion-year gap in our record of organic S across Earth history. We find that OM sulfurization was particularly extensive in shales from a relatively shallow-water section (Alpha River) but less extensive in shales deposited in deeper water (Shale Valley), which is consistent with models that infer sulfidic ‘wedges’ or O 2 -minimum-zone-type structures on shelf margins at least intermittently at this time. At the shallower site, organic S and pyrite are similarly 34 S-depleted and thus likely formed at roughly the same time near the sediment–water interface under conditions previously interpreted to have been ferruginous to intermittently sulfidic. In contrast, at the deeper-water site, large S-isotope differences between pyrite and organic S along with low apparent OM sulfurization intensities indicate that pyrite formation was favored over OM sulfurization during early sedimentation under variable but primarily ferruginous conditions. Although Mesoproterozoic biomass can be substantially sulfurized, indicators of sulfurization intensity are not correlated with OM concentrations, and therefore sulfurization does not appear to have been the primary driver of enhanced OM concentrations in Arctic Bay Formation shales. The link between sulfurization and total OM preservation may have been modulated during the deposition of Arctic Bay Formation shales by differences in iron availability, nutrient cycling, and particle dynamics in the Mesoproterozoic.

58 GEOSCIENCES↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Developing Methods to Assess Changes in Mechanical Properties of Shale Modified by Engineered Mineral Precipitation

Fractures in subsurface shale formations serve multiple purposes, for example, in the recovery of resources in hydraulic fracturing or as potential harmful leakage passages through caprocks that may contribute undesired fluids to the atmosphere or functional groundwater aquifers. A proposed method to seal or influence fracture properties is Ureolysis-Induced Calcium Carbonate Precipitation (UICP), a bio-mineralization technology driven by the enzymatic hydrolysis of urea, resulting in the formation of calcium carbonate. Sporosarcina pasteurii is a common microbe used as the source of the urease enzyme that catalyzes the chemical reaction. The resulting calcium carbonate can bridge the gaps in fractured shale and reduce fluid flow through fractures. However, there is little information on how this process affects the mechanical properties of the resulting biomineralized shale. This study represents the first step toward determining the influence of UICP treatment on shale material and its subsequent mechanical strength properties. This methods development study aims to determine the effect that temperature has on the tensile strength of intact, unfractured shale cores (2.54 cm (1 in) diameter, 5.08 cm (2 in) long). Tensile strength was determined indirectly using a modified Brazilian test where the splitting tensile strength is attained by applying a compressive load onto the core. Shale cores from Eagle Ford and Wolfcamp formations were tested at both room temperature and 60°C to determine if increased temperature influences the tensile strength of the rock. This data will help to assess the necessity of testing biomineralized cores at temperature. Though 60°C may not mimic subsurface temperatures of the shales used in this study, it was chosen due to limitations of the UICP process while still approaching temperatures of shale formations. This project aims to evaluate what effect temperature has on the mechanical properties of intact shale cores so that engineered or natural rock fractures that are sealed by biomineralization can be better understood.

Bedey, Kayla↗
DOE OSTI2020

Reactive Transport Modeling of Shale–Fluid Interactions after Imbibition of Fracturing Fluids

Injection of hydraulic fracturing fluid (HFF) into shale formations for unconventional oil/gas production results in chemical reactions in the shale matrix. Our recent experimental study determined the depths to which different types of reactions between the shale matrix and the HFF extended. In the present study, we built continuum-scale reactive transport models to understand the coupling of chemical reactions and the transport of aqueous species in these shale–HFF systems. Calibration of the model with our previous experimental results reveals that it takes hours to months for the shale matrix to completely neutralize the imbibed acids, depending primarily upon the carbonate content of the shale. Both the HFF pH and pore pH affect the location of barite precipitation, resulting in unique barite precipitation profiles extending millimeters into calcite-rich Eagle Ford shale but only tens of micrometers into low-carbonate Marcellus shale. In addition, dissolved oxygen and extracted bitumen are key to reproducing the experimental observation of Fe(III) (oxyhydr)oxide formation in the shale matrix as a result of pyrite dissolution in the shales. A comparison between the modeling results of porosity in the present study to experimentally measured permeabilities in our previous study suggests that chemical reactions occurring at a greater depth than the observable reaction zone might have impacted permeability. Our model serves as a benchmark for efficiently modeling water–rock interactions in similar systems where bulk rock samples react with a solution in batch reactors. Important reactive transport processes were ascertained via modeling, which allows for quantitative prediction of shale–HFF interactions in shale matrices given the shale and HFF compositions.

58 GEOSCIENCES↗
DOE OSTI2021

Toxicity of hydraulic fracturing wastewater from black shale natural-gas wells influenced by well maturity and chemical additives

Hydraulic fracturing of deep shale formations generates large volumes of wastewater that must be managed through treatment, reuse, or disposal. Produced wastewater liberates formation-derived radionuclides and contains previously uncharacterized organohalides thought to be generated within the shale well, both posing unknown toxicity to human and ecological health. Here, we assess the toxicity of 42 input media and produced fluid samples collected from four wells in the Utica formation and Marcellus Shale using two distinct endpoint screening assays. Broad spectrum acute toxicity was assessed using a bioluminescence inhibition assay employing the halotolerant bacterium Aliivibrio fischeri, while predictive mammalian cytotoxicity was evaluated using a N-acetylcysteine (NAC) thiol reactivity assay. The acute toxicity and thiol reactivity of early-stage flowback was higher than later produced fluids, with levels diminishing through time as the natural gas wells matured. Acute toxicity of early stage flowback and drilling muds were on par with the positive control, 3,5-dichlorophenol (6.8 mg L -1 ). Differences in both acute toxicity and thiol reactivity between paired natural gas well samples were associated with specific chemical additives. Samples from wells containing a larger diversity and concentration of organic additives resulted in higher acute toxicity, while samples from a well applying a higher composition of ammonium persulfate, a strong oxidizer, showed greater thiol reactivity, predictive of higher mammalian toxicity. Both acute toxicity and thiol reactivity are consistently detected in produced waters, in some cases present up to nine months after hydraulic fracturing. These results support that specific chemical additives, the reactions generated by the additives, or the constituents liberated from the formation by the additives contribute to the toxicity of hydraulic fracturing produced waters and reinforces the need for careful consideration of early produced fluid management.

03 NATURAL GAS↗
DOE OSTI2022

Tuscaloosa Marine Shale: Seal or Source? Petrophysical Comparative Study of Wells in SE Louisiana and SW Mississippi

The Tuscaloosa Marine Shale (TMS) is a versatile Late Cretaceous shale formation present in central and SE Louisiana and SW Mississippi, which drew attention because of the various roles played within the Tuscaloosa Group. In this paper, it is debated whether the Tuscaloosa Marine Shale can act as a source, reservoir, or seal all throughout the shale play or only in certain areas. Well log and core data from Adams County, Mississippi, are compared to data from East Feliciana Parish in Louisiana. Conclusions were drawn based on the results of well log analysis, X-ray Diffraction (XRD), porosity–permeability measurements, programmed pyrolysis, and fracture analysis. It was shown that the Tuscaloosa Marine Shale interval in SE Louisiana consists of important amounts of calcite, exhibits multiple natural fractures, has porosity values as high as 9.3%, and shows a TOC content of up to 2.8 wt%. On the other hand, samples from a well at the Cranfield field, MS, are characterized by considerably lower TOC values of around 0.88 wt%, porosities between 0.33% and 4%, and no serious fracturing. The formation demonstrates better reservoir and source potential in SE Louisiana and reliable CO2 sealing capacity in SW Mississippi. The analysis presented in this paper represents a holistic approach to the characterization of shale formations, is applicable to other plays around the world, and can be used as an integral part of CO2 sequestration or hydraulic fracturing programs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Aperture size distribution, length, and preferential location of bed-parallel veins in shale

Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.

58 GEOSCIENCES↗
DOE OSTI2020

Assessing the thermal maturity of black shales using vitrinite reflectance: Insights from Devonian black shales in the eastern United States

Thermal maturity of source rocks indicated by vitrinite reflectance (R o ) is an important and reliable parameter to determine the petroleum potential of sedimentary basins. Originally, R o was used to characterize the degree of coalification of coals and thus works well for coal seams, but it has limitations when used to determine the thermal maturity of black shale successions. A suite of six coal-shale pairs from the Upper Devonian black shale formations in the eastern United States (New Albany Shale, Ohio Shale, Dunkirk Shale, and Rhinestreet Shale) was selected to study the applicability of R o in assessing the thermal maturity of black shale successions. The results show that vitrinite in the studied coal samples is dominated by collotelinite, whereas vitrinite in black shales occurs as small dispersed particles (~5 μm) in the mineral matrix. When comparing the size and morphology of dispersed vitrinite particles in shales and collotelinite fragments in coals, vitrinite in shales and adjacent coals should have the same origin. The measured mean random R o of vitrinite in coals ranges from 0.51–0.68%, and is 0.04–0.11% (average 0.07%) lower than that of dispersed vitrinite particles in enclosing shales. This observation contrasts with previously reported R o suppression in liptinite-rich black shales. A high standard deviation of R o measurements indicates a highly heterogeneous nature of vitrinite in black shales. Three mechanisms can be envisioned to contribute to the higher R o of dispersed vitrinite particles in black shales relative to enclosed coals. First, small vitrinite particles may become more oxidized because small particles are more frequently suspended during transport to the site of deposition and reworking of bottom sediments after deposition than large pieces of driftwood that will turn into enclosed coal lenses during burial diagenesis. Second, misidentification of zooclast (e.g., chitinozoan) fragments as vitrinite due to loss of diagnostic morphology. Third, various macerals in the vitrinite group (e.g., corpogelinite and collotelinite) may have had different original reflectance values. Based on the petrographic characteristics of vitrinite and vitrinite-like particles, caution should be applied when using the R o values of dispersed vitrinite and vitrinite-like particles in black shales as indicators of thermal maturity.

58 GEOSCIENCES↗
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 OSTI2022

Influence of sequential stimulation practices on geochemical alteration of shale

Water-based hydraulic fracturing fluids (HFFs) can chemically interact with formation shale, resulting in altered porosity and permeability of the host rock. Experimental investigations of spatial and temporal shale-HFF interactions are helpful in interpreting chemical compositions of the injectate, as well as predicting alteration of hydraulic properties in the reservoir due to mineral dissolution and precipitation. Most bench-top experiments designed to study shale-HFF chemical interactions, either using batch reactors or flow-through setups, are carried out assuming that the acid spearhead has already become mixed with neutral HFFs. During operations, however, HFFs are typically injected according to a sequenced pumping schedule, starting with a concentrated acid spearhead, followed by multiple additions of near-neutral pH HFFs containing chemical amendments and proppant. In this study, we use geochemical modeling to consider whether this pre-mixed experimental protocol provides results directly comparable to a sequential discrete fluid-shale interaction protocol. Our results show that for the batch system, the transient evolution in major ion concentrations is faster with the sequential procedure. After 2 h of reaction time, the two protocols converge to the same aqueous concentrations. In a flow-through geometry, the pre-mixed model predicts extensive chemical alteration close to the injection point but negligible alteration downstream. In contrast, the sequential model predicts mineral reactions over hundreds of meters along the flow path. The extent of shale alteration in the sequential model at a given location depends on shale mineralogy and where the acid spearhead resides during the shut-in period. The predictive model developed in this study can help experimentalists to design bench-top tests and operators to better translate the results of laboratory experiments into practical applications.

Li, Qingyun↗