Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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

The Thermal Response of a Packed Bed Thermal Energy Storage System upon Saturated Steam Injection Using Distributed Temperature Sensing

The effectiveness of a thermal energy storage (TES) system is typically characterized with the help of thermal stratification or temperature gradients along the direction of heat injection, which is typically the flow direction of heat transfer fluid. The steepness of temperature gradients are a direct indicator of the effectiveness or efficiency of the heat storage or dispatch process. The temperature gradient evolution along the packed bed of ceramic particles upon saturated steam injection is presented in this work. Distributed temperature sensing based on optical frequency domain reflectometry was deployed in a packed bed of ceramic particles to capture the thermal front evolution in the axial direction. The physical processes accompanying steam injection in packed beds are complex due to phase change, transitioning two-phase flow, and changes in condensate accumulation. Therefore, the variation of thermal response of the TES system for various steam injection flow rates was experimentally studied using a high-resolution distributed temperature sensing system in a chemically inert alumina particle-packed bed. Distinct zones of different heat transfer modes were observed during the steam injection experiments. A distinct conduction zone, evident from diffuse thermal fronts, was observed at low flow rates, and these thermal gradients became sharper as the flow rate increased. The diffuse thermal fronts in the heat storage media suggest a low exergy efficiency of the TES system, as energy losses started initiating before a significant fraction of the bed was saturated with steam.

25 ENERGY STORAGE↗

The effectiveness of D 2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign

Abstract Edge localized modes (ELMs) in H-mode plasmas erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes D 2 pellet injection on the DIII-D fusion experiment used as an ELM mitigation technique applied during the 2016 tungsten Metal Rings Campaign to reduce W erosion during ELMs. The 400.9 nm photon wavelength line emission intensity of tungsten atoms (WI) filterscope channels and Langmuir probes were used to infer the gross erosion rate of tungsten-coated tiles installed in the divertor of DIII-D. D 2 mass injection rates ranging from 34 to 41 arbitrary units (A.U.) and no D 2 injection resulted in a similar total W erosion rate during ELMs (intra-ELM). On average, results show a 29% increase in the total gross W erosion rate with intermediate mass injection rates (∼13–23 A.U.) compared to the no pellets and the highest injection rate cases. On average, the fast D 2 mass injection rate cases had 15% less erosion in the inter-ELM phase than the case with no pellets. Generally, higher D 2 mass injection rates increased the ELM frequency, and the highest injection rates reduced the average erosion per ELM and fractional carbon impurities at the top of the pedestal by nearly 40% when compared to the no-pellet case. As expected, a higher D 2 pellet injection rate led to a higher plasma density and lower plasma temperature in the divertor. Additionally, an increasing divertor inter-ELM plasma electron density directly correlated to more frequent pellet injection and a decrease in both the average gross intra-ELM W erosion and the total gross intra-ELM W erosion rate. Simulations of intra-ELM erosion using the ‘free-streaming plus recycling model’ (FSRM) underestimate W erosion during pellet injection by about 30% on average. The discrepancies between the experimental measurements and the FSRM intra-ELM W erosion predictions are postulated to be due to C/W material mixing. A simple analytic mixed-material model is presented and results in better agreement with the experimental data. These results highlight the importance of incorporating the effects of a mixed-material layer in the analysis of PFC erosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a High-Pressure Fuel System for Use with Dimethyl Ether

The paper documents the modeling and experimental work on a common rail fuel injection system for Dimethyl Ether, a potential diesel substitute with a low carbon intensity signature. The DME fuel system is deployed on a light duty 2.2L compression ignition engine. The paper describes the injector optimization to shift to higher flows to account for the lower heating value and density of the DME when compared to diesel. The type of the injection system used for the DME application is an advanced rendering of the Common rail noted for a one-piece piston-needle injector construction and a solenoid driven spill valve featuring a pressure balanced poppet. A dedicated high-pressure fuel pump designed to pressurize DME is used. The design results in a fast acting open and close injection event, reduced leakage, with reduced cavitation in the fuel injector volume. Design parameters for system optimization included fill and spill orifices, needle lift, bias spring, and injector hole size. The design model provides good correlation of the instantaneous rates of injection with experiments across a wide range of pressure and injection timings. Proposed performance milestones for the design included similar DME injection duration to the diesel counterpart for same fuel energy injected into the cylinder to retain high engine cycle efficiency. Here, the dedicated DME design provided reduced hydraulic delays of 50%. Tests demonstrated sustained operation at pressures of 1000 bar, with capability to reach 1500bar. Durability tests showed no cavitation-deterioration over a 200-hour test cycle by means of spray imaging and hardware inspection.

De Ojeda, William [WM International Engineering, L↗

Evaluating the Aftershock Duration of Induced Earthquakes

As the carbon sequestration community prepares to scale up the number and size of commercial operations, the need for tools and methods to assess and mitigate risks associated with these operations becomes increasingly important. One outstanding question is whether aftershocks of induced events decay quickly after injection operations cease or if aftershock activity persists for hundreds of years before returning to background levels more akin to tectonic events (Stein and Liu, 2009). Appropriate estimates of the aftershock duration impact several operational management decisions including mitigation strategies and post-injection monitoring for seismic activity. It is hypothesized that induced earthquake rates may diminish more quickly after injection is stopped, owing to higher stressing rates from injected fluids. Alternatively, it is plausible that only the first event in the sequence is induced by increased fluid overpressures, whereas subsequent events (e.g., aftershocks) respond to the stored tectonic stresses and static and dynamic stress changes due to the mainshock (Keranen et al., 2013). Further, if the aftershock duration can be linked to stressing rates due to injection, then it follows that operational strategies to reduce seismic hazard by reducing injection rates or volumes may be successful. However, if aftershocks of induced events are relieving stored tectonic stresses, then altering injection volumes may not alleviate ongoing seismic activity. Furthermore, knowledge of an aftershock duration could aid in the determination of an appropriate post-injection monitoring period for ongoing seismicity, which is a factor in overall operational costs. In this study, we model induced seismicity sequences in Oklahoma with a coupled Coulomb rate–state earthquake rate model (Dieterich, 1994; Kroll et al., 2017) to estimate aftershocks durations. Results for the current study indicate that elevated rates of aftershock activity following induced mainshocks return to background seismicity rates in less than five years, contrary to the tens to hundreds of years observed for tectonic aftershocks.

58 GEOSCIENCES↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗

Experimental Study of Transport Behavior of Swellable Microgel Particles in Superpermeable Channels for Conformance Control

Summary Gel treatment is an effective way to attack excessive water production during oil development. The transport behavior of gel materials in reservoirs is of crucial importance to the effectiveness of gel treatments. The aim of this paper is investigating the transport behavior of swellable micrometer-sized preformed particle gels (PPGs, or microgels) through superpermeable (super-K) channels. Sandpacks with permeabilities ranging from 27 to 221 darcies were used to mimic the super-K channels. Multiple pressure sensors were applied along the sandpack models to monitor the propagation behavior of the microgels. The tested microgel particles could transport through the super-K channels, and a higher driving pressure gradient was required when the particle/pore size ratio was larger. The pressure gradient distribution along the super-K channels was relatively uniform when the particle/pore ratio was low (less than 1.3). However, the inlet section would show increasingly higher pressure gradients as the particle/pore ratio was increased, indicating increased difficulty in propagation. The propagation of the gel particles was significantly slower compared with the carrying fluid. The delayed propagation behavior was more pronounced when the particle/pore ratio was higher. The injection pressure was much less sensitive to the injection flow rate compared with a Newtonian fluid. The gel dispersion exhibited an apparent shear thinning (pseudoplastic) behavior when transporting through the porous channels. Breakage of the gel particles was observed especially at high superficial velocities. The particle breakage was partially responsible for the apparent shear thinning behavior. The breakage phenomenon was in favor of deep placement of the gel particles. The channel permeabilities were significantly reduced by the microgels, bringing sufficient resistance to subsequent waterflooding (more than 99.5%). At given matching size conditions, softer gels were more likely to establish in-depth placement and uniform water blocking capacity in the channels. The microgel particles exhibited salinity-responsive behavior to the post-brine flush. The gel particles could shrink and reswell according to the salinity of the injected water. Possibilities were discussed to use this salinity-responsive behavior. Also, the microgels exhibited a particular disproportionate permeability reduction (DPR) effect. After gel injection, the channel permeability to water flow was reduced by more than 20 to 92 times of the permeability to oil flow. This work provides important support to understand the transport behavior of gel particles in super-K channels. The achievements are helpful for gel product selection and gel treatment design.

Engineering↗

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗

Developing a roadmap for carbon capture, and storage in Oklahoma by assessing the viability of stacked storage

Abstract The Intergovernmental Panel on Climate Change concludes that CO 2 capture and storage (CCS) is critical for climate‐stabilizing energy transitions. In CCS, captured CO 2 is sequestered in saline aquifers within sedimentary basins. The CO 2 storage capacity and the rate of injection are functions of the geology of the saline aquifer, which is uncertain. To minimize impacts of this uncertainty, CCS projects could include backup plans, such as co‐locating geologic CO 2 storage (GCS) sites with or near existing CO 2 ‐enhanced oil recovery (CO 2 ‐EOR) operations. These “stacked storage” projects could hedge against uncertainty in the saline formation performance because captured CO 2 could be injected into either location in the event of unexpected events (e.g., the injectivity decreases). Here, we investigate the possibility and ramifications of developing CCS networks in Oklahoma that are amendable to stacked storage. We find that stacked storage is possible in Oklahoma but the counties with the lowest‐cost saline storage resources do not have existing CO 2 ‐EOR operations. At the systems level, we find it is slightly more expensive (e.g., $1/tCO 2 to $5/tCO 2 ) to site GCS in counties with CO 2 ‐EOR projects. This increased expense is largely due to increased CO 2 transportation costs because hundreds of km of additional pipeline is required to capture CO 2 from the lowest‐cost sources. Overall, our results suggest that it is optimal to build more pipelines and avoid injecting CO 2 in some of the lowest‐cost saline storage resources, to enable capturing CO 2 from the least‐cost sources. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

42 ENGINEERING↗

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗

The impact of chromium coatings on Zircaloy cladding deformation behavior under reactivity-initiated accident-like mechanical loading conditions

A reactivity-initiated accident (RIA) occurs when a control rod ejection or control blade drop causes an increase in the fission rate. The injection of energy results in an increase in fuel temperature which in turn causes rapid thermal expansion of the fuel pellet. This thermal expansion may result in pellet-cladding mechanical interaction (PCMI) in which the fuel imparts a mechanical strain to the cladding. PCMI may cause the cladding to fail, and thus, the mechanical response of cladding due to PCMI must be investigated when characterizing new cladding materials. Chromium-coated Zircaloy-4 is a near-term accident-tolerant fuel cladding that exhibits improved high-temperature oxidation resistance. Modified burst testing was utilized to experimentally simulate the effects of PCMI on both uncoated and chromium-coated Zircaloy cladding samples at hot zero power conditions. Samples were coated using either cold spraying or physical vapor deposition to understand the differences in behavior that the coating application method may cause. Digital image correlation was used to analyze images of the deforming specimens to extract the in-situ strain behavior of the cladding. The uncoated specimens burst at hoop strains ranging from 8.8 % to 17.2 %. The cold-spray chromium-coated Zircaloy specimens burst at hoop strains of 7.0 % to 11.0 %. The physical vapor deposition coated tubes burst at hoop strains of 9.1 % to 11.5 %. These results indicate that the chromium coating causes a loss in the ductility of the cladding. Furthermore, the higher burst hoop strains of the physical vapor deposition-coated samples relative to the cold-spray samples indicate that the cold-spraying technique causes more of a loss in ductility than physical vapor deposition. All samples burst at higher hoop strains than those expected to occur in an RIA for fresh cladding.

Accident-tolerant fuel↗

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using the Shallow Strain Tensor to Characterize Deep Geologic Reservoirs

Abstract Storing and recovering water, carbon, and heat from geologic reservoirs is central to managing resources in a changing climate. We tested the hypothesis that the strain tensor caused by injecting or producing fluids can be measured at shallow depths and interpreted to advance understanding of underlying deep aquifers or reservoirs. Geodetic‐grade strainmeters were deployed at 30 m depth overlying the Bartlesville Formation, a 500‐m‐deep sandstone near Tulsa, OK. The strainmeters are 220 m east of injection well 9A completed in a permeable lens at the base of the Bartlesville Formation. Water was injected into well 9A at approximately 1.0 L/s during four tests that ranged in duration from a few hours to a few weeks. The horizontal strain increased (tension) and the circumferential strain was a few times larger than the radial strain. The vertical strain decreased (compression) during injection. Strain rates were approximately 100 nε/day during the first few hours, but the rates decreased and were approximately 10 nε/day during most of the tests. Four independent methods of poroelastic simulation and inversion predict reservoir properties and geometries that are similar to each other and consistent with independent information about the reservoir. All strain interpretations predict that a boundary to the permeable lens occurs beneath the vicinity of the strainmeters, which is consistent with core data from the site. The boundary of the permeable lens is located by matching the vertical, radial and circumferential strains, which demonstrates the value of measuring the strain tensor.

Murdoch, Lawrence C.↗

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗

Analytical Functions for 200 West Pump-and-Treat SCADA Sensor Data

Historical operations at the U.S. Department of Energy’s Hanford Site included disposal of waste fluids to the subsurface in the 200 West Area on the Hanford Central Plateau. Subsequent infiltration of fluids has resulted in groundwater contamination with carbon tetrachloride, nitrate, uranium, technetium-99, and other contaminants. A pump-and-treat (P&T) system, with an extraction/injection well network and an aboveground treatment plant, was implemented as part of interim and final remedies in the 200 West area. The HYPATIA single-page web application (part of the SOCRATES suite) is being developed to provide access to and analysis of chemistry and treatment facility sensor data for this 200 West P&T system. For the web application, analytical algorithms were developed to perform summing, differencing, smoothing, outlier detection, change-point detection, mass flow rate, and injectivity calculations on the data. Candidate algorithms were identified and tested, with the best-performing algorithms then assembled for implementation in HYPATIA. Because HYPATIA is hosted on the Amazon Web Services (AWS) cloud computing platform, algorithms were implemented in a back-end AWS Lambda function that can be called by the HYPATIA front end. The Lambda function applies the requested data processing to specified data via functions written in R, Python, and JavaScript. Development, testing, and review of the data analysis algorithms was completed under an NQA-1 quality program. This new HYPATIA functionality will provide information to support site decisions regarding P&T system performance and optimization.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

CO 2 -Responsive Fracturing Fluids for Enhanced Geothermal Systems (Final Report)

Our group has recently developed StimuFrac, a non-toxic stimuli-responsive fracturing fluid consisting of a CO 2 -reactive polymer which has shown at the lab-scale to consistently fracture rock cores at significantly lower net pressures in a range of representative geothermal pressure/temperature conditions. However, until now the mechanism/s responsible for more effective fracturing, of critical importance to optimize fracturing performance as well as strategize injection methodologies for field deployment, was not understood. In this document, we report (1) on the two main mechanisms responsible for fracturing rock at lower net pressures with StimuFrac; (2) the phase behavior of StimuFrac/CO 2 under geothermal wellbore conditions; and (3) based on high-temperature true triaxial stimulations, detailed evidence that StimuFrac/CO 2 is the best performing stimulation fluid under EGS T/P conditions as compared to water, waterless CO 2 , and CO 2 /water fracturing fluids. This is because i) it requires significantly lower volumes of CO 2 due to its reduced leak off into the formation by the crosslinked polymer solution; ii) large fractures can be generated reproducibly at both low and high CO 2 injection flow rates, and iii) the reversible (previously reported) viscosity increase of StimuFrac could be beneficial to transport proppants when they become available for EGS. These results were particularly evident for hot nearly dry rock as well as partially and nearly fully water-saturated granitic rock. GTO requires StimuFrac to be evaluated in fully water-saturated rock to determine whether the above-described performance applies under these conditions. Since (1) GTO considers there is no enough evidence for Sa different StimuFrac formulation to work below 90C (where full water-saturation in an open system is possible) and (2) no polyaxial loading frame larger than a few centimeters that can do hydraulic fracturing tests while maintaining the rock sample fully saturated with water at 200 C exists; PNNL concludes that the only way to determine if StimuFrac represents an advanced fracturing fluid alternative for EGS, is to perform a stimulation in an actual EGS reservoir.

15 GEOTHERMAL ENERGY↗