Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid injection”

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 55 records · Page 3

Fault Slip and Fluid Flow: Seismic Source Analysis to Assess Role of Multiple Slip Patches in Fault Permeability

The relationship between fault reactivation, microearthquakes (MEQs), and permeability evolution during fluid injection plays a critical role in energy harvesting and waste disposal. Recent studies have demonstrated the possibility of predicting fault permeability using cumulative seismic moments of MEQs quantitatively. To understand the underlying physical processes, we conduct fault reactivation experiments using Utah FORGE granitoid and analyze acoustic emission (AE) signals generated during stepwise increases in fluid injection pressure. Frequency analysis of thousands of calibrated AE signals reveals that fault reactivation produces multiple AE source patches with millimeter-scale radii—smaller than the sample fault radius. The cumulative area of the reactivated patches covers the fault multiple times over (∼10x–50x area) for each pressure step. These findings provide mechanistic insight that measured permeability enhancement is not driven by a single large slip event, but by the sequential and interacting activation of multiple slip patches that create a continuous flow pathway.

Nurshal, M. E. M. [Pennsylvania State University, ↗

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

Deep learning forecasts the spatiotemporal evolution of fluid-induced microearthquakes

Microearthquakes generated by subsurface fluid injection record the evolving stress state and permeability of reservoirs. Forecasting their spatiotemporal evolution is therefore critical for applications such as enhanced geothermal systems, carbon dioxide sequestration and other geoengineering applications. Here we propose a transformer neural network model that ingests hydraulic stimulation history and prior microearthquake observations to forecast four key quantities: cumulative microearthquake count, cumulative logarithmic seismic moment, and the 50th- and 95th-percentile extents of the microearthquake cloud. Applied to the EGS Collab Experiment 1 dataset, the model achieves R2 > 0.98 for the 1-s forecast horizon and R2 > 0.88 for the 15-s forecast horizon across all targets, and supplies uncertainty estimates through a learned standard deviation term. These accurate, uncertainty-quantified forecasts enable real-time inference of fracture propagation and permeability evolution, demonstrating the strong potential of deep-learning approaches to improve seismic-risk assessment and guide mitigation strategies in future fluid-injection operations.

Chung, Jaehong↗

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗

Pore-scale CFD simulations of clay mobilization in natural porous media due to fresh water injection

The present work investigates mechanisms of permeability impairment as a result of low-salinity fluid injection into brine-saturated porous media containing dispersible clays. Here, we present a computational fluid dynamics model at the pore-scale to simulate detachment, migration and straining of fine particles in porous media. The model uses an immersed boundary method to simulate the motion of clay fines in a fluid. In addition to the hydrodynamic forces, we model the Derjaguin-Landau-Verwey-Overbeek forces (DLVO) between clay fines and grains. Our simulations show the impact of the injected fluid's salinity and velocity on the concentration of clay fines retained on the grain surface. We see clay particles dislodging from the grain surface in clusters of up to 12 particles. Our simulation results also demonstrate clogging of narrow pore spaces by the detached particles from the grain surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Abstract for Joint Work Statement 5 (AGMT-1005) between National Energy Technology Laboratory (NETL) and ExxonMobil Research and Engineering Company (EMRE) to Umbrella CRADA (AGMT-0922)

NETL and EMRE will focus on fundamental questions aimed at understanding the controlling factors on tight oil production flow paths and how those flow paths are impacted by the hydraulic fracturing fluids and processes. Key questions for the Parties to explore include the following: Are there geochemical reactions between hydrocarbon producing tight reservoir rocks and water from the hydraulic fracturing process? To what extent does this alteration impact oil production? By further characterizing the flow paths near the fracture-matrix interface, what is the relative impact of geochemical process compared to geomechanical changes on production? This research focuses on rock types found in the Permian Basin. The research drives along two main thrusts: 1) Inorganic geochemistry studies involving representative rock types, connate fluids (hydrocarbons and water), and injected fluids (fracturing water) will be conducted in order to bound the kinetics of key mineral reactions; 2) Computed tomography (CT) for characterization of flow paths as a function of geochemical and geomechanical processes associated with hydraulic fracture creation and reservoir production. The objective of this characterization is to ultimately develop new insights into the role geochemical and geomechanical effects play on production. The work will include studying mechanical damage near the rock-fracture interface and how it impacts flow properties.

02 PETROLEUM↗

Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks

Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.

47 OTHER INSTRUMENTATION↗

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗

Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites

Presentation slides on “Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. Accurate reservoir characterization is fundamental for the safe and efficient operation of carbon capture, utilization, and storage (CCUS) projects. This process is not only critical during initial site selection but also plays a key role in continuously updating reservoir geomodels throughout injection operations to effectively manage subsurface storage sites. Our primary objective is to develop a comprehensive understanding of fracture networks within CO2 injection reservoirs, which are critical for predicting the behavior of injected fluids, ensuring storage integrity post-injection, and mitigating risks such as leakage or induced seismicity.

artifical intelligence / machine learning (AI/ML)↗

Constraining maximum event magnitude during injection-triggered seismicity

Abstract Understanding mechanisms controlling fluid injection-triggered seismicity is key in defining strategies to ameliorate it. Recent triggered events (e.g. Pohang, Mw 5.5) have exceeded predictions of average energy release by a factor of >1000x, necessitating robust methodologies to both define critical antecedent conditions and to thereby constrain anticipated event size. We define maximum event magnitudes resulting from triggering as a function of pre-existing critical stresses and fluid injection volume. Fluid injection experiments on prestressed laboratory faults confirm these estimates of triggered moment magnitudes for varied boundary conditions and injection rates. In addition, observed ratios of shear slip to dilation rates on individual faults signal triggering and may serve as a measurable proxy for impending rupture. This new framework provides a robust method of constraining maximum event size for preloaded faults and unifies prior laboratory and field observations that span sixteen decades in injection volume and four decades in length scale.

58 GEOSCIENCES↗

Alternative, Less Water-Intensive, Fracturing Fluids for Enhanced Geothermal Systems

Four fracturing fluids; water, CO2, CO2 with water, and CO2 with an aqueous solution of a CO2-reactive polymer, poly(allylamine) (PAA, 1wt%) were evaluated using a high-temperature true-triaxial fracturing apparatus and ½ foot side granite cubic samples. All three CO2-based fracturing fluids, CO2, CO2 with water, and CO2 with aqueous PAA fractured granite at higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When partially saturating the rock sample with PAA solution followed by fracturing with CO2, the volume expansion caused by CO2-induced cross-linking of PAA leads to a faster pressure increase due to the associated stress generated and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. It is concluded that CO2 (when injected in HDR) and CO2/PAA are the fluids of choice for stimulation of granitic rock samples under the studied geothermal P/T conditions. CO2/PAA, however, offer the following three additional advantages; 1) it requires significantly lower volumes of CO2, 2) fracture permeability is independent of injection flow rate, and 3) the reversible (previously reported) viscosity increase is beneficial to transport proppants if they are ever developed enhanced geothermal systems.

Jian, Guoqing↗

High Fidelity Modeling of WELL58-32 DFITS

Between 2017 and 2019, a series of injection tests were performed at pilot vertical well 58-32 at the Utah FORGE site. These injection activities were conducted in three distinct zones, isolated with packers and bridge plugs during the tests. Zone 1 consists of an open-hole section at the toe of the well, whereas Zones 2 and 3 are three-meter long cased and perforated sections. Several injection cycles were conducted in each zone. A detailed description of the geology of the FORGE site, of well 58-32 and of all injection cycles can be found in (Xing et al., 2020). The injection cycles included Diagnostic Fracture Injection Tests (DFITs). DFITs comprise two primary stages: an injection (pump-in) phase, where controlled fluid injection is employed to create a small hydraulic fracture in the surrounding rock, followed by a shut-in phase, where injection ceases, allowing observation of the formation’s pressure response. Utilizing methodologies such as the G-function analysis, DFITs provide valuable insights, including fracture closure pressure and the Instantaneous Shut-in Pressure (ISIP) and can be employed to estimate the minimum horizontal stress.

42 ENGINEERING↗

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying Direct SP -Converted Waves Constrains Local Induced Earthquake Depths

Abstract Seismicity in southern Kansas and northern Oklahoma in the past decade has been associated with fluid injections. In southcentral Kansas, the Wellington earthquake catalog is primarily composed of local, low-magnitude events. Approximately 22% of recorded earthquakes over a 2.5 yr period exhibit a seismic phase arriving between the direct P phase and direct S phase with particle motion similar to the P wave. This intermediate phase was identified as an S to P conversion (SP phase) occurring in the sedimentary rocks instead of the hypothesized basement to sedimentary section transition. We exploit the SP-converted phases to improve the depth accuracy of shallow earthquakes and to constrain VP/VS. The revised depth calculations further confirm that these local induced earthquakes are occurring in the shallow crystalline basement, below the sedimentary section in which fluids are injected.

Geochemistry & Geophysics↗

Insights into the hydro-mechanical behavior of a decimeter-scale fracture using the mini-SIMFIP probe

ABSTRACT: Understanding hydro-mechanical couplings in fractured rocks is essential for predicting the rock mass response during high-pressure fluid injection, including the stimulation of enhanced geothermal systems. However, fluid-driven fracture dislocations are challenging to measure due to the need for local displacement data at high fluid pressures. In this study, fluid-driven displacement across a decimeter-scale laboratory fracture was investigated using the mini-SIMFIP (step rate injection method for fracture in-situ properties) probe, which is a smaller version of the SIMFIP tool (Guglielmi et al., 2014). The mini-SIMFIP probe is able to resolve the full 3D displacement vector of a fracture. The probe was installed in one of two boreholes across a decimeter-scale saw-cut granite fracture. Two pressure step injection tests were conducted under the same isotropic stress conditions. By varying injection between the boreholes, we estimated the aperture profile across the fracture. The results show a consistent pressure-dependent opening as long as steady-state flow was maintained. At a certain pressure step, steady-state conditions were no longer achievable and the fracture opened rapidly. The pressure-opening relationship diverged between the two tests and indicated a homogenization of the aperture profile and an increase in the non-linearity of the flow regime.

Osten, J↗

Review: Induced Seismicity During Geoenergy Development—A Hydromechanical Perspective

The basic triggering mechanism underlying induced seismicity traces back to the mid-1960s that relied on the process of pore-fluid pressure diffusion. The last decade has experienced a renaissance of induced seismicity research and data proliferation. An unprecedent opportunity is presented to us to synthesize the robust growth in knowledge. The objective of this article is to provide a concise review of the triggering mechanisms of induced earthquakes with a focus on hydro-mechanical processes. Four mechanisms are reviewed: pore-fluid pressure diffusion, poroelastic stress, Coulomb static stress transfer, and aseismic slip. For each, an introduction of the concept is presented, followed by case studies. Diving into these mechanisms sheds light on several outstanding questions. For example, why did some earthquakes occur far from fluid injection or after injection stopped? Our review converges on the following conclusions: (a) Pore-fluid pressure diffusion remains a basic mechanism for initiating inducing seismicity in the near-field. (b) Poroelastic stresses and aseismic slip play an important role in inducing seismicity in regions beyond the influence of pore-fluid pressure diffusion. (c) Coulomb static stress transfer from earlier seismicity is shown to be a viable mechanism for increasing stresses on mainshock faults. (d) Multiple mechanisms have operated concurrently or consecutively at most induced seismicity sites. (e) Carbon dioxide injection is succeeding without inducing earthquakes and much can be learned from its success. Future research opportunities exist in deepening the understanding of physical and chemical processes in the nexus of geoenergy development and fluid motion in the Earth’s crust.

58 GEOSCIENCES↗

Impact of Acid–Base Stimulation Sequence on Mineral Stability for Tight/Impermeable Unconventional Carbonate-Rich Rocks: A Delaware Basin Case Study

We report mineral precipitation due to reactions with injected fluids during unconventional fracture stimulation is a well-recognized problem. The goal of this study is to evaluate secondary mineral precipitation and permeability attenuation under chemical injection scenarios specific to the Delaware basin. Whole cylindrical cores (2.54 cm diameter and 2.54 cm height) and ground shale (150–250 μm) from the carbonate-rich Bone Spring Formation, Delaware Basin TX (Leonardian), were reacted at 80 °C and 85 bar using a hydraulic fracturing fluid (HFF) recipe and an injection sequence typical of the Delaware Basin. The reacted shales and solutions were analyzed using a variety of laboratory- and synchrotron-based techniques to characterize both the chemical and spatial distributions of secondary mineral precipitation and identify changes in permeability and mineralogy. This carbonate-rich shale (>84% calcite) rapidly neutralized the acidic HFF. Synchrotron-based X-ray fluorescence mapping coupled with X-ray absorption spectroscopy (both bulk and micro) showed that most of the iron was in an oxidized form prior to exposure to HFF and that almost all iron(II) became fully oxidized after the reaction. Scanning electron microscopy images of the ground shale samples primarily identified iron(oxyhydr)oxide microcrystals on grain surfaces. A few small isolated iron-rich areas also contained sulfur, suggesting that some pyrite was preserved when isolated within a calcite crystal but that most was oxidized. The rapid neutralization of the acid spearhead in these carbonate-rich samples demonstrates that the acid spearhead is useful for initiating fractures in extremely calcite-rich rocks but does little to enhance rock permeability. This suggests that the impact of the acid spearhead is significantly smaller for carbonate-rich shales compared to clay-rich shales, which has broad implications for acidizing in carbonate-rich shale formations and iron transformations within these shales.

58 GEOSCIENCES↗