Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗
DOE OSTI2022

Fracture behavior of a melt-infiltration-processed SiC/Si composite at 900°C in argon, air, and steam

The fracture behavior of a melt-infiltration-processed SiC/Si composite, used to mimic the matrices of industrial fiber-reinforced ceramic composites, was examined in different atmospheres and temperatures. Specimens tested in four-point bending at 900 °C in oxygen-gettered argon, dry air, or steam-rich atmospheres exhibited higher average fracture strengths than specimens tested at 25°C. Higher mean fracture strength values were obtained for specimens tested in dry air or in a steam-rich atmosphere at 900°C than for specimens tested in high-purity, oxygen-gettered argon at this temperature. The increased fracture strengths obtained in air and in steam-rich atmospheres coincided with increased specimen oxidation and apparent oxide filling and blunting of flaws in these composites. Further, a transition in the location of catastrophic failure, from sites of preexisting damage created by Vickers indentations for tests in argon to other locations for tests in air or steam-rich atmospheres, was also consistent with such apparent oxide filling/blunting of indentation-induced flaws.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effect of loading direction on the fracture behaviors of cortical bone at a dynamic loading rate

Here, the dynamic cracking processes in porcine cortical bone were visualized in real-time using the high-speed synchrotron X-ray phase-contrast imaging (PCI) technique in three osteon orientations: in-plane transverse, out-of-plane transverse and in-plane longitudinal. The dynamic flexural loading applied on the pre-notched bone specimens was introduced by a modified Kolsky compression bar. High-speed X-ray images of the entire loading events were documented with a high-speed camera. Three-dimensional X-ray micro-computed tomography was conducted to examine the intact microstructures and obtain the basic material properties of the bone material used for mechanical characterizations. The onset location, where crack initiated, and the subsequent direction, along which the incipient crack propagated, were measured quantitatively using the high-speed X-ray images and the latter was found dependent on the osteon direction significantly. The crack propagation velocities were dependent on crack extension over the entire crack path significantly for all the three directions while the initial velocity for in-plane longitudinal direction was lower than the other two directions. Straight-through crack paths were observed for in-plane longitudinal specimens while the cracks were deflected and twisted in the in-plane transverse direction. For out-of-plane transverse direction, the cracks follow paths with tortuosity fall in between the other two directions, showing a mixed mode of fractures of the former two extreme cases. The toughening mechanisms, visualized by the high-speed X-ray images, and the corresponding fracture toughness, evaluated in terms of fracture initiation toughness and crack growth resistance curve (R-curve), were also found significantly different among the three osteon directions, suggesting an overall transition from brittle to ductile-like fracture behaviors at the dynamic displacement rate (5.4 m/s) as the osteon orientation varies from in-plane longitudinal to out-of-plane transverse, and to in-plane transverse eventually.

42 ENGINEERING↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2020

Visco-inertial gas flow through wellbore cement fractures

Understanding the nature of fluid flow through fractured wellbore cement is fundamental for evaluating the leakage potential and risk assessments of leaky wellbores. In this study, the conditions that require considering visco-inertial flow for describing the gas flow through wellbore cement fractures were investigated. Nitrogen gas flow tests were conducted on fractured cement samples under varying pressure conditions and flow rates, covering both viscous and visco-inertial flow regimes. The data substantially deviated from Darcy's law at higher flowrates and were well-fit to Forchheimer's equation for visco-inertial flow. The inertial coefficient and critical Reynolds number were expressed as a function of the hydraulic aperture. Finally, the empirical function obtained from the experiments was used as an input to numerical simulations which showed the significant role of visco-inertial flow in wellhead pressure build-up and leakage rates, and demonstrated the importance of visco-inertial flow when modeling gas flow through wellbore cement fractures.

03 NATURAL GAS↗
DOE OSTI2020

Numerical study of the effect of propped surface area and fracture conductivity on shale gas production: Application for multi-size proppant pumping schedule design

Hydraulic fracturing is a technique extensively used in the oil and gas industry, where water, proppant (sand) and additives are injected into unconventional reservoirs to enhance the recovery of shale hydrocarbon. Although some previous studies have developed pumping schedules that maximize gas production for a single-size proppant, there are very few studies that consider the effect of varying proppant diameters across pumping stages on shale gas production. Motivated by this, we carried out an extensive sensitivity analysis to determine the effect of different proppant diameters on the average fracture conductivity (FC), average propped surface area (PSA) and cumulative shale gas production volume. We found out that the cumulative shale gas production volume depends on both the average PSA and average FC. We also found out that small-diameter proppant resulted in higher average PSA and lower average FC, whereas large-diameter proppant resulted in lower average PSA and higher average FC. Hence, we designed a multi-size proppant pumping schedule considering both of these parameters into account for simultaneously propagating multiple fractures to maximize shale gas production from unconventional reservoirs. Since the size of injected proppant particles determines the average PSA and average FC for the propped hydraulic fractures, we developed a novel framework called Sequentially Interlinked Modeling Structure (SIMS) to predict the average PSA, average FC and cumulative shale gas production volume at the end of 10 years for a given pumping schedule. Then, we used this SIMS framework to obtain a multi-size proppant pumping schedule that maximizes shale gas production. In closing, we demonstrated that the obtained pumping schedule gives a cumulative shale gas production volume greater than the values obtained from the existing pumping schedules.

03 NATURAL GAS↗
DOE OSTI2024

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Influence of fatigue precracking and specimen size on Master Curve fracture toughness measurements of EUROFER97 and F82H steels

EUROFER97 and F82H steels are two leading reduced-activation ferritic-martensitic (RAFM) steels for fusion first wall and blanket applications. Exposure to the harsh environment of fusion reactors can result in severe degradation of fracture toughness. Thus, the post-irradiation evaluation of fracture toughness is critical for understanding the material behavior. Due to the space constraints of irradiation facilities and challenges in controlling a uniform irradiation condition for large size specimens, the development of small specimen test techniques (SSTT) is indispensable to evaluate the performance of irradiated materials. In this study, we evaluated specimen size effects on the Master Curve fracture toughness of EUROFER97 and F82H steels. A wide variety of specimens, including 0.5 T compact tension (C(T)) specimens, 0.16 T mini-compact tension (miniC(T)) specimens, and 1.65 mm miniature bend bar specimens, were tested. The testing methodology was based on the Master Curve method in the ASTM E1921 standard. No specimen size effect was observed in 0.5 T C(T) and 0.16 T miniC(T) specimens on the Master Curve reference temperature T0, while 1.65 mm miniature bend bar specimens yielded a higher T 0Q . A strong effect of fatigue precracking on T 0 for 0.5 T C(T) and 0.16 T miniC(T) specimens was observed, such that testing on specimens with skewed fatigue precrack fronts resulted in lower T 0 than for specimens with ASTM standard qualified straight fatigue precrack fronts. The results highlight the importance of experimental quality control in developing SSTT for Master Curve fracture toughness testing. Lastly, we also evaluated and provided recommendations on the minimum number of specimens needed for each specimen type for yielding reliable T 0Q values.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of irradiation embrittlement effect on fatigue life of a pressurized-water reactor pressure vessel using the fracture toughness master curve approach

The reactor pressure vessel (RPV) is a critical structural component in pressurized-water reactors, and it is designed to withstand extreme conditions, such as high pressures, elevated temperatures, and prolonged radiation exposure. Ensuring RPV integrity is essential for the safe and reliable long-term operation of nuclear power plants, especially as aging mechanisms such as fatigue and irradiation embrittlement pose increased risks. Fatigue, caused by cyclic thermal and mechanical loading, can lead to crack initiation in localized high-stress regions. Simultaneously, neutron irradiation, particularly in the beltline region, progressively reduces fracture toughness, increasing susceptibility to brittle fracture. These combined effects of fatigue and irradiation embrittlement potentially impact the RPV structural integrity, necessitating fitness-for-service assessments. This study applies the fracture toughness master curve approach to evaluate the impact of irradiation embrittlement on RPV fatigue life. A coupled thermo-mechanical stress analysis identifies critical stress locations under normal service transients, pinpointing regions most vulnerable to fatigue crack initiation and growth. Stress intensity factors for postulated flaws at these locations are calculated, enabling an assessment of fatigue life under irradiated and unirradiated conditions. The results indicate that neutron irradiation embrittlement accelerates the conditions in which a critical crack can form and lead to failure, particularly at lower temperatures. The failure occurs where reduced fracture toughness limits the material’s resistance to crack growth. Axial cracks at nozzle corners are the most life limiting without irradiation, while circumferential cracks demonstrate longer fatigue lives. The findings highlight the importance of incorporating irradiation effects into fatigue life predictions to ensure the long-term structural integrity of RPVs.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗
DOE OSTI2021

High temperature fracture toughness of single-layer CrAlN and CrAlSiN hard coatings

Hardness and toughness are among of the most important physical properties of hard coatings and it is important to analyze them at temperatures present during coating applications. CrAlN and CrAlSiN are hard coating materials frequently used in high temperature applications. Here, while both possess high oxidation resistance and thermal stability, they differ in their microstructure. In this work, the fracture toughness of magnetron-sputtered single-layer CrAlN and CrAlSiN hard coatings was determined up to 700 °C by in-situ scanning electron microscopy micro-cantilever deflection tests. The results show a decreasing trend of the fracture toughness with increasing temperature. While both coatings' fracture toughness trends are similar up to 400 °C, they diverge at higher temperatures. In general, the evolution of the fracture toughness with increasing temperature shows similarities to the observed changes in hardness in the same temperature range.

36 MATERIALS SCIENCE↗
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 OSTI2021

Employing Fracture Statistics to Track Cell Reliability Through Module Fabrication

In this work we demonstrate how fracture statistics may be employed to track cell reliability through module fabrication steps. We were able to detect how cell cutting equipment could reduce cell fracture strain by 25 % as well as increase its variability. A smaller reduction was found to result from the ribbon soldering operation, and all cells measured were weaker when strained in an orientation perpendicular to the gridlines. Cell fracture measurements were found to be consistent with the number of cracked cells experienced in full sized modules exposed to mechanical loading. A structural mechanics model was created to elucidate how the requisite fracture strain is developed in a cell encapsulated within a PV module and demonstrates that this strain decreases more than 50 % from its center to edge when placed within the inner loading span of four-point flexure. The associated implications on the apparent cell strength is examined with an effective volume treatment.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2026

Seismic Elastic Double-Beam Characterization of Faults and Fractures for CO₂ Storage Site Selection

Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.

02 PETROLEUM↗
DOE OSTI2021

Effects of Carbonate Minerals on Shale-Hydraulic Fracturing Fluid Interactions in the Marcellus Shale

Natural gas extracted from tight shale formations, such as the Marcellus Shale, represents a significant and developing front in energy exploration. By fracturing these formations using pressurized fracturing fluid, previously unobtainable hydrocarbon reserves may be tapped. While pursuing this resource, hydraulic fracturing operations leave chemically complex fluids in the shale formation for at least two weeks. This provides a substantial opportunity for the hydraulic fracturing fluid (HFF) to react with the shale formation at reservoir temperature and pressure. In this study, we investigated the effects of the carbonates on shale-HFF reactions with a focus on the Marcellus Shale. We performed autoclave experiments at high temperature and pressure reservoir conditions using a carbonate-rich and a decarbonated or carbonate-free version of the same shale sample. We observed that carbonate minerals buffer the pH of the solution, which in turn prevents clay dissolution. Carbonate and bicarbonate ions also scavenge reactive oxidizing species (ROS), which prevents oxidation of shale organic matter and volatile organic compounds (VOCs). Carbonate-free samples also show higher pyrite dissolution compared to the carbonate-rich sample due to chelation reactions. This study demonstrates how carbonate minerals (keeping all other variables constant) affect shale-HFF reactions that can potentially impact porosity, microfracture integrity, and the release of heavy metals and volatile organic contaminants in the produced water.

58 GEOSCIENCES↗