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

Metagenome-assembled genomes provide insight into the metabolic potential during early production of Hydraulic Fracturing Test Site 2 in the Delaware Basin

Demand for natural gas continues to climb in the United States, having reached a record monthly high of 104.9 billion cubic feet per day (Bcf/d) in November 2023. Hydraulic fracturing, a technique used to extract natural gas and oil from deep underground reservoirs, involves injecting large volumes of fluid, proppant, and chemical additives into shale units. This is followed by a “shut-in” period, during which the fracture fluid remains pressurized in the well for several weeks. The microbial processes that occur within the reservoir during this shut-in period are not well understood; yet, these reactions may significantly impact the structural integrity and overall recovery of oil and gas from the well. To shed light on this critical phase, we conducted an analysis of both pre-shut-in material alongside production fluid collected throughout the initial production phase at the Hydraulic Fracturing Test Site 2 (HFTS 2) located in the prolific Wolfcamp formation within the Permian Delaware Basin of west Texas, USA. Specifically, we aimed to assess the microbial ecology and functional potential of the microbial community during this crucial time frame. Prior analysis of 16S rRNA sequencing data through the first 35 days of production revealed a strong selection for a Clostridia species corresponding to a significant decrease in microbial diversity. Here, we performed a metagenomic analysis of produced water sampled on Day 33 of production. This analysis yielded three high-quality metagenome-assembled genomes (MAGs), one of which was a Clostridia draft genome closely related to the recently classified Petromonas tenebris. This draft genome likely represents the dominant Clostridia species observed in our 16S rRNA profile. Annotation of the MAGs revealed the presence of genes involved in critical metabolic processes, including thiosulfate reduction, mixed acid fermentation, and biofilm formation. These findings suggest that this microbial community has the potential to contribute to well souring, biocorrosion, and biofouling within the reservoir. Our research provides unique insights into the early stages of production in one of the most prolific unconventional plays in the United States, with important implications for well management and energy recovery.

natural gas↗

Fluid-structure interaction solver for transient dynamics of fracturing media

Computer implemented methods include providing a numerical simulation solid mesh comprising a plurality of solid nodes and a numerical simulation fluid mesh comprising a plurality of fluid nodes, and performing a numerical simulation with the solid mesh and fluid mesh including: determining a fluid-solid interaction between the fluid mesh and the solid mesh at a time step by (i) immersing boundaries of the solid and fluid meshes, (ii) computing current corrected fluid and solid velocities using an interaction parameter applied to an intermediate fluid-solid relative velocity, and (iii) computing a fluid-solid interaction force based on the current corrected fluid and solid velocities.

Rougier, Esteban↗

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗

Opening-mode fracturing and cementation during hydrocarbon generation in shale: An example from the Barnett Shale, Delaware Basin, West Texas

Relative timing of fracturing is a key input for predictive fracture models, but timing information for fractures is commonly difficult to obtain. In this study, we used crosscutting relations and fluid inclusion assemblage temperatures from fracture cements from a few well-documented sampled fractures, combined with a one-dimensional burial history model, to establish timing for three generations of opening-mode fractures in a Barnett Shale core from the southern part of the Delaware Basin, Pecos County, West Texas. A burial history model is presented for the cored well and matched to measured vitrinite reflectance in samples from the core, and bottomhole temperature in the well. Here, the earliest fractures (group 1) likely formed due to early fluid-expulsion events (ca. 300 Ma) and were folded during host-rock compaction. Later group 2 fractures are sealed with fibrous barite containing primary, liquid hydrocarbon inclusions (mean homogenization temperature [T h ] = –9°C) and aqueous fluid inclusions (mean T h = 108.1°C). Group 2 fractures likely formed in response to fluid overpressure associated with cracking of type II kerogen to oil. Group 3 vertical fractures are up to 2 m in height with kinematic apertures ranging from less than 0.05 to 1.4 mm, partly open, and strike dominantly 010°–020°. Sequentially trapped aqueous fluid inclusions in fracture-spanning quartz cement bridges (mean T h = 110°C in crack-seal texture and 128°C in post-crack-seal fracture cement) record fracture opening under increasing temperature, inferred to reflect increasing burial, with continued overpressuring during the Triassic to Late Cretaceous. Some group 3 fractures may have continued to fill during Cenozoic uplift.

02 PETROLEUM↗

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗

Deterministic Discrete Fracture Network (DFN) Model for the EGS Collab Project on the 4850 Level of the Sanford Underground Research Facility (SURF)

The EGS Collab is conducting hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF) located in Lead, South Dakota. A total of eight ~60m-long subhorizontal boreholes were drilled from the 4850 Level (~1.5 km below the ground surface) into the crystalline rock of this former mine. Six of these holes are used for geophysical monitoring, one is used for hydraulic fracture stimulation, and the remaining hole was designed as a production borehole that receives water from the injection well via the induced and natural fracture system. The primary goal of creating the discrete fracture network model is to show that these modeling methods are critical for the development of enhanced geothermal systems (EGS). This includes the prediction of rock behavior during fracturing and during an extended period of water flow between the parallel injection and production boreholes. Understanding the results from the induced fracturing and flow is complicated by the presence of significant natural fractures that interact with the stimulation and/or flow pathways. The delineation and characterization of natural fractures is thus an important part of the project, and therefore a model of the Discrete Fracture Network (DFN) was developed on a deterministic basis. The DFN was populated using observations and interpretations integrated from drift (horizontal passageways that allow access in the underground) fracture mapping, analysis of core recovered from the eight boreholes, borehole televiewer logs and videos, and observations of flow between and within boreholes and in the drift. The natural fracture system is dominated by a pervasive northwest-trending, steeply dipping shear system that is identifiable in the drifts and the core. Hydraulic fracture stimulation, flow/tracer circulation tests, and geophysical monitoring revealed that the behavior of the injected water, and perhaps the growth of induced fractures, has been significantly influenced by the existing fractures identified in the DFN.

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

Large-Volume Stimulation of Rock for Greatly Enhanced Fluids Recovery using Targeted Seismic-Assisted Hydraulic Fracturing (Final Technical Report)

This project has developed and demonstrated a new technology for large-volume and targeted comminution of rock in low permeability formations to enhance recovery from unconventional oil and gas (UOG) resources. The technology is based on a strategically designed interaction of multiple induced seismic pulses that assist the hydraulic fracturing process to enhance shear and multi-planar crack formation. This greatly increased rock stimulation, through bulk comminution, is expected to cause significant increase in permeability leading to enhancement of recovery factors for sub- surface fluids. The proposed technology is especially applicable for enhanced recovery in emerging UOG plays, such as ductile shales that are resistant to opening-mode fracturing by conventional hydraulic fracturing processes. The project combines an integrated experimental and computation approach to develop and demonstrate a modular technology that can be easily implement in the field to augment current practices. The effort integrates a fundamental scientific understanding of dynamic material response under constraint and damage-induced permeability and porosity enhancements at multiple length scales, along with models of comminution due to the local release of kinetic energy associated with high shear strain rate of dynamic deformation. The results will be validated through small-scale experiments and then implemented in a lab-scale field test to demonstrate the developed technology. The project involves collaboration between a solid mechanician, a materials scientist, and a petroleum engineering geo-mechanician. It builds on a successful history of success and collaboration that includes both fundamental scientific exploration and technology development.

04 OIL SHALES AND TAR SANDS↗

Geochemical reactions and alteration of pore architecture in saturated shale after injection of stimulation fluid

Pore architecture regulates fluid flow between unconventional shale reservoir and hydraulically-induced fractures. Imbibition of stimulation fluid may change this architecture and alter hydrocarbon flow. Hydrothermal experiments were conducted at reservoir conditions (125 °C, 45 MPa) to test two hypotheses: 1) Shale, not stimulation fluid, dominates the geochemistry of an unconventional reservoir containing formation water; and 2) Mineral dissolution and precipitation induced by stimulation fluid is transient and manifests across micro-, meso-, and macro-scales of pore architecture. Experiments replicated a shut-in well in the Cretaceous Baxter Shale, Green River Basin, Wyoming USA. Stimulation fluid was injected into one experiment after formation water (I = 0.35 mol/kg, pH = 6.4) and core reacted for 48 days. This novel approach equilibrated rock and formation water and saturated pores with formation water before introducing stimulation fluid. The second experiment served as a control. Trends of aqueous calcium, silica and aluminum in the injection experiment suggest transient dissolution of calcite and feldspar and/or quartz as well as clay and barite precipitation; mineralogic evidence was limited to calcite dissolution and barite precipitation. The rock maintained reducing conditions (Eh = +0.08 to -0.16 V) despite injection of oxidizing stimulation fluid (Eh = 1.1 V). pH of the stimulation fluid-formation water mixture evolved from 2.3 to pre-injection values (~6) within 24 h. The results indicate that mineral dissolution and precipitation manifest in macropores with barely detectable alteration to micro- and mesopores. Formation water or organic matter in pores may have inhibited access of stimulation fluid to micro- and mesopores

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗

Abstract for CRADA between National Energy Technology Laboratory and Shell International Exploration & Production, Inc

Introducing CO₂ into geothermal systems as a working fluid in reservoirs can enhance geothermal conductivity, production, and pressure maintenance. A cross-disciplinary interaction of geothermal reservoir stimulation and CO₂ utilization satisfies renewable energy demands and operations that support sustainable energy infrastructure. Challenges to implementing CO₂-stimulated geothermal enhancement (CS-GE) include (1) accurately characterizing and imaging CO₂-stimulated geothermal reservoirs; (2) quantitatively inferring CS-GE evolution under current and future engineered conditions for cost effective operations; and (3) monitoring resources by improving observational methods to advance the understanding of complex geothermal systems for sweep efficiency and, ultimately, cost effectiveness. NETL and Shell will collaborate under this CRADA to develop software to image key features, including CO₂-stimulated fracture imaging, CO₂ fluid sweep imaging and heat transfer and exchange imaging by leveraging available datasets and applying advanced Artificial Intelligence/Machine Learning (AI/ML), multi-level data analytics and data/information fusion to better understand the comprehensive mechanisms of CO₂-stimulated geothermal systems.

15 GEOTHERMAL ENERGY↗

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗

The evolution of seafloor venting from hydrate-sealed gas reservoirs

Here we use 3D seismic data to show that three rows of seafloor gas mounds can be traced downward to leak points that lie at the hydrate-gas contact within three individual dipping coarse-grained sand bodies in the Terrebonne Basin, Walker Ridge block 313, northern Gulf of Mexico. We predict the overpressure within the sand bodies by assuming that the gas pressure at the vent points equal the least principal stress. We interpret that free gas accumulates at the base of the hydrate stability zone, trapped by the overlying sand which has a high methane hydrate saturation. The free gas accumulates until the gas pressure at the base of the hydrate stability zone reaches the least principal stress in the overlying mudrock, whereupon hydraulic fractures form and fluids are vented to the surface. The warm rising fluids and perhaps localized exothermic formation of hydrate raise the local salinity and temperature. This process progressively shifts the base of the hydrate stability zone to shallower depths and dissociates the hydrate seal within the sand, which creates new leak points and results in the observed migration of the seafloor vents. Within the southwest corner of the Terrebonne Basin, this process has repeated multiple times within the Blue, Orange, and Green sands. Here this study shows how the hydrate stability zone can influence the location of fluid expulsion and in turn be affected by the warm, rising, saline fluids sealed by hydrate.

58 GEOSCIENCES↗

Temperature-responsive smart tracers for field-measurement of inter-well thermal evolution: Heterogeneous kinetics and field demonstration

Temperature-responsive smart tracers enable advanced warning of a “premature thermal breakthrough” and may therefore improve reservoir management and reduce financial risk. A successful calculation of inter-well fluid temperatures requires that temperature-dependent kinetics are known, which typically result from homogeneous batch reactor experiments. However, recent meso-scale field experiments at the Altona Field Laboratory involving fluid flow in a discretely-fractured reservoir suggest that silica-fluid interfaces may accelerate hydrolysis kinetics. Here, the Arrhenius parameters of phenyl acetate hydrolysis are investigated under heterogeneous reaction conditions in packed-bed column experiments. The breakthrough curve of the reaction product is compared to an inert reference tracer (carbon-cored nanoparticles). Temperatures experienced during field testing ranging from 10 to 40 °C were studied using a phosphate buffer solution to maintain a near-neutral pH, as found in the reservoir. The empirically-determined pre-exponential factor and activation energy values are subsequently used in a first-order kinetic model to improve calculations of effective reservoir temperatures for these meso-scale field tests. Here, the results suggest that hydrolysis reactions can directly monitor inter-well reservoir temperatures, but only if the influence of solid-fluid interactions are carefully accounted for.

15 GEOTHERMAL ENERGY↗