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 361 records · Page 20

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗

Evolution of Pore Structure and Permeability of Rocks Under Hydrothermal Conditions (Final Report)

The physical and transport properties of porous rocks can be altered by a variety of diagenetic, metamorphic, and tectonic processes, and the changes that result are of critical importance to such industrial applications as resource recovery, carbon dioxide sequestration, and waste isolation in geologic formations. The associated interrelationships between rocks, pore fluids and deformation are also key to understanding many natural processes, including dynamic metamorphism, fault mechanics, fault stability and pressure solution creep. In this project, we investigated the changes of permeability and pore geometry owing to inelastic deformation by solution-transfer, brittle fracturing and dislocation creep. In particular, we studied the coupling between pore fluids and deformation in fluid-filled quartz, calcite, mudstones and ultramafic rocks and examine the effects of coupled mechanical and chemical processes on the evolution of porosity and permeability under hydrothermal conditions. The investigations used a combination of techniques, that included triaxial laboratory experiments; uniaxial compressive loading with in situ observations of microstructure; permeability under hydrostatic compression; observations of microstructure using optical and electron microscopes, and micro-CT imaging; and numerical calculations. Laboratory experiments were designed to provide mechanical and transport data under conditions that isolate the particular mechanisms responsible for the changes. The data obtained were used to quantify changes in surface roughness, porosity, pore dimensions and their spatial fluctuations. The results of the experiments and data from image analyses were compared to the results of network, finite-difference and other numerical models to verify the validity of experimentally established relations between permeability and other rock properties. A bibliography for the period from 2015-2019 is appended below. New results obtained in the final year and the unfunded extension period are detailed below.

58 GEOSCIENCES↗

Utah FORGE 1-2551: Multi-Stage Fracturing System and Well Tractor to Enable Zonal Isolation During Stimulation and EGS Operations in Horizontal Wellbores - Workshop Presentation and Report

Included here are a presentation recording, slides, and report on the Multi-Stage Fracturing System and Well Tractor to Enable Zonal Isolation During Stimulation and EGS Operations in Horizontal Wellbores project by Colorado School of Mines, presented by Dr. William W. Fleckenstein. One goal of the project is to develop and demonstrate the use of low-cost, rapid multistage fracture stimulation technology with cemented casing frac sleeves that eliminates packers used in both conventional stimulation and injection flow control. The second goal is to develop and demonstrate the devices to effectively detect and control the flow of heat-carrier fluid solely from a long-reach injector through the network of induced and existing fractures and produced from an open-hole long-reach well to improve heat recovery. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

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

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

58 GEOSCIENCES↗

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗

Experimental flow-through a single fracture with monolayer proppant at reservoir conditions: A case study on Caney Shale, Southwest Oklahoma, USA

Global energy systems are undergoing a crucial transition to clean energy sources, thereby moving away from fossil-fuel based energy with the ultimate aim of reducing carbon dioxide emissions. It is essential to continue to advance in exploring and developing shale gas resources as this will underpin the global supply of natural gas that is needed to support the energy transition. The challenge hereby is that the slow fluid flow of the shale matrix due to its very low permeability requires large-scale volume fracturing to produce a conductive fracture network that enhances the connectivity between the shale formation and the wellbore. This work examines the influence of a thin proppant layer on a single fracture using a flow-through experiment, and also explores whether proppant embedment in the Caney shale is affected by the rock mineralogy, surface roughness, fluids, confining stress, time, temperature and bedding. The experiment was conducted for a duration of nine days (216 h), experimental temperature was varied from ambient temperature to reservoir temperature of 125 °C (257 °F) and confining stress was varied from 367 psi (2.53MPa) to a maximum of 4011.82 psi (27.66 MPa). Here, we find that the conductivity of the fracture is primarily influenced by the layer of proppant used, surface roughness, mineralogy, fluids, temperature and closure stress.

Caney shale↗

Experimental study on rheological and settling properties of shape memory polymer for fracture sealing in geothermal formations

This article studies the rheology and annular flow of a smart lost circulation material (LCM) at 160 °C, using a high-temperature flow loop. The smart LCM is a shape memory polymer (SMP) that activates at high temperatures, and its particle size increases to seal fractures in geothermal formations. SMP was mixed with the base fluid in two different concentrations, 1.0 and 3.0 wt%, to study rheology, wellbore hydraulics, activation process, and settling behavior under different testing conditions. The results of this study showed that the SMP could be activated at high temperatures. An increase of 80–100% in the particle size was observed at 160 °C. The mud samples showed a high shear-thinning behavior at the two concentrations with a Power-law flow index (n) ranging between 0.025 and 0.101. No additional frictional pressure losses were observed when SMP was added to the base fluid. SMP particles showed an excellent suspension at 1.0 wt% while, at 3.0 wt%, a bed was formed at a low flow rate and without pipe rotation. Increasing the drill pipe rotational speed or flow rate effectively removed the bed and homogeneously dispersed the SMP particles, ensuring a better sealing efficiency. SMP particle dispersion in inclined wells was better than in horizontal wells. Moreover, the findings of this study help optimize the lost circulation treatment by considering a wide range of operating parameters that can further be extended to different systems and geometries.

02 PETROLEUM↗

Utah FORGE: Fault Reactivation Through Fluid Injection Induced Seismicity Laboratory Experiments

Included are results from shear reactivation experiments on laboratory faults pre-loaded close to failure and reactivated by the injection of fluid into the fault. The sample comprises a single-inclined-fracture (SIF) transecting a cylindrical sample of Westerly granite. All experiments are conducted at ambient temperature and follow a similar protocol: (i) application of confining stresses (3MPa) on the fault fully saturated with DI water, (ii) shear-mobilization through the increase of axial loading at a constant displacement rate until a post-peak steady-state condition is reached, (iii) reduction of axial loading and related shear stress to a prescribed fraction of the peak steady-state frictional strength (typically 60% to 90%, representing intermediate to high magnitudes) and (iv), fault reactivation triggered by a stepwise increase of pore pressure on the fault in 0.1 MPa increments held constant for 1-5 minutes. Mechanical data from three ISCO pumps connected to a Temco pressure vessel measure axial, confining, and fault-related parameters, including fluid pressure (kPa), fluid flow rate (mL/min), and axial displacement (mm). See included code for initial data analysis and visualization for select experiments. Resource names represent experiment numbers found in the "Read Me" file, which describes each experimental setup and parameters.

15 GEOTHERMAL ENERGY↗

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

Fracture Adjacent Matrix Permeability: Insights from a Direct Experimental Approach

Description of a measurement technique, using a triple-pressure tap Hassler Style Core holder, that facilitates pressure measurement perpendicular to a fracture. Through the measurement of permeability changes at multiple distances from the fracture, this system enables an understanding of how end effects and localized ‘skin effects’ can impact the ability of fluids to move to and from the bulk matrix. Measurements of this type may lead to improved input values for reservoir simulations involving discrete fractures and/or dual-permeability.

47 OTHER INSTRUMENTATION↗

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

Fracture, Dissolution, and Cementation Events in Ordovician Carbonate Reservoirs, Tarim Basin, NW China

Ordovician carbonate rocks of the Yijianfang Formation in the Tabei Uplift, Tarim Basin, contain deeply buried (>6000 m), highly productive oil and gas reservoirs associated with large cavities (>10 m). Previous workers inferred that large cavities are paleocaves (paleokarst) formed near the surface and subsequently buried. Alternately, caves may have formed by dissolution at depth along faults. Using 227 samples from 16 cores, we document textures and cement compositions bearing on cavity histories with petrographic, high-resolution scanning electron microscopy (SEM), isotopic, and fluid inclusion microthermometric observations. Results show that dissolution occurred at depth and was caused by (1) acidic fluids derived from Middle-Late Silurian and/or Devonian-Permian hydrocarbon generation and maturation, (2) high-temperature fluids, of which some were associated with Late Permian igneous activity, and (3) Mg-rich fluids that accompanied Jurassic-Cretaceous deformation and the formation of partially open fractures and stylobreccias (fault breccias). The relative paragenetic sequence of the structure-related diagenesis suggests seven stages of fracturing, dissolution, and cementation. Mottle fabrics in the Yijianfang Formation contain argillaceous carbonate-rich silt and are bioturbation features formed within the marine environment. Those mottled fabrics differ from clearly karstic features in the overlying Lianglitage Formation, which formed by near-surface dissolution and subsequent infilling of cavities by allochthonous sediment. Mottle fabrics are crosscut by compacted fractures filled with phreatic-vadose marine cements and followed by subsequent generations of cement-filled fractures and vugs indicating that some fractures and vugs became cement filled prior to later dissolution events. Calcite cements in fractures and vugs show progressively depleted values of δ18O documenting cement precipitation within the shallow (~220 m), intermediate (~625 m), and deep (~2000 m) diagenetic environments. Deep (mesogenetic) dissolution associated with fractures is therefore the principal source of the high porosity-permeability in the reservoir, consistent with other pieces of evidence for cavities localized near faults.

58 GEOSCIENCES↗

Minimizing exposure to legacy wells and avoiding conflict between storage projects: Exploring area of review as a screening tool

Elevated pressure from large-volume injection is a key driver of risk and project cost. If transmissive features (e.g., non-isolating wells or fracture systems) are present, increased injection-zone pressure can drive fluids from depth toward protected freshwater resources. In US Carbon Capture and Storage (CCS) law, the area at risk is known as the Area of Review (AoR). The size and number of potentially transmissive features to be evaluated and possibly remediated or managed is a function of the size and location of the AoR. The size of the AoR depends on several variables, including properties of the injection zone, properties of protected resources, and injection rate and duration. Evaluation of the intersection of these variables across a portfolio of sites highlights the injection zone depth and boundary conditions as top-level controls. Deep injection, use of multiple stacked injection zones, reduced injection rate and choice of injection well location can all be used to minimize AoR and the number of potentially transmissive features within it. Here, we introduce the concept of pressure space (defined as connected pore volume times pressure) as the key subsurface commodity for CO 2 storage and we suggest that it forms a more robust basis for leasing and regulation than pore space alone.

58 GEOSCIENCES↗

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗