Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing fluid”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2022

Subsurface stress criticality associated with fluid injection and determined using machine learning

Machine-learning methods and apparatus are disclosed to determine critical state or other parameters related to fluid-driven failure of a terrestrial locale impacted by anthropogenic activities such as hydraulic fracturing, hydrocarbon extraction, wastewater disposal, or geothermal harvesting. Acoustic emission, seismic waves, or other detectable indicators of microscopic processes are sensed. A classifier is trained using time series of microscopic data along with corresponding data of critical state or failure events. In disclosed examples, random forests and artificial neural networks are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to live data from a fluid injection locale in order to assess a frictional state, assess seismic hazard, assess permeability, make predictions regarding a future fluid-driven failure event, or drive engineering solutions for mitigation or remediation. Variations are disclosed.

Johnson, Paul Allan↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Distributed Acoustic Sensing as a Distributed Hydraulic Sensor in Fractured Bedrock

Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three-dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO 2 sequestration, and remediation of groundwater in fractured bedrock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Small-angle Neutron Scattering (SANS) Characterization of Clay- and Carbonate-rich Shale at Elevated Pressures

Unconventional oil and gas from shale formations have emerged as some of the fastest growing energy resources in the United States, providing both cleaner energy to consumers and reducing the nation’s reliance on energy imports. To properly harness these important natural resources, the nanopore structure of associated shales must be fully understood, particularly under hydraulic fracturing conditions, where they are exposed to both overburden compressive and hydrostatic fluid pressures. The current study uses small-angle neutron scattering (SANS) to characterize pore structure, including porosity, pore accessibility, and pore size distribution, in the 1–100 nm regime at elevated pressures for mineralogically distinct clay- and carbonate-rich shales from the Permian Basin. Unlike typical porosity measurement techniques, SANS is uniquely capable of characterizing both open and closed porosity, allowing measurement of how pore accessibility changes with pressure and determination of the size range of accessible versus inaccessible pores. The porosity of the clay-rich shale was 7.7%, compared to 0.51% for the carbonate-rich shale. However, only 2.6% of the nanopores in the carbonate-rich shale were inaccessible to water at 8 kPSI (55.1 MPa) compared to 7.8% for the clay-rich shale. Furthermore, our analyses indicated that the closed pores fall within distinct size ranges, likely corresponding with the chemical nature of the pore host material. These results provide valuable insight into the effects of shale petrophysical properties on hydrocarbon extraction from unconventional reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗
DOE OSTI2021

Comparative geochemistry of flowback chemistry from the Utica/Point Pleasant and Marcellus formations

Flowback/Produced fluid samples were collected from several wells from two Utica/Point Pleasant (UPP) sites (UPPW and UPPS) in Ohio, and one Marcellus (Marcellus Shale Energy and Environment Laboratory (MSEEL)) site in West Virginia over a period of approximately two years. Although these formations have different ages, depositional environments, diagenetic histories, and geochemical and mineralogical compositions (i.e. the UPP is significantly more carbonate rich than the Marcellus which is more siliceous), analysis of trends in fluid species over time shows that, overall, the TDS and major solubilized elements (Na, Ca, Cl) in the UPP and Marcellus brines are remarkably similar. Total dissolved solutes (TDS) in these brines ranged from approximately 40 to 250 g/L salt, and in general, concentrations increased with time elapsed since natural gas well completion and stimulation. The behavior of Na, Br, and Cl suggests that the produced water signatures from these formations are largely derived from the native formational brines which display evidence of originating from evaporated seawater. There is a strong correlation between Cl and Br, indicating that both species behave conservatively, and the similarity among each of these brines suggests no appreciable contribution of salt from halite dissolution because Br is excluded from the halite structure. Cl/Br ratios in the brines range from ~80 to 120 (mg/L/mg/L). Other elements, such as K, which readily reacts between fluids and ion exchange sites on clays, generally exhibit conservative behavior for an individual site, but show significant variations among each of the different well pads. The concentrations of Sr and Ba vary dramatically among well sites, and increase with respect to Cl- over time, suggesting increasing solubilization, presumably from desorption from clay minerals or dissolution of carbonates or sulfates from the source formation(s). The UPPW well site has very low Ba due to high-sulfate input fluid, which resulted in precipitation of barite/celestite in the brines. In contrast the UPPS well site had elevated Sr (~ 3500 mg/L), presumably due to the use of Sr-rich recycled brine used in hydraulic fracturing. The Marcellus site had the highest Ba concentrations (up to 10 g/L) and highest Ba/Sr ratios in the fluids, due to the high concentration of barium in the Marcellus target (~ 1000 ppm, as compared to ~200 ppm in the UPP). These observations suggest that solutes in the FP fluids are derived from native brines, water-rock interactions that have occurred over geologic time scales, as well as some contribution from contemporaneous reactions in the subsurface. The results also show that the composition of the injected fluid can influence flowback fluid chemistry and possibly production efficiency.

58 GEOSCIENCES↗
DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗
DOE OSTI2022

Modeling gas breakthrough and flow phenomena through engineered barrier systems using a discrete fracture approach

Compacted clays are being considered to build engineered barrier systems (EBS) intended for the safe isolation of high-level nuclear waste (HLW) and spent nuclear fuel (SNF). The corrosion of the metallic canister containing the HLW/SNF will lead to the generation and buildup of the gas pressure in the more internal part of the clay buffer. This phenomenon would eventually trigger the formation and propagation of fractures in the clay barrier, jeopardizing its safety functions. Here, in this work, we propose to use the fragmentation technique (MFT) to model evolving fractures in clays triggered by gas pressurization. The MFT has been successfully used to model the formation of fractures in concrete, drying cracks in soil, hydraulic and thermo-fractures in rocks. In this work, we extend the MFT to deal with multiphase fluid flow in deformable porous media, and we upgraded a fully coupled computer finite element code using the extended technique. The proposed approach is first verified against analytical solutions and is then applied to model gas breakthrough experiments in clays. A very satisfactory performance of the method is observed in all the analyses, showing the potential of the MFT to tackle multiphase flow problems in deformable porous media with evolving discontinuities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Relocating microseismicity from downhole monitoring of the Decatur CCS site using a modified double-difference algorithm

SUMMARY The injection of CO2 at the Decatur carbon capture and storage site has generated significant microseismic activity, which occurs in distinct spatial clusters up to approximately 2.2 km from the primary injection well. Accurate and precise event locations are vital for the characterization of the microseismicity to help understand the reservoir response to the CO2 injection, whilst enabling the identification of minor faults and fractures below the resolution of conventional active seismic imaging. However, microseismic monitoring of fluid injection sites, such as Decatur, is often performed using a network of borehole sensors often from a single well. While these downhole sensors have excellent detection capabilities, their poor azimuthal coverage limits the ability to precisely determine event locations. We have developed a modified double-difference relocation algorithm suitable for both 1-D and 3-D velocity models, and which incorporates differential back azimuth observations to allow the benefits of the original double-difference algorithm to be applicable to a downhole microseismic monitoring setting. Applying the modified double-difference algorithm to the microseismicity at Decatur, we have successfully relocated 4293 events. The relocation included over 59 million observations for 757 285 event pairs, split across seven geographic regions. Despite the majority of observations being recorded in only two boreholes, with an almost identical azimuthal coverage, the results have shown to be reliable with significantly reduced residuals and low uncertainties associated with the final locations. We have analysed the residuals in terms of their association with each geographic region, data type, station and individual events, to fully appreciate their influence in the inversion and the fit of the data to the final set of event locations. For each region, the relocated seismicity has become less diffuse with improved clustering, and with newly visible linear features often orientated in a NE–SW direction. These results show the potential improvements that can be made to microseismic event locations recorded by a borehole network with a limited and variable azimuthal distribution.

58 GEOSCIENCES↗
DOE OSTI2021

Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA

The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.

58 GEOSCIENCES↗
DOE OSTI2025

Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites

Presentation slides on “Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. Accurate reservoir characterization is fundamental for the safe and efficient operation of carbon capture, utilization, and storage (CCUS) projects. This process is not only critical during initial site selection but also plays a key role in continuously updating reservoir geomodels throughout injection operations to effectively manage subsurface storage sites. Our primary objective is to develop a comprehensive understanding of fracture networks within CO2 injection reservoirs, which are critical for predicting the behavior of injected fluids, ensuring storage integrity post-injection, and mitigating risks such as leakage or induced seismicity.

artifical intelligence / machine learning (AI/ML)↗
DOE OSTI2022

Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (Final Report)

This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.

02 PETROLEUM↗
DOE OSTI2020

Multi-physics modeling of injected nanoparticles effect on remediation of CO 2 leakage through cracks

One potential risk in CO 2 sequestration is the leakage of carbon dioxide, which can result in contamination of underground water, creating potential threats to existing ecosystems. The common leakage pathway is through the pre-existing fractures or discontinuities within cement in the wellbore, incurred by the environmental conditions imposed on the cement. Injecting nanoparticles into pre-existing cracks is one of the most recently proposed ideas for mitigating fracture propagation in cement CO 2 sequestration. To demonstrate the feasibility of this new technology, a numerical approach was taken in this work, as it is challenging to investigate it in a laboratory setting. We proposed a coupled ALE (Arbitrary Lagrangian–Eulerian)–DEM (Discrete Element Method)–peridynamic modeling strategy within the LS-Dyna package to investigate the intertwined interaction among the CO 2 fluid flow, native fluid (e.g., brine), particle clusters, and cracks within the cement. The numerical results demonstrate that injected nanoparticles can effectively reduce the pressure exerted on the crack surface. Accordingly, the potential fracture propagation at the crack tip would be reduced compared to corresponding cases without nanoparticels as pressurized by fluid flow. This result verifies the effectiveness of proposed nanoparticle injection technology. Finally, using this established modeling strategy, the effect of filling particles on the fracture mitigation for different crack geometries (e.g. particle cluster patterns, aspect ratio of crack aperture and length) and CO 2 reservoir pressure are examined. The result shows that injected particles successfully reduce the fracture propagation in these scenarios.

42 ENGINEERING↗
DOE OSTI2023

Autonomous Fracture Conductivity Using Expandable Proppants in Enhanced Geothermal Systems

Summary Early thermal breakthrough in enhanced geothermal systems (EGS) due to the presence of preferential flow channels is a major challenge that endangers efficient and economic heat extraction in such systems. Previous studies mainly focused on adjusting circulation rates of the working fluid, which still leaves significant amounts of untapped heat behind. Currently, there is a lack of technologies for altering flow distribution within the fracture network to achieve uniform heat sweeping in the reservoir. This work presents a novel concept for making proppants to autonomously control fracture conductivity based on the surrounding temperature. Here, proppants with negative thermal expansion coefficients have demonstrated the capability for appropriate fracture conductivity adjustment as a function of temperature to achieve uniform flow across the fracture network. Particle-particle interactions governing such functions are explicitly modeled, and then the Lattice Boltzmann methods (LBM) is used to determine the potential impact of closure stress and temperature changes on the permeability of the proposed proppant packs. Microscale analyses are further used to determine the required material properties to achieve a certain improvement in the permeability of the proppant pack. Our analyses show an enhancement in permeability and the associated fracture conductivity by half of their initial values. Field-scale analysis further confirms the effectiveness of the proposed concept as 31.4% more heat can be extracted from EGS over 50 years of production when the proposed proppants are used. Such novel proppants may effectively delay thermal breakthrough, sweep heat from larger rock volumes, and elongate the life span of the EGS project.

Engineering↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗
DOE OSTI2025

Microseismicity Modulation Due To Changes in Geothermal Production at San Emidio, Nevada, USA

Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.

Guo, Hao [University of Wisconsin-Madison, WI (Uni↗