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 577 records · Page 32

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗

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.↗

Projected Seismic Activity at the Tiger Stripe Fractures on Enceladus, Saturn, From an Analog Study of Tidally Modulated Icequakes Within the Ross Ice Shelf, Antarctica

Abstract Dissipation of tidal energy is expected to generate seismicity on icy‐ocean worlds; however, the distribution and timing of this seismic activity throughout an orbital cycle is not known. We used new observations from an icy‐ocean‐world analog environment on Earth to examine the relationship between tidally driven tensile stress and seismic activity within an ice shell. We investigated a pair of rifts within Antarctica's Ross Ice Shelf which are tidally stressed in a manner analogous to the orbital cycle of tidal stress experienced by Enceladus' Tiger Stripe Fractures. We found that seismic activity at the Antarctic rifts is sensitive to both the amplitude and the rate of tensile stress across the rifts. We combined these findings with calculated stress values along Enceladus' Tiger Stripe Fractures to predict seismic‐activity levels expected along the ice‐shell fractures. We predict a peak in seismicity along the four Tiger Stripe Fractures when Enceladus is 90°–120° past pericenter in its orbit around Saturn, at which point tensile stresses would reach ∼2/3 of their maximum value. We also used the magnitude distribution of icequakes along Antarctic rifts to investigate implications for the likely size of stick‐slip rupture patches along icy faults on Enceladus. Our findings predict that the Tiger Stripe Fractures should produce sustained, low‐magnitude seismic events that involve rupture along discrete portions of each fracture's total length. We predict that seismicity would fall to 50% of peak levels when stresses across the Tiger Stripe Fractures are dominantly compressional.

Olsen, Kira G.↗

Aperture and roughness govern iron oxide passivation in olivine fractures during carbon mineralization

Efficient carbon dioxide mineralization in silicate-rich formations depends on reaction rates within fracture networks, yet the role of fracture microstructure remains poorly constrained. We investigated how fracture aperture and surface roughness influence mineralization by reacting synthetic forsteritic olivine fractures with controlled apertures (0.5–1.6 millimeters) and roughness (6–16 micrometers) with supercritical carbon dioxide at 13.7 megapascals and 185 degrees Celsius. Optical microscopy and Raman spectroscopy showed that rougher fracture surfaces promoted iron oxide precipitation, which passivated olivine surfaces, limiting magnesium availability for carbonate formation while removing iron as a reactant. These effects were more pronounced in smaller-aperture fractures. A geochemical model incorporating surface topography and passivation reproduced the observed trends, confirming that increased roughness can reduce reactivity despite greater surface area. Our results demonstrate that fracture microstructure strongly influences iron oxide passivation and carbon dioxide mineralization, with important implications for predicting long-term performance of engineered carbon dioxide storage systems.

54 ENVIRONMENTAL SCIENCES↗

Effect of hook characteristics on the fracture behaviour of dissimilar friction stir welded aluminium alloy and mild steel sheets

Under tensile shear loading, fracture modes of dissimilar lap welds produced by friction stir scribe technology were studied. Three fracture modes were observed. For zone A fracture, the initial crack was restrained, and the joint ultimately fractured in the base mild steel. For zone B fracture, the initial crack progressed through the aluminium sheet just above the Al/steel interface. For zone C fracture, the initial crack proceeded along the steel hook to the aluminium sheet surface. Fracture mode and joint strength were greatly influenced by steel hook size, and the steel hook size was affected by welding parameters and tool scribe height. In this study, the experimental joint strength achieved the calculated joint load limit.

Wang, Tianhao↗

Quantified fracture (joint) clustering in Archean basement, Wyoming: application of the normalized correlation count method

We demonstrate statistically significant self-organized clustering over a length scale range from 10 -2 to 10 1 m for north-striking opening-mode fractures (joints) in Late Archean Mount Owen Quartz Monzonite. Spatial arrangement is a critical fracture network attribute that until recently has only been assessed qualitatively. We use normalized fracture intensity plots and the normalized correlation count (NCC) method of Marrett et al. to discriminate clustered from randomly placed or evenly spaced patterns quantitatively over a wide range of length scales and to test the statistical significance of the resulting patterns. We propose a procedure for interpreting cluster patterns on NCC diagrams generated by the freely available spatial analysis software CorrCount. Results illustrate the efficacy of NCC to measure fracture clustering patterns in texturally homogeneous Archean granitic rock in a setting distant (>2 km) from folds or faults. In their current geological setting, these regional fractures are conduits for water flow and their patterns – and the NCC approach to defining clusters – may be useful guides to the spatial arrangement style and clustering magnitude of conductive fractures in other, less accessible fractured basement rocks.

02 PETROLEUM↗

Numerical Modeling of CO 2 and Brine Leakage through Open Fracture in a Fault Zone: Open Channel Flow or Darcy Flow

Understanding fluids migration and leakage risk along the fault zone is necessary to guarantee the safety of CO 2 geological storage. The validity of Darcy’s law gets challenged in dealing with the flow in open fractures since the occurring of turbulence flow. In this study, we develop a 2D model with usage of T2Well, an integrated wellbore-reservoir simulator, to investigate the leakage problem along open fractures which are embedded in a fault zone from the deep injection reservoir to shallow aquifers. The results record a positive feedback of gas expansion and pressure response in fracture, which causes a quick downward propagation of highly gas saturated zone from the top of fracture and an easy gas breakthrough in the shallower aquifers. The decreasing of aperture size of fracture significantly enhances the leakage rates in fracture, but with less influences as aperture increases. In comparison, the Equivalent Porous Media models show a good approximation with the momentum model of large apertures but poor for the small one. Nevertheless, the differences are small in terms of final CO 2 distribution among various aquifers, suggesting that Darcy’s law may be still “effective” in solving flow problem along fractures in a constant injection system at a large time scale.

58 GEOSCIENCES↗

Prediction of the Shear Failure of Opened Rock Fractures and Implications for Rock Slope Stability Evaluation

Rock slope commonly fails due to the shear failure of rock fractures. Shear strength of rock fractures are reduced substantially once the fracture surfaces are mismatched or opened. We propose a new criterion to predict the shear strength of rock fractures in different opening states. The degree of interlocking representing the true asperity contact area is incorporated into the modified model of Saeb and Amadei. The effect of fracture opening on asperity dilation and degradation is separately considered. The transitional stress that is a critical parameter involved in the model is analytically determined based on energy consideration. The new model is validated with experimental results from direct shear tests on synthetic fractures with regular‐shaped asperities. Good agreement between the analytical solution and the experimental data confirms the capacity of the proposed model. Therefore, the model has great potential for assessing the stability of rock slopes where fractures are often opened due to stress relief and engineering disturbances.

Li, Yingchun (ORCID:0000000201083704)↗

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↗

Surface Seismic Monitoring of Hydraulic Fracturing Test Site (HFTS) in the Midland Basin, Texas

Low-frequency seismic events of long duration were recorded by a surface array of 48, 3-component geophones during the fracture stimulation of two offset horizontal wells in the Middle Wolfcamp Formation (8000-8100 ft. below the geophone array) at the Hydraulic Fracturing Test Site (HFTS) in Reagan County, TX. Each recorded event was approximately 5-10 seconds long, with no clear first wave arrival; the seismic energy was largely visible on the vertical component suggesting that the events are dominated by P waves. A finite temporal moveout was observed across the seismic network in event arrival times indicating a causal relationship with the local source of deformation. The low-frequency seismic events observed in this study can be distinguished from the more commonly known microseismic events by their: 1) emergent waveform, 2) lack of clear P wave and S wave arrivals, 3) longer duration (few seconds), and 4) lower-frequency energy content (10-60 Hz). Previous researchers have observed similar low-frequency, long duration events in conjunction with the hydraulic fracturing of hydrocarbon reservoirs and have referred to them as long-period, long-duration (LPLD) events in the literature. LPLD events may be important to a comprehensive understanding of fracture stimulation because these events are thought to represent: 1) movement along pre-existing discontinuities such as bedding planes and fractures/faults, 2) deformation of ductile rocks with high clay content, and 3) the jerky tensile opening of hydraulic fractures. LPLD events resulting from slow shear-slip on faults and fractures of relatively larger dimension or rupture in the clay rich rock are collectively referred to as slow-slip seismicity, which was first observed in tectonic and volcanic settings, areas with high pore pressure like hydraulically fractured reservoirs. Because LPLDs are thought to represent deformation mechanisms not captured by microseismic analysis (which dominantly indicate shear failure of brittle rock), we believe that incorporating LPLD events into a new methodology for delimiting stimulated reservoir volume might improve predictions of well production.

Kumar, Abhash↗

Hydraulic Fracturing Test Site (HFTS) – Project Overview and Summary of Results

The Hydraulic Fracture Test Site (HFTS) is a field-based hydraulic fracturing research experiment performed in the West Texas Permian (Midland) Basin. The HFTS includes $25 million of hydraulic fracturing research that is centralized around eleven horizontal wells fractured with over 400 stages in the Upper and Middle Wolfcamp formations as well as two recompleted legacy wells. As part of the HFTS experiment, and in addition to the comprehensive field data that was collected, approximately 600 feet of excellent quality, whole core was obtained by drilling a high-angle core well through the induced hydraulic fractures at the test site. Based upon observations of the acquired core, the understanding of hydraulic fracture propagation and effectiveness of proppant placement is challenging current thinking. Well interference testing and in-situ reservoir pressure measurements using bottom hole gauges in producing wells and discrete multilevel pressure gauges in an offset monitor well aid in understating fracture connectivity and conductivity over time.

Ciezobka, Jordan↗

Laboratory Evaluation of a Natural Gas–Based Foamed Fracturing Fluid

Abstract Foams have been used as hydraulic fracturing fluids to reduce water usage and minimize the potentially deleterious impact on water-sensitive formations. Traditionally, carbon dioxide (CO2) and nitrogen (N2) have been used as the internal phase in these foamed fluids. Hydraulic fracturing with natural gas (NG) is a relatively inexpensive option, particularly if NG produced from the wellhead can be used without significant processing. In an ongoing program sponsored by the US Department of Energy (DOE), an alternative fracturing process is being developed that uses NG-based foam. Previously, the optimal thermodynamic pathway was identified to transform wellhead NG into pressurized NG suitable for use as the internal phase in a foamed fracturing fluid. Recent work has focused on preparing a NG-based foam at surface conditions typically encountered in hydraulic fracturing and measuring the stability and rheological properties of the foam. In addition, the transient response of the foam during fracture initiation was simulated using a fast-acting solenoid valve. A single base-fluid mixture was prepared by combining a commercially available viscosifier and foaming surfactant with water. The base fluid was then injected into a tee using a water pump. Simultaneously, liquefied natural gas (LNG) was pressurized using a cryogenic pump, vaporized using a heat exchanger, and injected into the tee to mix with the base fluid and generate foam. The foam then flowed through approximately 300 ft of 0.312-in. inside diameter (ID) tubing equipped with pressure transducers at several locations. The test fixture included a sight glass to visually inspect the quality of the foam. This paper reports on findings related to foam stability and rheology and compares these results to previous studies on foamed fracturing fluids.

Beck, Griffin↗

Application of Phase-Field Techniques to Hydraulically- and Deformation-Induced Fracture.

Phase-field techniques provide an alternative approach to fracture problems which mitigate some of the computational expense associated with tracking the crack interface and the coalescence of individual fractures. The technique is extended to apply to hydraulically driven fracture such as would occur during fracking or CO 2 sequestration. Additionally, the technique is applied to a stainless steel specimen used in the Sandia Fracture Challenge. It was found that the phase-field model performs very well, at least qualitatively, in both deformation-induced fracture and hydraulically-induced fracture, though spurious hourglassing modes were observed during coupled hydralically-induced fracture. Future work would include performing additional quantitative benchmark tests and updating the model as needed.

42 ENGINEERING↗

A Numerical Tool for the Design of Enhanced Geothermal Reservoirs based on Modeling of Thermo-Hydro-Mechanical Processes and Explicit Representation of 3D Discrete Fracture Networks

This project has developed a coupled thermo-hydro-mechanical numerical modeling tool for prediction of energy production from Enhanced Geothermal Reservoirs (EGS). The tool is unique because it can explicitly represent thousands of pre-existing natural fractures. The capabilities developed in this project were built on the already existing expertise of Itasca arising from thirty years of research and development in discrete element modeling of fractured rocks. A thermal module suitable for analysis of different heat transfer mechanisms in fractured rocks percolated by fluids was developed and coupled with the already existing hydro-mechanical calculation in 3DEC, a commercial software for analysis of jointed rock based on the discrete element method (DEM). The heat transfer module simulates heat conduction, heat advection by moving fluid in rock fractures, and convective heat exchange at the fracture walls between the rock and the moving fluid. Techniques for spatial discretization of heat transfer in EGS-type problems were evaluated, and a boundary layer meshing approach was developed. Also, given the significant differences in time scale of physical processes involved, temporal integration schemes for coupled advection-diffusion processes and coupled thermal-hydro-mechanical processes were developed by taking advantage of both implicit and explicit time integration and also optimizing the coupling intervals. These improvements were necessary for modeling EGS reservoirs. Also, the DFN capacities in 3DEC were extended to include a larger number of fractures. Research was carried out on approaches to simplify the DFN so that the model can meet computational demands, while preserving thermo-hydro-mechanical characteristics of the DFN. The developed modules were coupled with the already existing capabilities of 3DEC. The results were presented as a submitted conference paper (Riahi et al., 2018a). A series of benchmark and test examples showed the developed capabilities. Finally, the developed tool was tested for a large-scale reservoir model for longer periods of injection and production. The thermal draw-down curves were generated and presented. The project successfully met all the milestones laid out in the proposal. Results demonstrated the feasibility of developing a fast and reliable software for the EGS community to use for assessment and prediction of reservoir performance. Some of its applications include assessment of stimulated volume and produced temperature and power. It can also be used for improved engineering, design and optimization of EGS reservoirs as well as other unconventional resources such as shale gas. Finally, the 3D analysis presented in this study showed that simulating EGS reservoirs with tens of thousands of fractures is feasible, and confirmed that the developed tool can serve toward meeting the goals and objectives of Geothermal Technologies Office in overcoming technical barriers associated with commercial EGS development.

15 GEOTHERMAL ENERGY↗

Modeling Flow and Transport in Fractured Rock using Machine Learning [Slides]

Fractured systems in the subsurface play a role in many natural and engineered applications such as geologic carbon sequestration, hydraulic fracturing and underground nuclear test detection. Structural information (fracture size, orientation, etc.) plays a key role in governing the dominant physics for these systems but can only be known statistically. Traditional modeling approaches either ignore or idealize structural information at these larger scales because we lack a computational framework that utilizes it in its entirety. The work presented here integrates computational physics, machine learning and graph theory to make a paradigm shift from computationally intensive high-fidelity models to coarse-scale graphs without loss of critical structural information. We exploit the underlying discrete structure of fracture networks in systems considering flow through fractures and fracture propagation. We demonstrate that compact graph representations require significantly fewer degrees of freedom (dof) to capture micro-fracture information and further accelerate these models with Machine Learning.

58 GEOSCIENCES↗

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

Clinical effects of 3D printing-assisted posterolateral incision in the treatment of ankle fractures involving the posterior malleolus

Objective To explore the clinical outcomes of a 3D printing-assisted posterolateral approach for the treatment of ankle fractures involving the posterior malleolus. Methods A total of 51 patients with ankle fractures involving the posterior malleolus admitted to our hospital from January 2018 to December 2019 were selected. The patients were divided into 3D printing group (28 cases) and control group (23 cases). 3D printing was performed for ankle fractures, followed by printing of a solid model and simulation of the operation on the 3D model. The operation was then performed according to the preoperative plan, including open reduction and internal fixation via the posterolateral approach with the patient in the prone position. Routine x-ray and CT examinations of the ankle joint were performed, and ankle function was evaluated using the American Foot and Ankle Surgery Association (AOFAS) ankle-hindfoot score. Results All patients underwent x-ray and CT examinations. All fractures healed clinically, without loss of reduction or failure of internal fixation. Good clinical effects were achieved in both groups of patients. The operation time, intraoperative blood loss and intraoperative fluoroscopy frequency in the 3D printing group were significantly less than those in the control group ( p < 0.05). There was no significant difference between the two groups in the anatomical reduction rate of fractures or the incidence of surgical complications ( p > 0.05). Conclusion The 3D printing-assisted posterolateral approach is effective in the treatment of ankle fractures involving the posterior malleolus. The approach can be well planned before the operation, is simple to perform, yields good fracture reduction and fixation, and has good prospects for clinical application.

Zheng, Hongming↗

Fundamental mechanisms of tensile fracture in aluminum sheet undirectionally reinforced with boron filament

Results are presented from an experimental study of the tensile-fracture process in aluminum sheet unidirectionally reinforced with boron filament. The tensile strength of the material is severely limited by a noncumulative fracture mechanism which involves the initiation and sustenance of a chain reaction of filament fractures at a relatively low stress level. Matrix fracture follows in a completely ductile manner. The minimum filament stress for initiation of the fracture mechanism is shown to be approximately 1.17 GN/sq m (170 ksi), and appears to be independent of filament diameter, number of filament layers, and the strength of the filament-matrix bond. All the commonly observed features of tensile fracture surfaces are explained in terms of the observed noncumulative fracture mechanism.

Herring, H. W.↗