Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracturing”

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 451 records · Page 25

Quantitative Assessment of Parent Well Effect on Hydraulic Fracture Propagation at HFTS2: Insights from Cross-Well Strain Measurements and Microseismic Data

Understanding fracture propagation behavior is essential for optimizing hydraulic fracturing in unconventional reservoirs. This study demonstrates the value of integrating Low-Frequency Distributed Acoustic Sensing (LF-DAS) and microseismic data, which together provide a more complete picture of fracture growth. Using data from Hydraulic Fracturing Test Site 2 (HFTS2), we identify stress changes in depletion zones induced by parent wells as a key factor influencing fracture propagation. This result is shown by new measurements of in-situ fracture propagation velocity and fracture-hit volume (fluid volume at fracture hit?) from LF-DAS and event density from microseismic. These findings highlight the importance of considering parent well effects, well spacing, and stimulation sequencing in completion design to improve reservoir development and production efficiency.

depletion zones↗

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗

Experimental and numerical investigation of proppant embedment and conductivity reduction within a fracture in the Caney Shale, Southern Oklahoma, USA

The current worldwide energy supply is insufficient to meet the rising demand. As a result, the energy prices are expected to keep soaring despite the recent increases in a variety of renewable energy resources. Although not renewable, shale oil and gas — “unconventional” hydrocarbon resources are relatively clean forms of energy resources, which still hold a vast share of the energy market. For many oil and gas companies, meeting profitable production goals from shale reservoirs is sometimes challenging, due to the loss of fracture conductivity and premature declines in the production. In this paper we investigate the stress-dependent changes in the hydraulic conductivity of proppant-filled fractures and mechanical fracture–proppant interactions in Caney Shale, a calcareous, organic-rich mudrock, through laboratory experiments and numerical modeling. American Petroleum Institute (API) fracture conductivity tests were conducted using 2% KCl on five locations within the Caney Shale that consisted of selecting three brittle (reservoir) zones and two ductile zones. Confining pressures ranged from 1,000 psi to 12,000 psi at 210 °F. Conductivity, permeability as well as embedment were measured during the test. Also, an additional, laboratory in-situ visualization test was conducted to examine the detailed proppant-shale matrix interaction under elevated stress (3,920 psi effective stress) and temperature (252 °F), with a synthetic reservoir fluid. Our experimental results have confirmed that improved fracture conductivity is attributed to proppant size, and that the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress.

04 OIL SHALES AND TAR SANDS↗

Assessment of CO 2 Storage Potential in Naturally Fractured Reservoirs With Dual–Porosity Models

Naturally Fractured Reservoirs (NFR’s) have received little attention as potential CO 2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO 2 tends to be nonwetting in target formations and capillary forces will keep CO 2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO 2 plume. Numerical simulations are a powerful tool to understand the physics behind brine-CO 2 flow in NFR’s. Dual-porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual-porosity models are based on crude approximations of the matrix-fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO 2 transfer to the matrix. This work presents an assessment of CO 2 storage potential in NFR’s using dual-porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO 2 in the matrix blocks and the maximum CO 2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual-porosity simulations. Furthermore, we then analyze how injection rates should be limited in order to avoid early spill of CO 2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO 2 .

54 ENVIRONMENTAL SCIENCES↗

Nonsteroidal Anti-Inflammatory Drug Prescriptions Are Associated With Increased Stress Fracture Diagnosis in the US Army Population

Stress fractures are common in military personnel and endurance athletes, and nonsteroidal anti-inflammatory drug (NSAID) use is widespread in these populations. NSAIDs inhibit prostaglandin synthesis, which blunts the anabolic response of bone to physical activity and could therefore increase risk of stress fracture. The objective of this study was to determine whether prescribed NSAIDs were associated with stress fracture diagnoses among US Army soldiers. We also aimed to establish whether acetaminophen, an analgesic alternative to NSAIDs, was associated with stress fracture risk. A nested case-control study was conducted using data from the Total Army Injury and Health Outcomes Database from 2002 to 2011 (n = 1,260,168). We identified soldiers with a diagnosis of stress fracture (n = 24,146) and selected 4 controls per case matched on length of military service (n = 96,584). We identified NSAID and acetaminophen prescriptions 180 to 30 days before injury (or match date). We also identified soldiers who participated in basic combat training (BCT), a 10-week period of heightened physical activity at the onset of Army service. Among these individuals, we identified 9088 cases and 36,878 matched controls. Conditional logistic regression was used to calculate incident rate ratios (RR) for stress fracture with adjustment for sex. NSAID prescription was associated with a 2.9-fold increase (RR = 2.9, 95% confidence interval [CI] 2.8–2.9) and acetaminophen prescription with a 2.1-fold increase (RR = 2.1, 95% CI 2.0–2.2) in stress fracture risk within the total Army population. The risk was more than 5-fold greater in soldiers prescribed NSAIDs (RR = 5.3, 95% CI 4.9–5.7) and more than 4-fold greater in soldiers prescribed acetaminophen (RR= 4..4, 95% CI 3.9–4.9) during BCT. Our results reveal an association between NSAID and acetaminophen prescriptions and stress fracture risk, particularly during periods of heightened physical activity. Prospective observational studies and randomized controlled trials are needed to support these findings before clinical recommendations can be made.

60 APPLIED LIFE SCIENCES↗

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗

Numerical and Experimental Studies of Particle Settling in Real Fracture Geometries

In this study, proppant is a vital component of hydraulic stimulation operations, improving conductivity by maintaining fracture aperture. While correct placement is a necessary part of ensuring that proppant performs efficiently, the transport behavior of proppant in natural rock fractures is poorly understood. In particular, as companies pursue new propping strategies involving new types of proppant, more accurate models of proppant behavior are needed to help guide their deployment. A major difficulty with simulating reservoir-scale proppant behavior is that continuum models traditionally used to represent large-scale slurry behavior loose applicability in fracture geometries. Particle transport models are often based on representative volumes that are at the same scale or larger than fractures found in hydraulic fracturing operations, making them inappropriate for modeling these types of flows. In the absence of a first-principles approach, empirical closure relations are needed. However, even such empirical closure relationships are difficult to derive without an accurate understanding of proppant behavior on the particle level. Thus, there is a need for experiments and simulations capable of probing phenomena at the sub-fracture scale. In this paper, we present results from experimental and numerical studies investigating proppant behavior at the sub-fracture level, in particular, the role of particle dispersion during proppant settling. In the experimental study, three-dimensional printing techniques are used to accurately reproduce the topology of a fractured Marcellus shale sample inside a particle-flow cell.

15 GEOTHERMAL ENERGY↗

Effect of patterned inclusions on the fracture behavior of ceramic composites

Developments in additive manufacturing have enabled the reinforcement of ceramic composites across the scales. Understanding the fracture behavior of such composites is critical for structural design and engineering. However, systematic investigations on the relationship between the patterning characteristics and the effective fracture properties have been scarce. In this work, we propose a mixed-mode phase field model to investigate the effects of elastic modulus and pattern characteristics of circular inclusions on the fracture behavior of a ceramic composite beam. The modeling results have shown that the stiff and compliant inclusions can either repel or attract the crack. The extent of repelling and attraction can be controlled by altering the pattern characteristics. In addition, utilizing the relationship between the pattern characteristics and resulting fracture path, fracture growth can be controlled with appropriate inclusion properties. The targeted enhancement in fracture toughness can also be predicted using an effective fracture model linked with the underlying fracture path computed by the proposed mixed-mode phase field model.

36 MATERIALS SCIENCE↗

Dating fractures using luminescence

Rock fracturing (cracking) is a universal process that drives and limits chemical degradation, sediment production and erosion, and deterioration of infrastructure. Despite extensive research gains in rock mechanics on one hand and geochronology on the other, there remains a glaring gap in our ability to understand the long term evolution of natural, in situ fractures. Here we develop a novel fracture exposure dating technique, grounded in modern advances in luminescence geochronology. We apply our new dating method to a granitic boulder from a glacial outwash terrace in California, US. We conclude that the longest, clast-splitting E-W fracture appeared shortly after the boulder's deposit, whereas the secondary N-S fracture appeared 5 ka after the deposition, approximately correlating with the Last Glacial Maximum and Younger Dryas periods of the region, respectively. However, dating of the third fracture (<< 50 µm width) which does not fully split the rock, is ambiguous due to negligible daylight penetration and poor determination of fracture width. The fracture dating method presented herein brings with it the potential to decipher relationships that are crucial for the interpretation and modeling of, for example, long-term landscape and atmospheric evolution relating rock weathering to climate change and erosion.

58 GEOSCIENCES↗

Distributed Fiber Optic Sensing for in-well hydraulic fracture monitoring

This study presents the results from in-well hydraulic fracture monitoring within a horizontal well in an unconventional reservoir utilizing Distributed Fiber Optic Sensing (DFOS). An in-house-developed Brillouin-based Distributed Strain Sensing (DSS) interrogator was deployed to obtain strain measurements, complemented by a commercial Raman-based Distributed Temperature Sensing (DTS) interrogator for temperature measurements and a commercial Rayleigh-based Low-Frequency Distributed Acoustic Sensing (LF-DAS) interrogator for strain-rate measurements. Examined over a ten-day period, the spatio-temporal distribution of temperature-compensated strain obtained from DSS and DTS revealed distinct signatures of the multi-stage hydraulic fracturing process. These signatures were analyzed with respect to fracture width growth and closure, residual strain effects, and fracture conductivity near the wellbore. Fracture widths within the fracture zone were estimated for individual stages. The findings were assessed with LF-DAS measurements for further evaluation. This work integrates DFOS-measured strain, temperature, and strain-rate data for monitoring in-well hydraulic fracturing, with the goal of supporting future studies in interpreting DFOS measurements for improved understanding of hydraulic fracturing in unconventional reservoirs.

58 GEOSCIENCES↗

CO 2 mineral trapping in fractured basalt

Fractures in basalt can provide substantial surface area for reactions, and limited mass transfer in fractures can allow accumulation of cations to form carbonate minerals in geologic carbon sequestration. In this study, flood basalt and serpentinized basalt with engineered fractures were reacted in water equilibrated with 10 MPa CO 2 at 100 °C or 150 °C for up to 40 weeks. Carbonation in basalt fractures was observed as early as 6 weeks, with Mg- and Ca-bearing siderite formed in both basalts reacted at 100 °C and Mg-Fe-Ca carbonate minerals formed in the flood basalt reacted at 150 °C. X-ray μCT segmentation revealed that precipitates filled 5.4% and 15% (by volume) of the flood basalt fracture after 40 weeks of reaction at 100 °C and 150 °C, respectively. Zones of elevated carbonate abundance did not completely seal the fracture. Limited siderite clusters (<1% volume fraction) were found in localized areas in the serpentinized basalt fracture. A 1-dimensional reactive transport model developed in CrunchTope examined how geochemical gradients drive silicate mineral dissolution and carbonate precipitation in the fracture. The model predicts that siderite will form as early as 1 day after the addition of CO 2 . In conclusion, the predicted location of maximum siderite abundance is consistent with experimental observations, and the predicted total carbonate volumes are comparable to estimates derived from CT segmentation.

54 ENVIRONMENTAL SCIENCES↗

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗

The effect of ball-milling time and annealing temperature on fracture toughness of Ni-3 wt.% SiC using small punch testing

Highlights: • The fracture energy of Ni-SiC is strongly dependent on the grain size of Ni matrix • The ball-milling time governs the fracture behaviour of Ni-3 wt.% SiC • A higher proportion of big SiC particles at grain boundary is responsible for a lower fracture toughness - Abstract: A set of dispersion strengthened Ni-3 wt.% SiC samples have been prepared via a powder metallurgy route for an application in molten salt reactor (MSR) systems. A mixture of 97 wt.% Ni and 3 wt.% SiC powder was prepared by varying the ball-milling time (8, 24, 36 and 48 h). The ball-milling process was followed by powder consolidation using spark plasma sintering (SPS) and rapid cooling to room temperature. The samples were then annealed at different temperatures (300 °C, 700 °C and 850 °C). The fracture behaviour of the Ni-3 wt.% SiC samples was investigated using the Small Punch Test (SPT). A scanning electron microscope (SEM) with energy dispersive spectrometer (EDS) and electron back-scatter diffraction (EBSD), as well as transmission electron microscopy (TEM) were employed to examine the microstructure of the tested Ni-3 wt.% SiC samples. The obtained results show that the SPT fracture energy is strongly dependent on the grain size of Ni matrix, which is affected by ball-milling times and annealing temperature. It has been further found that the segregation of large SiC particles at grain boundaries promotes intergranular fracture and thus significantly reduces the fracture energy. The results indicate that the distribution and size of SiC particles at grain boundary is mainly affected by the ball-milling time. Hence, the ball-milling time governs the fracture toughness of Ni-3 wt.% SiC samples.

36 MATERIALS SCIENCE↗

Mineral Reactions in Shale Gas Reservoirs: Barite Scale Formation from Reusing Produced Water As Hydraulic Fracturing Fluid

Hydraulic fracturing for gas production is now ubiquitous in shale plays, but relatively little is known about shale-hydraulic fracturing fluid (HFF) reactions within the reservoir. To investigate reactions during the shut-in period of hydraulic fracturing, experiments were conducted flowing different HFFs through fractured Marcellus shale cores at reservoir temperature and pressure (66 °C, 20 MPa) for one week. Results indicate HFFs with hydrochloric acid cause substantial dissolution of carbonate minerals, as expected, increasing effective fracture volume (fracture volume + nearfracture matrix porosity) by 56-65%. HFFs with reused produced water composition cause precipitation of secondary minerals, particularly barite, decreasing effective fracture volume by 1-3%. Barite precipitation occurs despite the presence of antiscalants in experiments with and without shale contact and is driven in part by addition of dissolved sulfate from the decomposition of persulfate breakers in HFF at reservoir conditions. The overall effect of mineral changes on the reservoir has yet to be quantified, but the significant amount of barite scale formed by HFFs with reused produced water composition could reduce effective fracture volume. In conclusion, further study is required to extrapolate experimental results to reservoir-scale and to explore the effect that mineral changes from HFF interaction with shale might have on gas production.

03 NATURAL GAS↗

Revisiting the Analytical Solutions of Heat Transport in Fractured Reservoirs Using a Generalized Multirate Memory Function

Numerous analytical solutions have been developed for modeling thermal perturbations to underground formations caused by deep-well injection of fluids. Each solution has been derived for a specific boundary value problem and a simplified flow network with one set of parallel fractures. In this paper, new generalized solutions G*( x, s ) are developed using (existing) global transfer functions G$^*_0(x,s)$ and a new memory function g*( s ), where x and s are the space and Laplace variable. The memory function represents the solutions of conductive heat exchange between fractures and matrix blocks and between fractured aquifers and unfractured aquitards. The memory function is developed to account for multirate exchange induced by different shapes, sizes, properties, and volumetric fractions of matrix blocks bounded by multiple sets of orthogonal fractures with different spacing. The global transfer functions represent the fundamental solutions to convective, convective-conductive, and convective-dispersive heat transport in fractures (or aquifers) without exchange and are available for various (1-D linear, 1-D radial, 2-D dipole, and single-well injection-withdrawal) flow fields. The new solutions with exchange are developed using $G^*(x,s) = B^*(s)G^*_0{x, s[1 + Θg^*(s)]}$, thereby greatly simplifying solution development in a novel way, where Θ and B*( s ) are a fracture-matrix scaling factor and the boundary condition function. As a result, the new solutions are applied to several example problems, showing that heat transport in fractured aquifers is significantly impacted by (1) thermal dispersion in fractures that is rarely considered, (2) multirate heat exchange with a wide range of size and anisotropy of rectangular matrix blocks, and (3) heat exchange between aquifers and aquitards.

58 GEOSCIENCES↗

Characterizing the Influence of Fracture Density on Network Scale Transport

The topology of natural fracture networks is inherently linked to the structure of the fluid velocity field and transport therein. Here we study the impact of network density on flow and transport behaviors. We stochastically generate fracture networks of varying density and simulate flow and transport with a discrete fracture network model, which fully resolves network topology at the fracture scale. We study conservative solute trajectories with Lagrangian particle tracking and find that as fracture density decreases, solute channelization to large local fractures increases, thereby reducing plume spreading. Furthermore, in sparse networks mean particle travel distance increases and local network features, such as velocity zones where flow is counter to the primary pressure gradient, become increasingly important for transport. As the network density increases, network statistics homogenize and such local features have a reduced impact. We quantify local topological influence on transport behavior with an effective tortuosity parameter, which measures the ratio of total advective distance to linear distance at the fracture scale; large tortuosity values are correlated to slow-velocity regions. Additionally, these large tortuosity, slow-velocity regions delay downstream transport and enhance tailing on particle breakthrough curves. Finally, we predict transport with an upscaled, Bernoulli spatial Markov random walk model and parameterize local topological influences with a novel tortuosity parameter. Furthermore, Bernoulli model predictions improve when sampling from a tortuosity distribution, as opposed to a fixed value as has previously been done, suggesting that local network topological features must be carefully considered in upscaled modeling efforts of fracture network systems.

58 GEOSCIENCES↗

Fatigue as the missing link between bone fragility and fracture

The prevention of fragility fractures in bone-pathologic fractures resulting from daily activity and mostly occurring in the elderly population-has been a long-Term clinical quest. Recent research indicating that falls in the elderly might be the consequence of fracture rather than its cause has raised fundamental questions about the origin of fragility fractures. Is day-To-day cyclic loading, instead of a single-load event such as a fall, the main cause of progressively growing fractures? Are fragility fractures predominantly affected by bone quality rather than bone mass, which is the clinical indicator of fracture risk? Do osteocytes actively participate in the bone repair process? In the Perspective, we discuss the central role of cyclic fatigue in bone fragility fracture.

59 BASIC BIOLOGICAL SCIENCES↗

3D Experimental Measurement of Lattice Strain and Fracture Behavior of Sand Particles Using Synchrotron X-Ray Diffraction and Tomography

3D synchrotron X-ray diffraction (3DXRD) and synchrotron micro-computed tomography (SMT) techniques were used to measure and monitor the lattice strain evolution and fracture behavior of natural Ottawa sand particles subjected to 1D compression loading. The particle-averaged lattice strain within sand particles was measured using 3DXRD and then was used to calculate the corresponding lattice stress tensor. In addition, the evolution and mode of fracture of sand particles was investigated using high-resolution 3D SMT images. The results of diffraction data analyses revealed that the major principal component of the lattice strain or stress tensor increased in most of the particles as the global applied compressive load increased until the onset of fracture. Particle fracture and subsequent rearrangements caused significant variation and fluctuations in measured lattice strain/stress values from one particle to another and from one load step to the next one. SMT image analysis at the particle-scale showed that cracks in fractured sand particles generally initiate and propagate along the plane that connects the two contact points. Fractured particles initially split into two or three major fragments followed by disintegration into multiple smaller fragments in some cases. In conclusion, microscale analysis of fractured particles showed that particle position, morphology, the number and location of contact points play a major role in the occurrence of particle fracture in confined comminution of the sand assembly.

47 OTHER INSTRUMENTATION↗