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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2011

The Role of Low-Angle Extensional Tectonics, Flat Fracture Domains, and Gravity Slides in Hydrothermal and EGS Resources of the Western United States

The Steamboat Springs geothermal system provides the most dramatic example of subhorizontal thermal-fluid aquifers in crystalline rock in the Basin and Range, but this is by no means an isolated case. Similar but more diffuse subhorizontal permeability has been reported at Roosevelt Hot Springs and Cove-Fort Sulphurdale, Utah; and a km-scale gravity-slide block channels injectate at Dixie Valley, Nevada. During the course of this phase of the project 2543 reports including text, figures and large format enclosures, 1428 maps, and 698 well logs were scanned. The information is stored in a Microsoft Access Database on the Geothermal Server. Detailed geologic cross sections of the Desert Peak geothermal field were developed to identify the structural controls on the geothermal system and locate possible fluid flow paths. The results of this work were published by Lutz and others (2009, Appendix 1) in the Stanford Reservoir Engineering Conference Proceedings.

02 PETROLEUM↗

DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗

DOE OSTI2011

Intrinsic rotation generation in ELM-free H-mode plasmas in the DIII-D tokamak--Experimental observations

A detailed description is presented of the experiment reported in [S. H. Mueller et al., Phys. Rev. Lett. 106, 115001 (2011)], which reported the first measurements of fluid turbulent stresses in a tokamak H-mode pedestal. Mach probe measurements disclosed a narrow co-current rotation layer at the separatrix, which is also seen in some L-modes [J. A. Boedo et al., Phys. Plasmas 18, 032510 (2011)]. Independent evidence for the existence of the edge co-rotation layer is presented from main-ion rotation measurements by charge-exchange-recombination spectroscopy in comparable helium plasmas. The probe measurements are validated against density and electron temperature profiles from Thomson scattering and in terms of the measured turbulent particle transport, which is consistent with the global density rise. Non-diffusive non-convective angular momentum transport is required by two independent experimental observations: (1) A persistent dip in the rotation profile separates the edge layer from the evolving core region during intrinsic rotation development. (2) The rotation profiles with co- and counter-current neutral beam injection appear well described as the simple sum of a constant intrinsic part and the beam-driven part, also demonstrating the profile-independence of the intrinsic torque. Characteristics of the turbulent fluctuations composing the fluid turbulent stresses are discussed: Up to 0.5 cm inside the separatrix, the low amplitude of the Reynolds stress (<0.05 Nm of torque) is due to both a reduction of the fluctuation amplitudes at the peak of the edge co-rotation layer and weak correlations between the toroidal and radial velocity fluctuations. Further into the core, the correlations increase significantly up to a value of +0.75, resulting in an almost unidirectional character of the turbulent Reynolds stress, generating substantial counter-current torques up to -2 Nm. Additional mechanisms must be present to balance these torques and explain the co-current core-plasma spin-up at a rate of +0.3 Nm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Evaporation Sub-model Development for Volume of Fluid (eVOF) Method Applicable to Spray-Wall Interaction Including Film Characteristics with Validation at High Pressure and Temperature Conditions (Final Report)

Internal combustion engines have seen a great evolution over the last several decades through application of high pressure direct injection, multiple injections, and other technologies to reduced fuel consumption, NOx, and PM. Although combustion systems with advanced injection strategies have been studied extensively, there exists a significant fundamental knowledge gap on the fuel-spray interactions with the piston surface and chamber walls. Advanced computational codes validated with experimental techniques have to be developed for accurate representation of the drop impingement, fuel film formation, and vaporization. Current engine CFD spray models utilize a Lagrangian framework for modeling which lacks critical considerations of the physics pertaining to these interactions and thus requiring extensive parameterization, tuning and validation. The team from Michigan Technological University, University of Massachusetts Dartmouth, and Argonne National Laboratory is composed of experts in sprays, combustion, engines and CFD with a wide spectrum of knowledge including specific expertise in the area under consideration. In the proposed work, a VOF modeling approach has been adopted for the spray-wall interaction, film formation and spreading, and vaporization. With the inclusion of a vaporization submodel, a more predictive and accurate simulation of the spray-film was performed without extensive need of parameterization and tuning. Extensive experimentation of the spray-wall interaction under the range of conditions matching the thermodynamic and surface temperatures that occur in diesel and gasoline engines were conducted to validate the SWI submodels and for development of the evaporation sub-model, which has been implemented in the Converge software.

42 ENGINEERING↗

DOE OSTI2014

Geomechanical Evaluation of Thermal Impact of Injected CO 2 Temperature on a Geological Reservoir: Application to the FutureGen 2.0 Site

The impact of temperature variations of injected CO 2 on the mechanical integrity of a reservoir is a problem rarely addressed in the design of a CO 2 storage site. The geomechanical simulation of the FutureGen 2.0 storage site presented here takes into account the complete modeling of heat exchange between the environment and CO 2 during its transport in the pipeline and injection well before reaching the reservoir, as well as its interaction with the reservoir host rock. An ad-hoc program was developed to model CO 2 transport from the power plant to the reservoir and an approach coupling PNNL STOMP-CO 2 multiphase flow simulator and ABAQUS® has been developed for the reservoir model which is fully three-dimensional with four horizontal wells and variable layer thickness. The Mohr-Coulomb fracture criterion has been employed, where hydraulic fracture was predicted to occur at an integration point if the fluid pressure at the point exceeded the least compressive principal stress. Evaluation of the results shows that the fracture criterion has not been verified at any node and time step for the CO 2 temperature range predicted at the top of the injection zone.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2025

Investigation of fuel film formation and soot emissions in a GDI engine during cold-start with Split-injection strategies

This study investigates the impact of fuel-film formation on engine-out soot emissions in a gasoline direct injection (GDI) engine under cold-start conditions. Split-injection strategies were applied by varying the number of injections, injection duration, and total fuel quantity to affect wall wetting and control the average in-cylinder equivalence ratio. A combined experimental and numerical approach was employed to analyze fuel-film deposition, combustion efficiency, and engine-out soot and unburnt hydrocarbons (UHC) emissions. In particular, fuel film distribution estimated by means of non-reacting, 3-D, computational fluid dynamics (CFD) simulations, together with experimentally measured soot data, were used to investigate fuel film formation and its role in soot generation. In the experiments, a skip-fired engine control strategy was applied to mimic the transient nature of engine cold-start operation. The results indicate that, under the same total number of injections, increasing the average in-cylinder equivalence ratio through longer injection durations improves combustion stability, as indicated by the decrease of the coefficient of variation of IMEP n (Net Indicated Mean Effective Pressure) from 6.0 % to 0.6 %. However, this strategy leads to higher soot emissions, which increased by nearly an order of magnitude, primarily due to enhanced wall-film formation. In contrast, increasing the number of injections while maintaining a constant equivalence ratio significantly impacts fuel-film deposition and, consequently, soot emissions, with a fivefold reduction of the measured engine-out soot, decreasing from 3.5 mg to 0.7 mg. A soot-film correlation was developed and achieved a high coefficient of determination (r 2 = 0.95) and was further extended to account for spark timing effects. These findings confirm the effectiveness of split-injection for avoiding wall film formation and soot emissions, and the critical role of fuel film in soot generation, supporting the hypothesis that pool fires play a crucial role in contributing to soot formation under these cold-start conditions. In conclusion, the study also indicates the value of a predictive soot-film correlation for developing cold-start emission control strategies.

Engine Cold-start↗

DOE OSTI2017

Geologic influence on induced seismicity: Constraints from potential field data in Oklahoma

Abstract Recent Oklahoma seismicity shows a regional correlation with increased wastewater injection activity, but local variations suggest that some areas are more likely to exhibit induced seismicity than others. We combine geophysical and drill hole data to map subsurface geologic features in the crystalline basement, where most earthquakes are occurring, and examine probable contributing factors. We find that most earthquakes are located where the crystalline basement is likely composed of fractured intrusive or metamorphic rock. Areas with extrusive rock or thick (>4 km) sedimentary cover exhibit little seismicity, even in high injection rate areas, similar to deep sedimentary basins in Michigan and western North Dakota. These differences in seismicity may be due to variations in permeability structure: within intrusive rocks, fluids can become narrowly focused in fractures and faults, causing an increase in local pore fluid pressure, whereas more distributed pore space in sedimentary and extrusive rocks may relax pore fluid pressure.

Shah, Anjana K.↗

DOE OSTI2019

Effect of natural fractures on determining closure pressure

Abstract In Utah FORGE reservoir, eight pressure transient tests (microhydraulic fracturing and DFIT™) show natural fracture/pressure-dependent leakoff. This behavior may lead to misinterpretation of the closure pressure (proxy for minimum principal stress). The closure pressures obtained from DFIT™ (or microhydraulic fracturing) testing may reflect shear failure along natural fractures or discontinuities rather than tensile failure and lead to inaccurate estimates of the minimum principal stress. In pressure tests conducted at the Utah FORGE site showed that reactivation or opening of natural fractures intersecting the wellbore and were suggested by multiple closure events in DFIT™ tests. In addition, comparisons between the pre- and post-well-tests FMI logs show that there are populations of induced fracture and reactivated natural fractures which are mostly vertical and sub-vertical. In this study, DFIT™ (or microhydraulic fracturing) test analysis and numerical simulations were used to suggest that hydraulic shearing of critically stressed natural fractures can contribute to multiple closure signatures and possibly lead to incorrect determination of the minimum principal stress. As other authors have previously advocated, the DFIT™ test analyses and numerical simulations suggest that better estimations of minimum principal stress may be derived by injecting at relatively high rate and pressure and insuring that tensile breakdown is reached. In addition, extended shut-in period is required to determine accurate reservoir characteristics and fluid flow regime.

Nadimi, S.↗

DOE OSTI2018

Oxidative stress induced necroptosis activation is involved in the pathogenesis of hyperoxic acute lung injury

Highlights: • Necroptosis is involved in the pathogenesis of hyperoxic acute lung injury (HALI). • Hyperoxia induced oxidative stress may induce necroposis in the lung of HALI rats. • Edaravone may inhibit oxidative stress to reduce necroposis and protect the lung against HALI. Necroptosis has been found to be involved in the pathogenesis of some lung diseases, but its role in hyperoxic acute lung injury (HALI) is still unclear. This study aimed to investigate contribution of necroptosis to the pathogenesis of HALI induced by hyperbaric hyperoxia exposure in a rat model. Rats were divided into control group, HALI group, Nec-1 (necroptosis inhibitor) group and edaravone group. Rats were exposed to pure oxygen at 250 kPa for 6 h to induce HALI. At 30 min before hyperoxia exposure, rats were intraperitoneally injected with Nec-1 or edaravone, and sacrificed at 24 h after hyperoxia exposure. Lung injury was evaluated by histology, lung water to dry ratio (W/D) and bronchoalveolar lavage fluid (BALF) biochemistry; the serum and plasma oxidative stress, expression of RIP1, RIP3 and MLKL, and interaction between RIP1 and RIP3 were determined. Results showed hyperoxia exposure significantly caused damage to lung and increased necroptotic cells and the expression of RIP1, RIP3 and MLKL. Edaravone pre-treatment not only inhibited the oxidative stress in HALI, but also reduced necroptotic cells, decreased the expression of RIP1, RIP3 and MLKL and improved lung pathology. Nec-1 pretreatment inhibited necroptosis and improved lung pathology, but had little influence on oxidative stress. This study suggests hyperoxia exposure induces oxidative stress may activate necroptosis, involving in the pathology of HALI, and strategies targeting necroptosis may become promising treatments for HALI.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle, Alabama: Field testing of MBM

A modular borehole monitoring concept has been implemented to provide a suite of well-based monitoring tools that can be deployed cost effectively in a flexible and robust package. The initial modular borehole monitoring system was deployed as part of a CO 2 injection test operated by the Southeast Regional Carbon Sequestration Partnership near Citronelle, Alabama. The Citronelle modular monitoring system transmits electrical power and signals, fibre-optic light pulses, and fluids between the surface and a reservoir. Additionally, a separate multi-conductor tubing-encapsulated line was used for borehole geophones, including a specialized clamp for casing clamping with tubing deployment. The deployment of geophones and fibre-optic cables allowed comparison testing of distributed acoustic sensing. We designed a large source effort (>64 sweeps per source point) to test fibre-optic vertical seismic profile and acquired data in 2013. The native measurement in the specific distributed acoustic sensing unit used (an iDAS from Silixa Ltd) is described as a localized strain rate. Following a processing flow of adaptive noise reduction and rebalancing the signal to dimensionless strain, improvement from repeated stacking of the source was observed. Conversion of the rebalanced strain signal to equivalent velocity units, via a scaling by local apparent velocity, allows quantitative comparison of distributed acoustic sensing and geophone data in units of velocity. We see a very good match of uncorrelated time series in both amplitude and phase, demonstrating that velocity-converted distributed acoustic sensing data can be analyzed equivalent to vertical geophones. We show that distributed acoustic sensing data, when averaged over an interval comparable to typical geophone spacing, can obtain signal-to-noise ratios of 18 dB to 24 dB below clamped geophones, a result that is variable with noise spectral amplitude because the noise characteristics are not identical. With vertical seismic profile processing, we demonstrate the effectiveness of downgoing deconvolution from the large spatial sampling of distributed acoustic sensing data, along with improved upgoing reflection quality. finally, we conclude that the extra source effort currently needed for tubing-deployed distributed acoustic sensing vertical seismic profile, as part of a modular monitoring system, is well compensated by the extra spatial sampling and lower deployment cost as compared with conventional borehole geophones.

58 GEOSCIENCES↗

DOE OSTI2021

Dilute Combustion Control Using Spiking Neural Networks

Dilute combustion with exhaust gas recirculation (EGR) in spark-ignition engines presents a cost-effective method for achieving higher levels of engine efficiency. At high levels of EGR, however, cycle-to-cycle variability (CCV) of the combustion process is exacerbated by sporadic occurrences of misfires and partial burns. Previous studies have shown that temporal deterministic patterns emerge at such conditions and certain combustion cycles have a significant influence over future events. Due to the complexity of the combustion process and the nature of CCV, harnessing all the deterministic information for control purposes has remained challenging even with physics based 0-D, 1-D, and high-fidelity computational fluid dynamics (CFD) models. In this study, we present a data-driven approach to optimize the combustion process by controlling CCV adjusting the cycle-to-cycle fuel injection quantity. Readily available data from in-cylinder pressure was used to train a spiking neural network (SNN) which learns the optimal way to manage fuel injection in order to reduce CCV while maintaining acceptable levels of fuel consumption. SNNs are particularly well suited for powertrain control applications due to their ability to be deployed on FPGA-based neuromorphic hardware which are small, inexpensive, and have a low power demand. The high-performance computing (HPC) resources of Oak Ridge National Laboratory were used to run an evolutionary-based training approach for choosing the best SNN configuration that minimizes the size of the network while achieving the desired goal. The neuromorphic hardware with the optimized SNN deployed was connected to the rapid prototyping engine control system for real-time control implementation and tested on a single cylinder version of a GM LNF 4-cylinder engine. The results show a significant reduction of CCV with a small percentage of additional fuel used to stabilize the charge.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2019

Single-Fuel Reactivity Controlled Compression Ignition Combustion Enabled by Onboard Fuel Reformation (Final Report)

Reactivity-Controlled Compression Ignition (RCCI) is a promising new advanced, low temperature combustion mode that offers efficiencies equal to or greater than conventional diesel combustion, but significantly lower soot and NOx emissions. Additionally, RCCI has enhanced controllability over the combustion process compared to some of the other advanced combustion modes. Therefore, RCCI combustion has the potential to be of great benefit to the public by increasing the efficiency and reducing the harmful pollutant emissions from transportation. However, RCCI combustion achieves this favorable level of performance through the use of two distinct fuels (e.g., diesel and gasoline, diesel and natural gas, or diesel and ethanol). In a vehicle application, RCCI therefore requires two different fuel tanks and fuel injection systems, which presents a significant drawback to the commercialization RCCI. This project proposed a new concept to enable RCCI combustion from a single parent fuel stored onboard the vehicle in a single fuel tank. Specifically, the parent fuel would be stored in the fuel tank and one branch of the fuel system would deliver the parent fuel, unaltered, to the cylinder, while a second branch of the fuel system would route the parent fuel through a fuel reformer to chemically react the fuel into a gaseous fuel mixture, called “reformate”, that has different properties from the original parent fuel. In this way, there can be two fuels at the engine that have different properties, but only one parent fuel needs to be stored and refilled in the vehicle. This initial, exploratory investigation into the feasibility of this promising new concept included reacting three potential parent fuels (diesel, gasoline, and natural gas) in a fuel reformer and determining the composition of the gaseous reformate fuel that exited the fuel reformer. There are different types of reformation processes, but since this project was a proof-of-concept investigation into the potential of Single-Fuel RCCI and the combustion process in the cylinder, a catalytic partial oxidation (CPOX) reformation process was selected, which is more straightforward to perform and control, and would potentially be easier to implement in a vehicle application compared to other reformation processes. Once the three parent fuels were reformed to varying degrees and their reformate compositions were determined, their autoignition tendencies were experimentally determined by comparing advanced combustion results of the reformate mixtures with known reference fuels. This showed that diesel fuel and its reformate fuel mixture exhibited the largest reactivity separation, which is necessary to enable RCCI combustion in the cylinder. Using diesel fuel as the parent fuel, experimental testing was conducted to provide a better understanding of the operational considerations, benefits, and drawbacks of Single-Fuel RCCI with onboard fuel reformation. The experimental testing was accompanied by computational fluid dynamics and system-level modeling to help complement the experiments and further expand the understanding. The main conclusions are that Single-Fuel RCCI using diesel fuel as the direct injected, high reactivity fuel and diesel’s gaseous reformate mixture as the premixed, low reactivity fuel is viable. The efficiencies were equal to, or sometimes greater than, RCCI with conventional fuel pairs, and the emissions were comparable as well. A key consideration is that gaseous reformate fuel contains a large amount of diluents in the fuel which displace air entering the cylinder and act as EGR inside the cylinder. Additionally, the reactivity separation between diesel and its reformate fuel mixture is larger than diesel-gasoline, but less than diesel-natural gas, and the combustion process therefore responded in a manner that was between the two conventional fuel pairs. Finally, although achievable in the cylinder, there are practical drawbacks or other disadvantages of Single-Fuel RCCI with onboard fuel reformation. First, this project used a CPOX reformation process which is exothermic, meaning that a portion of the fuel’s energy was released in the reformer where it could not be converted into useful work. This constituted an efficiency penalty of the system of on the order of 10% or 4 percentage points. Future research should therefore focus on alternate reforming processes including steam-methane reforming or autothermal reforming where there might not be any energy released in the reformer, or there could even be an energy gain in the reformer (through an endothermic process with waste heat recovery). Practical drawbacks that exist, but were beyond the scope of this preliminary proof-of-concept research project, are the added cost of the fuel reformer, which is a catalyst made of precious metals, transient response of the system, and the lifetime of the fuel reformer.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2018

Effective Elastic and Neutron Capture Cross Section Calculations Corresponding to Simulated Fluid Properties from CO2 Push-Pull Simulations

The submission contains a .xls files consisting of 10 excel sheets, which contain combined list of pressure, saturation, salinity, temperature profiles from the simulation of CO2 push-pull using Brady reservoir model and the corresponding effective compressional and shear velocity, bulk density, and fluid and time-lapse neutron capture cross section profiles of rock at times 0 day (baseline) through 14 days. First 9 sheets (each named after the corresponding CO2 push-pull simulation time) contains simulated pressure, saturation, temperature, salinity profiles and the corresponding effective elastic and neutron capture cross section profiles of rock matrix at the time of CO2 injection. Each sheet contains two sets of effective compressional velocity profiles of the rock, one based on Gassmann and the other based on Patchy saturation model. Effective neutron capture cross section calculations are done using a proprietary neutron cross-section simulator (SNUPAR) whereas for the thermodynamic properties of CO2 and bulk density of rock matrix filled with fluid, a standalone fluid substitution tool by Schlumberger is used. Last sheet in the file contains the bulk modulus of solid rock, which is inverted from the rock properties (porosity, sound speed etc) based on Gassmann model. Bulk modulus of solid rock in turn is used in the fluid substitution.

15 GEOTHERMAL ENERGY↗

DOE OSTI2018

Devising a framework to Couple Electromagnetic Geophysical Measurements with Reservoir Simulations to Monitor CO2 Plume Movement

Monitoring CO2 plume movement and extent is critical to any CO2 utilization and storage (CCUS) project but can present significant challenges. Downhole measurements (such as cross-well seismic) can give the resolution needed to track the CO2 plume. That resolution comes at a significant expense in equipment and lost production time, and the resolution is only in the vicinity of that well. Surface seismic is less expensive but has decreasing resolution with depth and may not be sensitive to changes in the CO2 plume. Conversely, electrical conductivity has been shown to be correlated with changes in phase saturation associated with CCUS operations. A novel approach for subsurface monitoring that can be applied with minimal expense and operation disruption is charged wellbore casing controlled-source electromagnetics (CWC-CSEM). This technique injects electrical current into two legacy wells. Surface measurements of the induced electric and magnetic fields are obtained with distributed data loggers. Time-lapse measurements are sensitive to changes in CO2 plume movement, but inversion requires the application of prior knowledge from prior characterization efforts. We propose a methodology where simulations are used to inform CWC-CSEM inversion while simultaneously being validated. As such, uncertainty is reduced, and more robust flow models can be devised. A three-year study is underway at the Bell Creek Field in southeast Montana (DOE-NETL Project DE-FE0028320) to evaluate the effectiveness of CWC-CSEM on an active CO2-EOR operation. Determining the relationship between electrical resistivity and phase saturations, and coupling that to a reservoir model is the thrust of this project. Initially, the relative permeability/capillary pressure relationship and the compositional oil model were identified as model parameters with high uncertainty. To elucidate the impact of these uncertain parameters, simulations were run for ten years with a single five-spot pattern. Final water saturation distributions were converted into resistivity using Archie’s Law and available fluid sampling data to provide an exploration of the time lapse electromagnetic anomaly which can be anticipated using CWC-CSEM methods. This technique has promise as an inexpensive and accurate way to monitor CO2 injection in a CCUS site.

Moodie, Nathan↗

DOE OSTI2015

A computational study of an atomizing liquid sheet

Linear instability predictions of liquid sheets injected into a gas medium are well established in the literature. These analyses are often used in Lagrangian-Eulerian spray simulations, a prominent simulation method, to model the dynamics occurring in the near-nozzle region. In the present work, these instability predictions are re-examined by first generalizing the treatment of interfacial conditions and related assumptions with a two-phase Orr-Sommerfeld (OS) system, and second, by employing highly resolved-Volume-of-Fluid (VoF) simulations. After presenting some validation exercises for both the VoF and OS solvers, the OS predictions are compared to earlier studies from the literature leading to reasonable agreement in the limit as the boundary layer thickness tends to zero. Results from VoF simulations of liquid sheet injection are used to characterize the range of scales of the liquid structures immediately before atomization. The mean value in this range is found to be approximately two to three orders of magnitude larger than the corresponding predictions from previous studies. A two-phase mixing layer under the same physical conditions is used to examine this disparity, revealing that within the linear regime, relatively good agreement exists between the VoF and OS predicted instability mechanisms. However, the most unstable mode in the linear regime is too small to cause a fracture or atomization of the liquid sheet and hence cannot be directly responsible for the atomization. The generation of a much larger mode, which emerges well beyond the linear regime, is the one causing breakup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

Unraveling the Geometry Dependence of In-Nozzle Cavitation in High-Pressure Injectors

Cavitation is an intricate multiphase phenomenon that interplays with turbulence in fluid flows. It exhibits clear duality in characteristics, being both destructive and beneficial in our daily lives and industrial processes. Despite the multitude of occurrences of this phenomenon, highly dynamic and multiphase cavitating flows have not been fundamentally well understood in guiding the effort to harness the transient and localized power generated by this process. In a microscale, multiphase flow liquid injection system, we synergistically combined experiments using time-resolved x-radiography and a novel simulation method to reveal the relationship between the injector geometry and the in-nozzle cavitation quantitatively. We demonstrate that a slight alteration of the geometry on the micrometer scale can induce distinct laminar-like or cavitating flows, validating the multiphase computational fluid dynamics simulation. Furthermore, the simulation identifies a critical geometric parameter with which the high-speed flow undergoes an intriguing transition from non-cavitating to cavitating.

FLUID DYNAMICS↗

DOE OSTI2014

Coupled THM Modeling of Hydroshearing Stimulation in Tight Fractured Volcanic Rock

Here, we use the TOUGH-FLAC simulator for coupled thermo–hydro-mechanical modeling of well stimulation for an Enhanced Geothermal System (EGS) project. We also analyze the potential for injection-induced fracturing and reactivation of natural fractures in a porous medium with associated permeability enhancement. Our analysis aims to understand how far the EGS reservoir may grow and how the hydroshearing process relates to system conditions. We analyze the enhanced reservoir, or hydrosheared zone, by studying the extent of the failure zone using an elasto-plastic model, and accounting for permeability changes as a function of the induced stresses. For both fully saturated and unsaturated medium cases, the results demonstrate how EGS reservoir growth depends on the initial fluid phase, and how the reservoir extent changes as a function of two critical parameters: (1) the coefficient of friction, and (2) the permeability-enhancement factor. Furthermore, while well stimulation is driven by pressure exceeding the hydroshearing threshold, the modeling also demonstrates how injection-induced cooling further extends the effects of stimulation.

15 GEOTHERMAL ENERGY↗