Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15

In-Situ Stress Measurements at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) Site

A scientific injection campaign was conducted at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site in 2017 and 2019. The testing included pump-in/shut-in, pump-in/flowback, and step rate tests. Various methods have been employed to interpret the in-situ stress from the test dataset. This study focuses on methods to interpret the minimum in-situ stress from step rate, pump-in/extended shut-in tests data obtained during the stimulation of two zones in Well 58-32. This well was drilled in low-permeability granitoid. A temperature of 199 °C was recorded at the well’s total depth of 2297 m relative to the rotary Kelly bushing (RKB). The lower zone (Zone 1) consisted of 46 m of the openhole at the toe of the well. Fractures in the upper zone (Zone 2) were stimulated between 2123–2126 m measured depths (MD) behind the casing. The closure stress gradient variation depended on the depth and the injection chronology. The closure stress was found to increase with the pumping rate/volume. This stress variation could indicate that poroelastic effects (“back stress”) and the presence of adjacent natural fractures may play an important role in the interpretation of fracture closure stress. Further, progressively increasing local total stresses may, consequently, have practical applications when moderate volumes of fluid are injected in a naturally fractured or high-temperature reservoir. The alternative techniques that use pump-in/flowback tests and temperature signatures provide a valuable perspective view of the in-situ stress measurements.

15 GEOTHERMAL ENERGY↗

A Mathematical Model for the Multiphase Transport and Reaction Kinetics in a Ladle with Bottom Powder Injection

A computation fluid dynamics-population balance model-simultaneous reaction model (CFD-PBM-SRM) coupled model has been proposed to study the multiphase flow behavior and refining reaction kinetics in a ladle with bottom powder injection, and some new and important phenomena and mechanisms are presented. For the multiphase flow behavior, the effects of bubbly plume flow, powder particle motion, particle-particle collision and growth, particle-bubble collision and adhesion, and powder particle removal into top slag are considered. For the reaction kinetics, the mechanisms of multicomponent simultaneous reactions, including Al, S, Si, Mn, Fe, and O, at the multi-interface, including top slag-liquid steel interface, air-liquid steel interface, powder droplet-liquid steel interface, and bubble-liquid steel interface, are presented, and the effect of sulfur solubility in the powder droplet on the desulfurization is also taken into account. Model validation is carried out using hot tests in a 2-t induction furnace with bottom powder injection. The result shows that the powder particles gradually disperse in the entire furnace; in the vicinity of the bottom slot plugs, the desulfurization product CaS is liquid phase, while in the upper region of the furnace, the desulfurization product CaS is solid phase. The predicted sulfur contents by the present model agree well with the measured data in the 2-t furnace with bottom powder injection.

36 MATERIALS SCIENCE↗

Geochemical Assessment for Carbon Sequestration in the Conasauga Group, Northwest Georgia, USA

Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.

42 ENGINEERING↗

Experiments and Modeling of Proppant Embedment and Fracture Conductivity for the Caney Shale, Oklahoma, USA

ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.

Katende, A.↗

Validation of single-fluid and two-fluid magnetohydrodynamic models of the helicity injected torus spheromak experiment with the NIMROD code

We present a comparison study of 3-D pressureless resistive MHD (rMHD) and 3-D presureless two-fluid MHD models of the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI). HIT-SI is a current drive experiment that uses two geometrically asymmetric helicity injectors to generate and sustain toroidal plasmas. The comparable size of the collisionless ion skin depth d{sub i} to the resistive skin depth predicates the importance of the Hall term for HIT-SI. The simulations are run with NIMROD, an initial-value, 3-D extended MHD code. The modeled plasma density and temperature are assumed uniform and constant. The helicity injectors are modeled as oscillating normal magnetic and parallel electric field boundary conditions. The simulations use parameters that closely match those of the experiment. The simulation output is compared to the formation time, plasma current, and internal and surface magnetic fields. Results of the study indicate 2fl-MHD shows quantitative agreement with the experiment while rMHD only captures the qualitative features. The validity of each model is assessed based on how accurately it reproduces the global quantities as well as the temporal and spatial dependence of the measured magnetic fields. 2fl-MHD produces the current amplification (I{sub tor}/I{sub inj}) and formation time τ{sub f} demonstrated by HIT-SI with similar internal magnetic fields. rMHD underestimates (I{sub tor}/I{sub inj}) and exhibits much a longer τ{sub f}. Biorthogonal decomposition (BD), a powerful mathematical tool for reducing large data sets, is employed to quantify how well the simulations reproduce the measured surface magnetic fields without resorting to a probe-by-probe comparison. BD shows that 2fl-MHD captures the dominant surface magnetic structures and the temporal behavior of these features better than rMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Northwest Geysers EGS Demonstration Project, California. Pre-stimulation Modeling and Interpretation of the Stimulation

The Northwest Geysers Enhanced Geothermal System (EGS) demonstration project aims to create an EGS by directly and systematically injecting cool water at relatively low pressure into a known High Temperature (280–400 °C) Zone (HTZ) located under the conventional (240 °C) geothermal steam reservoir at The Geysers geothermal field in California. Here we report that , the results of coupled thermal, hydraulic, and mechanical (THM) analyses made using a model developed as part of the pre-stimulation phase of the EGS demonstration project is presented. The model simulations were conducted in order to investigate injection strategies and the resulting effects of cold-water injection upon the EGS system; in particular to predict the extent of the stimulation zone for a given injection schedule. The actual injection began on October 6, 2011, and in this paper a comparison of pre-stimulation model predictions with micro-earthquake (MEQ) monitoring data over the first few months of a one-year injection program is presented. The results show that, by using a calibrated THM model based on historic injection and MEQ data at a nearby well, the predicted extent of the stimulation zone (defined as a zone of high MEQ density around the injection well) compares well with observed seismicity. The modeling indicates that the MEQ events are related to shear reactivation of preexisting fractures, which is triggered by the combined effects of injection-induced cooling around the injection well and small changes in steam pressure as far as half a kilometer away from the injection well. Pressure-monitoring data at adjacent wells and satellite-based ground-surface deformation data were also used to validate and further calibrate reservoir-scale hydraulic and mechanical model properties. The pressure signature monitored from the start of the injection was particularly useful for a precise back-calculation of reservoir porosity. Ultimately, the first few months of reservoir pressure and surface deformation data were useful for estimating the reservoir-rock permeability and elastic modulus. Finally, although the extent of the calculated stimulation zone matches the field observations over the first few months of injection, the observed surface deformations and MEQ evolution showed more heterogeneous behavior as a result of more complex geology, including minor faults and fracture zones that are important for consideration in the analysis of energy production and the long-term evolution of the EGS system.

58 GEOSCIENCES↗

Inherent safety of minimum-burnup breed and burn reactors

Reactors that aim to sustain the breed and burn (B and B) mode of operation at minimum discharge burnup require excellent neutron economy, Minimum-burnup B and B cores are generally large and feature low neutron leakage probability and a hard neutron spectrum. While highly promising fuel cycles can be achieved with such designs, the very same features are pushing the limits of the core's ability to passively respond safely to unprotected accidents. Low leakage minimum-burnup sodium-cooled B and B cores have a large positive coolant void-worth and coolant temperature reactivity coefficient. In this study, the applicability of major approaches for fast reactor void-worth reduction is evaluated specifically for B and B cores. The design, shuffling scheme and performance of a new metallic-fueled, sodium-cooled minimum burnup B and B core, used as basis for the void-worth reduction analysis, is presented. The analysis shows that reactivity control systems based on passive {sup 6}Li injection during temperature excursions are the only option able to provide negative void-worth without significantly increasing the minimum burnup required for sustaining the B and B mode of operation. A new type of lithium expansion module (LEM) system was developed specifically for B and B cores and its effect on core performance is presented. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗

Foot-scale Evaluation of CO2-Responsive Polymer and CO2 Binary Fluid as an Alternative Fracturing Fluid for Enhanced Geothermal Systems

A foot-scale high-temperature true-triaxial fracturing apparatus was used to conduct hydraulic fracture tests on rock samples. The fracturing tests used different fluids including water, CO2, CO2 with water, and CO2 with a 1% polyallylamine (PAA)/water solution. The tests were conducted with both constant pressure injection and constant flow rate injection modes. CO2-based fracturing fluids were found to produce higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water-based fracturing fluids. 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 fracturing with CO2-based fluids in the presence of a PAA aqueous solution, for example if the rock is saturated with 1wt% PAA in water, the volume expansion caused by CO2-induced crosslinking of PAA leads to a faster pressure increase due to the associated volume expansion 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. Finally, the well-known low viscosity of CO2 phase prevents the efficient transport of proppants while the reacted CO2/PAA fluid system is known to form a high viscosity binary fracturing fluid system with the potential to carry proppant in addition of being a fine fracturing fluid.

Enhanced Geothermal Systems, polyallylamine, Co2-r↗

Analysis of ICRH antenna loading data in TEXTOR obtained during gas injection experiments

The possibility of increasing the coupling of ICRH power to plasmas on TEXTOR by gas injection was investigated, for a given voltage applied at the input of an antenna pair. The antenna pair was operated in the conjugated T mode, D(H) heating was used. Modeling of the antenna by transmission line theory shows that load resilience is maintained in a broad range of independent resistance variation for each strap. It is found that the absolute value of the reflection coefficient can be maintained below 0.2 for typical values of resistance measured in plasma shots (2-10 {Omega}/m) in good matching conditions. During gas injection, the loading resistance showed a clear increase with increasing line average electron density measured close to the plasma edge. Evidence of fast wave eigenmodes was also found, as characteristic resonant behavior of loading resistance and antenna self-inductance, due to poor absorption in the plasma caused by high H minority concentration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Revisiting the Analytical Solutions of Heat Transport in Fractured Reservoirs Using a Generalized Multirate Memory Function

Numerous analytical solutions have been developed for modeling thermal perturbations to underground formations caused by deep-well injection of fluids. Each solution has been derived for a specific boundary value problem and a simplified flow network with one set of parallel fractures. In this paper, new generalized solutions G*( x, s ) are developed using (existing) global transfer functions G$^*_0(x,s)$ and a new memory function g*( s ), where x and s are the space and Laplace variable. The memory function represents the solutions of conductive heat exchange between fractures and matrix blocks and between fractured aquifers and unfractured aquitards. The memory function is developed to account for multirate exchange induced by different shapes, sizes, properties, and volumetric fractions of matrix blocks bounded by multiple sets of orthogonal fractures with different spacing. The global transfer functions represent the fundamental solutions to convective, convective-conductive, and convective-dispersive heat transport in fractures (or aquifers) without exchange and are available for various (1-D linear, 1-D radial, 2-D dipole, and single-well injection-withdrawal) flow fields. The new solutions with exchange are developed using $G^*(x,s) = B^*(s)G^*_0{x, s[1 + Θg^*(s)]}$, thereby greatly simplifying solution development in a novel way, where Θ and B*( s ) are a fracture-matrix scaling factor and the boundary condition function. As a result, the new solutions are applied to several example problems, showing that heat transport in fractured aquifers is significantly impacted by (1) thermal dispersion in fractures that is rarely considered, (2) multirate heat exchange with a wide range of size and anisotropy of rectangular matrix blocks, and (3) heat exchange between aquifers and aquitards.

58 GEOSCIENCES↗

3-D Modeling of Induced Seismicity Along Multiple Faults: Magnitude, Rate, and Location in a Poroelasticity System

Understanding of the potential to injection–induced seismicity along faults requires the response of fault zone system to spatiotemporal perturbations in pore pressure and stress. In this study, three–dimensional (3–D) model system consisting of the caprock, reservoir, and basement is intersected by vertical strike–slip faults. We examine the full poroelastic behavior of the formation and perform the mechanical analysis along each fault zone using the Coulomb stress change. The magnitude, rate, and location of potential earthquakes are predicted using the spatial distribution of stresses and pore pressure over time. Rapid diffusion of pore pressure into conductive faults initiates failure, but the majority of induced seismicity occurs at deep fault zones due to poroelastic stabilization near the injection interval. Less permeable faults can be destabilized by either delayed pore pressure diffusion or poroelastic stressing. A two–dimensional (2–D) horizontal model, representing the interface between the reservoir and the basement, limits diffusion of pore pressure and deformation of the formation in the vertical direction that may overestimate or underestimate the potential of earthquakes along the fault. Lastly, our numerical results suggest that the 3–D modeling of faulting system including poroelastic coupling can reduce the uncertainty in the seismic hazard prediction by considering the hydraulic and mechanical interaction between faults and bounding formations.

3‐D modeling↗

Reaction Mechanism and Kinetics of Boron Removal from Molten Silicon via CaO-SiO{sub 2}-CaCl{sub 2} Slag Treatment and Ammonia Injection

To improve the boron-removal efficiency of metallurgical-grade silicon by increasing the reaction rate, a combined method with the 30 mol pct CaO-23.3 mol pct SiO{sub 2}-46.7 mol pct CaCl{sub 2} slag treatment and ammonia injection at 1723 K to 1823 K was proposed. For 1 hour and at 1823 K, the maximum removal efficiency of boron was 98 pct, and the final boron concentration in silicon decreased to 1.5 ppmw by the present method without the introduction of the iron catalyst. A kinetic model was also established to clarify the reaction mechanism and rate-limiting steps of this complicated boron-removal process. In this model, the rate-limiting step is the mass transfer of boron oxide at the interface between the slag and silicon phase.

36 MATERIALS SCIENCE↗

Silicate Carbonation in Supercritical CO2 Containing Dissolved H2O: An in situ High Pressure X-Ray Diffraction Study

Technological advances have been significant in recent years for managing environmentally harmful emissions (mostly CO2) resulting from combustion of fossil fuels. Deep underground geologic formations are emerging as reasonable options for long term storage of CO2 but mechanisms controlling rock and mineral stability in contact with injected supercritical fluids containing water are relatively unknown. In this paper, we discuss mineral transformation reactions occurring between supercritical CO2 containing water and the silicate minerals forsterite (Mg2SiO4), wollastonite (CaSiO3), and enstatite (MgSiO3). This study utilizes newly developed in situ high pressure x-ray diffraction (HXRD) and in situ infra red (IR) to examine mineral transformation reactions. Forsterite and enstatite were selected as they are important minerals present in igneous and mafic rocks and have been the subject of a large number of aqueous dissolution studies that can be compared with non-aqueous fluid tests in this study. Wollastonite, classified as a pyroxenoid (similar to a pyroxene), was chosen as a suitably fast reacting proxy for examining silicate carbonation processes associated with a wet scCO2 fluid as related to geologic carbon sequestration. The experiments were conducted under modest pressures (90 to 160 bar), temperatures between 35° to 70° C, and varying concentrations of dissolved water. Under these conditions scCO2 contains up to 3,500 ppm dissolved water.

carbon sequestration↗

Laboratory evaluation of cyclic underground hydrogen storage in the Temblor sandstone of the San Joaquin Basin, California

Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.

08 HYDROGEN↗

Stream lamination and rapid mixing in a microfluidic jet for X-ray spectroscopy studies

Microfluidic mixers offer new possibilities for the study of fast reaction kinetics down to the microsecond time scale, and methods such as soft X-ray absorption spectroscopy are powerful analysis techniques. These systems impose challenging constraints on mixing time scales, sample volume, detection region size and component materials. The current work presents a novel micromixer and jet device which aims to address these limitations. The system uses a so-called ‘theta’ mixer consisting of two sintered and fused glass capillaries. Sample and carrier fluids are injected separately into the inlets of the adjacent capillaries. At the downstream end, the two streams exit two micron-scale adjoining nozzles and form a single free-standing jet. The flow-rate difference between the two streams results in the rapid acceleration and lamination of the sample stream. This creates a small transverse dimension and induces diffusive mixing of the sample and carrier stream solutions within a time scale of 0.9 microseconds. The reaction occurs at or very near a free surface so that reactants and products are more directly accessible to interrogation using soft X-ray. We use a simple diffusion model and quantitative measurements of fluorescence quenching (of fluorescein with potassium iodide) to characterize the mixing dynamics across flow-rate ratios.

Huyke, Diego A. (ORCID:0000000283356613)↗

Modifications of Carbonate Fracture Hydrodynamic Properties by CO 2 -Acidified Brine Flow

Acidic reactive flow in fractures is relevant in subsurface activities such as CO{sub 2} geological storage and hydraulic fracturing. Understanding reaction-induced changes in fracture hydrodynamic properties is essential for predicting subsurface flows such as leakage, injectability, and fluid production. In this study, x-ray computed tomography scans of a fractured carbonate caprock were used to create three dimensional reconstructions of the fracture before and after reaction with CO{sub 2}-acidified brine (Ellis et al., 2011, Greenhouse Gases: Sci. Technol., 1:248-260). As expected, mechanical apertures were found to increase substantially, doubling and even tripling in some places. However, the surface geometry evolved in complex ways including ‘comb-tooth’ structures created from preferential dissolution of calcite in transverse sedimentary bands, and the creation of degraded zones, i.e. porous calcite-depleted areas on reacted fracture surfaces. These geometric alterations resulted in increased fracture roughness, as measured by surface Z{sub 2} parameters and fractal dimensions D{sub f}. Computational fluid dynamics (CFD) simulations were conducted to quantify the changes in hydraulic aperture, fracture transmissivity and permeability. The results show that the effective hydraulic apertures are smaller than the mechanical apertures, and the changes in hydraulic apertures are nonlinear. Overestimation of flow rate by a factor of two or more would be introduced if fracture hydrodynamic properties were based on mechanical apertures, or if hydraulic aperture is assumed to change proportionally with mechanical aperture. The differences can be attributed, in part, to the increase in roughness after reaction, and is likely affected by contiguous transverse sedimentary features. Hydraulic apertures estimated by the 1D statistical model and 2D local cubic law (LCL) model are consistently larger than those calculated from the CFD simulations. In addition, a novel ternary segmentation method was devised to handle the degraded zones, allowing for a bounding analysis of the effects on hydraulic properties. We found that the degraded zones account for less than 15% of the fracture volume, but cover 70% to 80% of the fracture surface. When the degraded zones are treated as part of the fracture, the fracture transmissivities are two to four times larger because the fracture surfaces after reaction are not as rough as they would be if one considers the degraded zone as part of the rock. Therefore, while degraded zones created during geochemical reactions may not significantly increase mechanical aperture, this type of feature cannot be ignored and should be treated with prudence when predicting fracture hydrodynamic properties.

roughness↗

Forsterite [Mg 2 SiO 4 )] Carbonation in Wet Supercritical CO 2 : An in Situ High-Pressure X-ray Diffraction Study

Technological advances have been significant in recent years for managing environmentally harmful emissions (mostly CO2) resulting from combustion of fossil fuels. Deep underground geologic formations are emerging as reasonable options for long term storage of CO2 but mechanisms controlling rock and mineral stability in contact with injected supercritical fluids containing water are relatively unknown. In this paper, we discuss mineral transformation reactions occurring with forsterite (Mg2SiO4) exposed to wet supercritical CO2. Forsterite was selected as it is an important olivine group mineral present in igneous and mafic rocks and has been the subject of a large number of aqueous dissolution studies that can be compared with non-aqueous fluid tests in this study. Transformation reactions were examined by in situ high pressure x-ray diffraction in the presence of supercritical carbon dioxide (scCO2) containing dissolved water at conditions relevant to carbon sequestration. Under modest pressures (90 bar) and temperatures (50°C), scCO2 saturated with water was found to convert >70 wt% forsterite to a hydrated magnesium carbonate, nesquehonite (MgCO3 •3H2O) and magnesite (MgCO3), after 72 hours of reaction. However, comparable tests with scCO2 at only partial water saturation (82%) showed a significantly slower carbonation rate with only ~30-39 wt% conversion to nesquehonite and no evidence of the anhydrous form (MgCO3). Further decreases in water content of the scCO2 continued to reduce the extent of carbonation, until a critical moisture threshold (~30%) was crossed where forsterite no longer reacted in the presence of the wet scCO2 to form crystalline carbonates. Increasing the temperature to 75°C produced anhydrous magnesium carbonate, magnesite (MgCO3), preceded by the intermediate phase, hydromagnesite [Mg(CO3)4(OH)2 •4H2O]. Measurements conducted during in situ IR experiments at 50°C and 30% saturation identified the presence of an amorphous carbonate phase as well as the formation of a thin liquid-like water layer on the forsterite surface. The presence of this water film appears to be critical for the mineral carbonation of forsterite exposed to water bearing scCO2. In contrast, our prior studies with the mineral brucite [Mg(OH)2] showed extensive carbonation in the absence of a condensed water layer on the mineral surface. The contrasts in reaction rate and products formed demonstrated by temperature and water-content dependence highlights the importance of these kinds of studies to help enable better predictions of the long term fate of geologically stored CO2.

54 ENVIRONMENTAL SCIENCES↗