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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular-scale mechanisms of CO 2 mineralization in nanoscale interfacial water films

The calamitous impacts of untethered carbon emission from fossil-fuel-burning energy infrastructure calls for accelerated development of large-scale CO 2 capture, utilization, and storage technologies that are underpinned by fundamental understanding of molecular-level chemical processes. In the subsurface, rocks rich in divalent metals can react with CO 2 , permanently sequestering it in the form of stable metal carbonate minerals, with the CO 2 -H 2 O composition of the post-injection pore fluid acting as a primary control variable. Herein, we compare mechanistic reaction pathways for aqueous-mediated carbonation with carbon mineralization occurring in nanoscale adsorbed water films. In the extreme of pores filled with a CO 2 -dominant fluid, carbonation reactions are confined to nanometer-thick water films coating mineral surfaces, which enable metal cation release, transport, nucleation, and crystallization of metal carbonate minerals. Though seemingly counterintuitive, laboratory studies have demonstrated facile carbonation rates in these low water environments, which in recent years has begun to be better understood in mechanistic detail. The overarching objective of this review is to delineate the unique underlying molecular-scale reaction mechanisms that govern CO 2 mineralization in these reactive and dynamic quasi-2D interfaces. Here. we highlight the importance of understanding unique properties in thin water films, such as how water dielectric properties, and consequently ion solvation/hydration behavior, can change under nanoconfinement. We conclude by identifying important frontiers for future work and opportunities to exploit these fundamental chemical insights for decarbonization technologies in the 21st century.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Study of Temperature Distribution and Thermal Stresses in a Hot Rock Due to Rapid Cooling

Thermal stresses may be induced in a hot dry rock when a cold fluid is injected in the well. To study this problem, we look at the thermoelastic response of a hot rock that is suddenly cooled. The cooling is assumed to be either at a constant temperature or at a constant heat flux per unit depth. Our approach is to nondimensionalize the equations and perform a parametric study and look at the temperature distribution and the induced-thermal stresses. The results indicate that depending on the extent of cooling and the cooling time, thermal stresses can be induced. Numerical simulations on sandstone, with an initial uniform temperature of 473 K, are also carried out. The results show that if the cooling is due to the surface temperature maintained at 463 K (10 °C lower than the initial temperature of the hot rock), thermal stresses that are larger than the rock tensile strength could be induced. When the cooling is due to a constant surface heat flux, this temperature can be reached after about 777 days of cooling with a minimum value of a heat flux of -20 W/m.

Tran, Phuoc↗
DOE OSTI2021

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗
DOE OSTI2022

EGS Collab Experiment 2: Continuous Broadband Seismic Waveform Data

Two broadband seismometers were installed on the 4100 level and recorded for the duration of EGS Collab Experiment #2. Inspired by published data from similar instruments installed in the Aspo Hard Rock Lab, these long-period instruments aimed to measure the tilting of the drift in response to the injection of fluid into the testbed. One instrument was installed underneath the wellheads in Site A (aka the "battery" alcove) and the other was installed along the east wall of the drift, south of Site B. Due to the feet of gravel (ballast) laid along the floor of the drift, we were unable to anchor the sensors directly to the rock. As a result, the coupling of the sensors to the experiment rock volume is likely poor. In addition, there are a number of noise sources that complicate the interpretation of the data. For example, sensor BBB is installed adjacent (within 3 ft) to the rail line that runs towards the Ross shaft. Trains (motors) run along this line almost daily and produce a large signal in these data. Careful extraction of periods of interest, as well as filtering for specific signals, is necessary. The sensors are Nanometrics Trillium Compact Posthole seismometers, sensitive down to 120 seconds period. They were installed as close to the drift wall and as deep as we could manually excavate (only about 1 ft or so). The holes were leveled with sand and the sensors were placed on a paver before backfilling with sand. The hole was then covered by a bucket filled with insulation to improve the sensor's isolation from daily temperature variations, which are minor but present due to drift ventilation from the surface. Data were recorded on Nanometrics Centaur digitizers at 100 Hz. The full response information is available in the StationXML file provided here, or by querying the sensors through the IRIS DMC (see links below). These instruments were provided free of charge through the IRIS PASSCAL instrument center. The network code is XP and the station codes are BBA and BBB. The waveform data can be queried through the IRIS FDSN server using any method the user likes. One convenient option is to use the Obspy python package: https://docs.obspy.org/packages/obspy.clients.fdsn.html

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Utah FORGE: Seismic DAS and Geophone Borehole Data Processing and 3D Imaging of Vp/Vs Ratio in the 2024 Stimulated Reservoir

This dataset includes a final report and a 3D velocity model derived from seismic DAS and geophone borehole data collected during the April 2024 stimulation of the reservoir at Utah FORGE. The report details the processing of over 50,000 P- and S-wave travel times used in a tomographic inversion to estimate the Vp/Vs ratio, revealing anomalies adjacent to wells 16A and 16B that may be associated with injected fracturing fluids. Additionally, Wadati analysis of more than 27,000 pairs of differential P- and S-wave travel times supports this interpretation with a high Vp/Vs estimate of 1.86, averaged over the dimensions of the earthquake cluster. The accompanying data file provides a 3D model of P- and S-wave velocities and Vp/Vs ratios, structured with spatial coordinates and velocity values.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A Project Lifetime Approach to the Management of Induced Seismicity Risk at Geologic Carbon Storage Sites

The geologic storage of carbon dioxide (CO 2 ) is one method that can help reduce atmospheric CO 2 by sequestering it into the subsurface. Large-scale deployment of geologic carbon storage, however, may be accompanied by induced seismicity. We present a project lifetime approach to address the induced seismicity risk at these geologic storage sites. This approach encompasses both technical and nontechnical stakeholder issues related to induced seismicity and spans the time period from the initial consideration phase to postclosure. These recommendations are envisioned to serve as general guidelines, setting expectations for operators, regulators, and the public. They contain a set of seven actionable focus areas, the purpose of which are to deal proactively with induced seismicity issues. Although each geologic carbon storage site will be unique and will require a custom approach, these general best practice recommendations can be used as a starting point to any site-specific plan for how to systematically evaluate, communicate about, and mitigate induced seismicity at a particular reservoir.

58 GEOSCIENCES↗
DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗
DOE OSTI2021

Deposition, Diagenesis, and Sequence Stratigraphy of the Pennsylvanian Morrowan and Atokan Intervals at Farnsworth Unit

Farnsworth Field Unit (FWU), a mature oilfield currently undergoing CO2-enhanced oil recovery (EOR) in the northeastern Texas panhandle, is the study area for an extensive project undertaken by the Southwest Regional Partnership on Carbon Sequestration (SWP). SWP is characterizing the field and monitoring and modeling injection and fluid flow processes with the intent of verifying storage of CO2 in a timeframe of 100–1000 years. Collection of a large set of data including logs, core, and 3D geophysical data has allowed us to build a detailed reservoir model that is well-grounded in observations from the field. This paper presents a geological description of the rocks comprising the reservoir that is a target for both oil production and CO2 storage, as well as the overlying units that make up the primary and secondary seals. Core descriptions and petrographic analyses were used to determine depositional setting, general lithofacies, and a diagenetic sequence for reservoir and caprock at FWU. The reservoir is in the Pennsylvanian-aged Morrow B sandstone, an incised valley fluvial deposit that is encased within marine shales. The Morrow B exhibits several lithofacies with distinct appearance as well as petrophysical characteristics. The lithofacies are typical of incised valley fluvial sequences and vary from a relatively coarse conglomerate base to an upper fine sandstone that grades into the overlying marine-dominated shales and mudstone/limestone cyclical sequences of the Thirteen Finger limestone. Observations ranging from field scale (seismic surveys, well logs) to microscopic (mercury porosimetry, petrographic microscopy, microprobe and isotope data) provide a rich set of data on which we have built our geological and reservoir models.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Flow assisted anti-fouling geometries for compact heat exchangers

A header for a heat exchanger and method for cleaning a heat exchanger in a loop without disconnecting loop components is provided. The header is in flow communication with the heat exchanger for distributing fluid through a plurality of adjacent channels. The header is connected between a main heat exchanger inlet nozzle and a channel flow distributor. A filter element is disposed within the header between the nozzle and channel flow distributor. Under normal operation, the filter element removes particulates and fouling material from the main flow stream before it enters the heat exchanger channels. During the cleaning process, fluid is injected on or through the filter element to remove particulates and fouling material through at least one outlet port. The header arrangement allows the filter element to be ‘cleaned in place’ without draining the system and disconnecting the heat exchanger or other components from the flow loop.

Carlson, Matthew David↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2022

Fluid mixing apparatus using liquid fuel and high- and low-pressure fluid streams

A fluid mixing apparatus includes mixing conduits that extend through a fluid plenum. The fluid plenum, which surrounds a first wall defining a main passage fluidly coupled to a low-pressure fluid source, is surrounded itself by a second wall defining a high-pressure plenum fluidly coupled to a high-pressure fluid source. An insulated tube disposed at the inlet of the first wall delivers a third fluid. The mixing conduits fluidly couple the high-pressure plenum to the main passage, where the high-pressure fluid is mixed with low-pressure fluid and the third fluid. Optionally, the fluid plenum may house a fourth fluid that is injected through injection holes in the mixing conduits. The fluid mixing apparatus may be used to mix one or more fuels with high- and low-pressure air in a gas turbine combustor. Alternately, the fluid mixing apparatus may mix a fluid with high- and low-pressure water streams.

Berry, Joseph Jonathan↗
DOE OSTI2021

Characterizing mineralization on low carbon steel exposed to aerated and degassed synthetic hydraulic fracture fluids

This study investigates the mineralization occurring at the interface between injected hydraulic fracturing fluids and steel production casing during the shut-in period of unconventional shale wells. Mild steel coupons were reacted with synthetic hydraulic fracturing fluids with varying base fluid and dissolved gas concentrations. Experiments were carried out in an autoclave system set at 50 °C and 2000 PSI for 2 and 14 days. Physical and mineralogical changes to the steel coupons’ surface were characterized using an FEI Quanta 600 field emission gun (FEG) scanning electron microscope (SEM) and a PANalytical XPert Pro diffractometer. Barite precipitation occurred ubiquitously during the first 48 hours in all experimental scenarios. Sulfate radicals created from degradation and reduction of ammonium persulfate breakers facilitate barite formation. Concomitant barite and corrosion related precipitates on the sample surfaces suggest iron hydroxy salts act as nucleation sites for mineral scale growth. These results suggest hydraulic fracturing fluids may have deleterious effects on the wellbore before production of the well begins.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2023

Fluid mixing apparatus using high- and low-pressure fluid streams

A fluid mixing apparatus includes mixing conduits that extend through a fluid plenum and that define injection holes therethrough. The fluid plenum, which surrounds a first wall defining a main passage fluidly coupled to a low-pressure fluid source, is surrounded itself by a second wall defining a high-pressure plenum fluidly coupled to a high-pressure fluid source. The mixing conduits fluidly couple the high-pressure plenum to the main passage, and the fluid from the fluid plenum is delivered with the high-pressure fluid to the main passage, where the fluids mix before being discharged from an outlet of the main passage. The fluid mixing apparatus may be used to mix one or more fuels with high- and low-pressure air in a gas turbine combustor. Alternately, the fluid mixing apparatus may mix a fluid with high- and low-pressure water streams.

Berry, Jonathan Dwight↗
DOE OSTI2020

Electrically conductive proppant and methods for detecting, locating and characterizing the electrically conductive proppant

Electrically conductive proppants and methods for detecting, locating, and characterizing same are provided. The electrically conductive proppant can include a substantially uniform coating of an electrically conductive material having a thickness of at least 500 nm. The method can include injecting a hydraulic fluid into a wellbore extending into a subterranean formation at a rate and pressure sufficient to open a fracture therein, injecting into the fracture a fluid containing the electrically conductive proppant, electrically energizing the earth at or near the fracture, and measuring three dimensional (x, y, and z) components of electric and magnetic field responses at a surface of the earth or in an adjacent wellbore.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Electrically conductive proppant and methods for detecting, locating and characterizing the electrically conductive proppant

Electrically conductive proppants and methods for detecting, locating, and characterizing same are provided. The electrically conductive proppant can include a substantially uniform coating of an electrically conductive material having a thickness of at least 500 nm. The method can include injecting a hydraulic fluid into a wellbore extending into a subterranean formation at a rate and pressure sufficient to open a fracture therein, injecting into the fracture a fluid containing the electrically conductive proppant, electrically energizing the earth at or near the fracture, and measuring three dimensional (x, y, and z) components of electric and magnetic field responses at a surface of the earth or in an adjacent wellbore.

Cannan, Chad↗
DOE OSTI2024

Cyclic Injection Leads to Larger and More Frequent Induced Earthquakes under Volume-Controlled Conditions

As carbon storage technologies advance globally, methods to understand and mitigate induced earthquakes become increasingly important. Although the physical processes that relate increased subsurface pore pressure changes to induced earthquakes have long been known, reliable methods to forecast and control induced seismic sequences remain elusive. Suggested reservoir engineering scenarios for mitigating induced earthquakes typically involve modulation of the injection rate. Some operators have implemented periodic shutdowns (i.e., effective cycling of injection rates) to allow reservoir pressures to equilibrate (e.g., Paradox Valley) or shut-in wells after the occurrence of an event of concern (e.g., Basel, Switzerland). Other proposed scenarios include altering injection rates, actively managing pressures through coproduction of fluids, and preinjection brine extraction. Here, in this work, we use 3D physics-based earthquake simulations to understand the effects of different injection scenarios on induced earthquake rates, maximum event magnitudes, and postinjection seismicity. For comparability, the modeled injection considers the same cumulative volume over the project’s operational life but varies the schedule and rates of fluid injected. Simulation results show that cyclic injection leads to more frequent and larger events than constant injection. Furthermore, with intermittent injection scenario, a significant number of events are shown to occur during pauses in injection, and the seismicity rate remains elevated for longer into the postinjection phase compared to the constant injection scenario.

58 GEOSCIENCES↗