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 199 records · Page 11

Machine learning-assisted upscaling analysis of reservoir rock core properties based on micro-computed tomography imagery

Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.

58 GEOSCIENCES↗

Mineralogical Controls on Frictional Strength, Stability, and Shear Permeability Evolution of Fractures

Massive fluid injection into the subsurface can induce microearthquakes by reactivating preexisting faults or fractures as seismic or aseismic slip. Such seismic or aseismic shear deformationsmay result in different modes of permeability evolution. Previous experimental studies have explored frictional stability-permeability relationships of carbonate-rich and phyllosilicate-rich samples under shear, suggesting that friction-permeability relationship may be primarily controlled by fracture minerals. We examine this relationship and identify the role of mineralogy (i.e., tectosilicate, carbonate, and phyllosilicate content) using direct-shear experiments on smooth saw-cut fractures of natural rocks and sintered fractures with distinct mineralogical compositions. These results indicate that the friction-permeability relationship is controlled by mineralogy. Frictional strength and permeability change upon reactivation decrease with phyllosilicate content but increase with tectosilicate content. In contrast, the reverse trend is observed for frictional stability (a-b). However, the permeability change decreases with carbonate content while both frictional strength and stability increase. The permeability change always decreases with an increase in frictional stability. This relationship implies a new mechanical-hydrochemical coupling loop via a linkage of frictional properties, mineralogy, and permeability.

58 GEOSCIENCES↗

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗

Hydromechanical Modeling of Fault Rupture in Geologic CO 2 Sequestration: A Comparison of Two Failure Criteria

Abstract The large‐scale implementation of geological carbon sequestration has raised concerns about potential fault activation and induced seismicity, which could compromise storage integrity and pose seismic risks. We theoretically compared two failure criteria, Mohr‐Coulomb (MC) and Modified Cam‐Clay criteria (MCC), to assess fault rupture during CO 2 storage. Both criteria characterize fault behavior in a specific stress regime but differ in reducing the complexity of fault rupture to a few key mechanisms. Using a coupled hydromechanical model, we demonstrate that the choice of a failure criterion and the physics of fault weakening associated with these criteria strongly condition fault response to a given fluid injection. MC mainly relates rupture to friction, while MCC relates rupture to fault poro‐plasticity. Our findings highlight that the selection of a failure criterion, being inherently subjective, can significantly alter the predicted fault behavior during CO 2 storage, thereby impacting the reliability of geomechanical risk assessments.

Cao, M. [Energy Geosciences Division Lawrence Berk↗

Heliostat Aiming Strategies for Arbitrary Flux Profiles for High Temperature Gas-Based Receivers

High temperature gas-based receivers are seen as one potential path toward achieving the next generation of concentrated solar power tower systems. A key challenge associated with high temperature, high pressure heat transfer fluids is the lifetime of receiver containment materials. Achieving the desired lifetime while also designing an efficient receiver requires implementation of unique flux profiles, tailored to place high fluxes at cold fluid injection points. An aimpoint selection methodology was developed in the present study to allow user-defined flux profiles to be prescribed for a receiver. The methodology was implemented in NREL's open-source software SolarPILOT, modifying the existing “image size priority” aimpoint strategy. A case study for a tubular gas-based receiver is presented, examining the effect of various parameters on the ability to replicate the desired flux profile.

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

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗

Natural time and nowcasting induced seismicity at the Groningen gas field in the Netherlands

The Groningen gas field in the Netherlands has a significant level of induced seismicity; this seismicity is attributed to fluid withdrawal not fluid injection. This is one of the most productive gas fields in Europe. Production began in 1963 and the first induced earthquake occurred in 1991. Increasing levels of seismicity led to the installation of a sensitive seismic network in 1995 with a minimum magnitude completeness of $M$ ≥ 1.5. The largest induced earthquake was a magnitude 3.6 in 2012, which caused thousands of damage claims and led to a reduction in gas production in 2014 January. In our work, we show that the induced earthquakes satisfy Gutenberg–Richter (GR) frequency–magnitude statistics over the range 1.5 ≤ $M$ ≤ 3.0. We utilize truncated GR scaling to introduce a possible maximum earthquake with $M$ = 3.85. We also demonstrate an almost complete absence of aftershocks following the largest induced earthquakes. To quantify the temporal evolution of the induced seismicity, we perform our nowcasting method. Nowcasting uses the catalogue of seismicity in the region. Two earthquake magnitudes are selected, one large say $M_λ$ ≥ 2.5, and one small say $M_σ$ ≥ 1.5. The method utilizes the number of small earthquakes that occur between pairs of large earthquakes. A major advantage of nowcasting is that it relies on ‘natural time’, earthquake counts, between events rather than clock time. Thus, the results are applicable to induced seismicity that varies in time. We count the rate of occurrence of small induced earthquakes to nowcast the probability of occurrence of larger earthquakes.

58 GEOSCIENCES↗

Semi-Analytical Leakage Solutions for Aquifers (SALSA) v1

SALSA is a semi-analytical modeling tool that can provide assessments of pressure perturbations and brine leakage because of fluid injection and extraction activities in sedimentary basins. Sedimentary basins typically contain vertical sequences of near-horizontal aquifers separated by less-permeable aquitards. In the context of geologic carbon sequestration (GCS), while aquitards are relied upon to contain injected buoyant free-phase CO2 and limit inter-aquifer brine flow, the large number of leaky wells (e.g., unused/unsuccessful exploration wells, or improperly plugged water or oil/gas wells) that are present in sedimentary basins creates concerns for potential leakage and groundwater contamination. The mathematical theory and solution method included in SALSA is presented in our recent manuscript under review. SALSA is the first to account for the coupled leakage through aquitards and leaky wells, with geologic pressurization sources. SALSA based on semi-analytical solutions can be computationally very efficient for problems containing many injection and leaky wells in multilayered sedimentary systems, compared to any numerical simulation software requiring local mesh refinement around each well to obtain accurate results. Therefore, stakeholders (e.g., regulators, operators) for site screening, injection and post-injection pressure behavior and assessing leakage risks can use SALSA as a fast-predictive tool in GCS applications in multilayered aquifer systems.

Cihan, Abdullah↗

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

DTS Raw Data Guelph, Ontario Canada

Unprocessed active distributed temperature sensing (DTS) data from 3 boreholes in the Guelph, ON Canada region. Data from borehole 1 was collected during a fluid injection while data from boreholes 2 and 3 were collected under natural gradient conditions in a lined borehole. The column labels/headers (in the first row) define the time since start of measurement in seconds and the row labels/headers (in the first column) are the object IDs that are defined in the metadata. Each object ID is a sampling location whose exact location is defined in the metadata file. Data in each cell are temperature in Celsius at time and sampling location as defined above.

15 GEOTHERMAL ENERGY↗

Geothermal Reservoir Simulation Results in support of Feasibility Study of Direct District Heating for the Cornell Campus Utilizing Deep Geothermal Energy

This dataset contains input data, code, ReadMe files, output data, and figures that summarize the results of a stochastic analysis of geothermal reservoir production from two potential geothermal reservoirs that were evaluated for the Cornell University Deep Direct-Use project. These potential reservoirs are the Trenton-Black River (TBR) from 2.27-2.3 km depth, and basement rocks from 3.0-3.5 km depth and 3.5-4.0 km depth. Several utilization scenarios consisting of different injection fluid temperatures and flow rates were evaluated for each reservoir. Uncertainty in geologic properties, thermal properties, economic costs, and utilization efficiencies were evaluated using a Monte Carlo analysis of the reservoir simulations. Some reservoir simulations of the TBR were completed using the TOUGH2 software, as implemented in PetraSIM. The PetraSIM run files and associated data are provided with this submission. All other reservoir simulations were completed using the GEOPHIRES software, with some modifications to complete the uncertainty analyses. ReadMe files that describe additions to GEOPHIRES, the GEOPHIRES input data, and the output data are all provided, and references are provided to the code repository. Figures that summarize the reservoir heat production, temperature drawdown, and the probability of meeting targeted building heating demands with the produced heat and fluid temperatures are provided.

15 GEOTHERMAL ENERGY↗

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Large Scale Grid

Mesh, properties, initial conditions, injection/withdrawal rates for modeling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, on large scale grid, with some results. First simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 1-5 pertain.

15 GEOTHERMAL ENERGY↗

EGS Collab Experiment 2: Continuous Active Source Seismic Monitoring (CASSM)

The dataset contains continuous active-source seismic monitoring (CASSM) data collected during EGS Collab Experiment 2, conducted from February to September 2022 at the Sanford Underground Research Facility in Lead, South Dakota. This experiment aimed to investigate enhanced geothermal systems through high-pressure fluid injections at depths of 1200-1500 meters. The seismic monitoring system included 16 three-component piezoelectric accelerometers and 24 hydrophones installed in boreholes around the injection zones, recording signals from piezoelectric seismic sources. Data were acquired using both continuous and triggered recording systems, with sampling rates of up to 100 kHz. The raw data are organized by timestamps and stored in .dat format, with accompanying log files. Calibration certificates for selected accelerometers are provided to aid in correcting sensor responses, though users are advised to consider possible effects of enclosures and installation on sensor performance. Users are strongly advised to consult the accompanying report, which outlines the experimental setup, data acquisition, sensor specifications, and recording systems.

15 GEOTHERMAL ENERGY↗

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗

Locating Microearthquakes In The Salton Sea Geothermal Field, California

We report the results of a 6-month microearthquake study of the Salton Sea Geothermal Field, California. This active geothermal field is also seismically active due to the flow of fluids through pre-existing fractures of rocks in the subsurface which can cause faults to slip. Locating these seismic events within the geothermal field enables researchers to determine whether the earthquakes are anthropogenic or tectonic from the nearby Brawley Seismic Zone. We analyzed November 2009 to April 2010 seismic waveform data retrieved from a network of continuously recording seismometers covering an area of 10 km 2 . 2,271 microearthquake hypocenters were processed out of 2,609 total recorded events and their locations suggest the majority could be at approximately the same depth as the fluid injection wells.

58 GEOSCIENCES↗

Mapping Diffuse Seismicity for Geothermal Reservoir Management with Matched Field Processing: Year 2 Report

Detecting and locating seismic events in managed underground geothermal reservoirs can become increasingly difficult during certain reservoir management operations. For example, during fluid injection procedures, large numbers of microearthquakes frequently occur during a relatively short amount of time. On seismic records, these events can overlap or have poor signal-to-noise ratios (SNR). This can obscure the onset of smaller signals and complicate the picking of seismic phases needed for earthquake location algorithms. To accurately illuminate the active areas of the underground reservoir through time however, a complete earthquake catalog down to very small magnitudes is necessary. To aid in the seismic characterization of reservoir fracture networks, we propose to complement traditional earthquake detection and location techniques with the empirical matched field processing (MFP) method. The principle of the empirical MFP method is to calculate propagation from source to receivers for pre-existing events, also called master templates. Then MFP matches the spatial structure of incoming seismicity observed by a network of sensors to master templates keyed to potential event locations. We demonstrate that empirical MFP can complement existing catalogs and techniques by using archived seismic information to increase earthquake catalog completeness. We identify and construct representative master templates using the Southern California Earthquake Data Center (SCEDC) earthquake catalog and archived waveform data of events originating within the Salton Sea Geothermal Field (SSGF). We identify 231 master events that have good quality seismic records at four or more seismic stations. We apply the empirical MFP method to continuous seismic data collected at the SSGF between November 2009 and December 2010. The MFP method successfully identified 5357 local events, while the original catalog only reported 1536 events. This increase in the number of events allows us to better determine the active fractures of the geothermal reservoir. We model the Coulomb stress change due to events between M2 – M3.5 that occurred during the January 2010 seismic swarm. Our results show a good correlation between regions of increased Coulomb stress change due to double-couple shear events on faults and the locations of subsequent microseismicity. In addition, we investigate the lower event detection magnitude threshold of our surface network using the MFP method. It clearly shows the event detection threshold is extended to include many smaller earthquakes with magnitudes down to M0.0. Therefore, we believe that the empirical MFP method, when combined with conventional methods, significantly improves network detection capabilities.

15 GEOTHERMAL ENERGY↗

High Fidelity Computational Analysis of CO2 Trapping at Pore Scales (Final Technical Report)

With an alarming rise in carbon dioxide (CO2) emission from anthropogenic sources, CO2 sequestration has become an attractive choice to mitigate the emission. Some popular storage media for CO 2 are oil reservoirs, deep coal-bed, and deep oceanic-beds. These have been used for the long term CO 2 storage. Due to special lowering viscosity and surface tension property of CO 2 , it has been widely used for enhanced oil recovery. The sites for CO 2 sequestration or enhanced oil recovery mostly consist of porous rocks. Lack of knowledge of molecular mobility under confinement and molecule-surface interactions between CO2 and natural porous media results in generally governed by unpredictable absorption kinetics and total absorption capacity for injected fluids, and therefore, constitutes barriers to the deployment of this technology. Therefore, it is important to understand the flow dynamics of CO 2 through the porous microstructures at the finest scale (pore-scale) to accurately predict the storage potential and long-term dynamics of the sequestered CO 2 . This report discusses about pore-network flow modeling approach using variational method and analyzes simulated results this method simulations at pore-scales for idealized network and using Berea Sandstone CT scanned images. Variational method provides a promising way to study the kinetic behavior and storage potential at the pore scale in the presence of other phases. The current study validates variational solutions for single and two-phase Newtonian and single phase non-Newtonian flow through angular pores for special geometries whose analytical and/or empirical solutions are known. The hydraulic conductance for single phase flow through a triangular duct was also validated against empirical results derived from lubricant theory.

54 ENVIRONMENTAL SCIENCES↗

A microseismic workflow for managing induced seismicity risk as CO 2 storage projects

It is well established that fluid injection has the potential to induce earthquakes—from microseismicity to large, damaging events—by altering state-of-stress conditions in the subsurface. While induced seismicity has not been a major operational issue for carbon storage projects to date, a seismicity hazard exists and must be carefully addressed. Two essential components of effective seismic risk management are (1) sensitive microseismic monitoring and (2) robust data interpretation tools. This report describes a novel workflow, based on advanced processing algorithms applied to microseismic data, to help improve management of seismic risk. This workflow has three main goals: (1) to improve the resolution and reliability of passive seismic monitoring, (2) to extract additional, valuable information from continuous waveform data that is often ignored in standard processing, and (3) to minimize the turn-around time between data collection, interpretation, and decision-making. These three objectives can allow for a better-informed and rapid response to changing subsurface conditions.

58 GEOSCIENCES↗