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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2020

Fracture toughness evaluation for Zr-4 clad tubing structure with pellet inserts

Applying fracture mechanics approach to spent nuclear fuel (SNF) system reliability investigation is warranted due to the inherent flaws and hydride structures existed in a SNF system after nuclear reactor operation. However, none of the existing fracture toughness data deal with fuel cladding specific geometry or spent fuel material conditions, such as cladding structure with the pellet-inserts and the associated pellet clad mechanical interactions induced mixed-mode damage mechanisms. Thus, the development of an intrinsic fracture mechanics approach that is suitable for SNF materials is needed. Due to thin wall and small dimension of clad tubing structure, the spiral notch torsion test (SNTT) method of small specimen approach was used to estimate the clad tubing structure fracture toughness. The estimated fracture toughness for Zr-4 cladding with alumina-pellet inserts are presented in this report. For SNTT samples with a short or medium crack length, between 5.4-mm and 8-mm, the estimate JIQ upon fracture initiation for the baseline Zr-4 cladding is at 50 kJ/m2 with 2-sigma uncertainty of 3.26 kJ/m2, and the associated KIQ is at 67.46 MPa√m. For SNTT samples with a long crack length, around 13-mm, the crack initialization is deviated from that of the Mode-I tensile fracture and appears to be a mixed-mode fracture of Mode I – tensile stress and Mode III – out of plane shear stress; the estimated JMQ is at 18.9 kJ/m2, the associated KMQ is at 41.4 MPa√m.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thin-Film Fracture Behavior for Diketopyrrolopyrrole Semiconducting Polymeric Films

Fracture energy, which quantifies a material’s resistance to the propagation of a pre-existing crack, is a key parameter for ensuring the mechanical reliability of stretchable organic electronic devices. However, most existing methods, such as a four-point bending fracture energy, utilized for measuring the fracture energy of semiconducting polymeric thin films are complicated by substrate effects, making it challenging to isolate the intrinsic behavior of the film from interfacial influences. In this study, we employed a pseudo free-standing pure shear method to systematically investigate the cohesive fracture energy of poly(diketopyrrolopyrrole-terthiophene) P(DPP-T)-based thin films to examine the effects of nanoconfinement, side chain length, degree of crystallinity, and strain rates. This method effectively eliminates substrate interference, enabling a direct assessment of the cohesive fracture energy of P(DPP-T) thin films. We found that thinner films and those with lower molecular weights exhibited significantly reduced fracture energies due to diminished chain entanglements. Additionally, films with shorter side chains displayed notably higher fracture energies, which were attributed to an increase in the degree of crystallinity. Finally, slower strain rates led to higher fracture energies, consistent with an enhanced stress relaxation. These insights offer practical guidelines for designing mechanically robust semiconducting polymers, contributing to the advancement of reliable, durable, flexible, and wearable electronic devices.

Zhang, Song [Univ. of Southern Mississippi, Hattie↗
DOE OSTI2022

Fracture toughness of schist, amphibolite, and rhyolite from the Sanford Underground Research Facility (SURF), Lead, South Dakota

The Cracked Chevron Notched Brazilian Disc (CCNBD) method was selected for Mode I fracture toughness tests on Poorman schist, Yates amphibolite, and rhyolite dikes from the EGS Collab site at the SURF in Lead, South Dakota. The effects of lithology, anisotropy, and loading rate were investigated. Fracture toughness was greatest in amphibolite, with schist and rhyolite having similar toughness values ($K$ amphibolite > $K$ rhyolite ≈ $K$ schist ). The effects of anisotropy on fracture toughness were investigated in the foliated schist samples. Schist samples were prepared in three geometries (divider, arrester, and short transverse) which controlled how the fracture would propagate relative to foliations. The divider geometry was strongest and short transverse geometry was the weakest ($K$ divider > $K$ arrester > $K$ shorttransverse ). Fracture toughness was observed to decrease with decreasing loading rate. Optical and SEM microscopy revealed that for the short transverse geometry, fractures tended to propagate along grain boundaries, whereas in arrester and divider geometries fractures tended to propagate through grains. In foliated samples, the tortuosity of the fracture observed in thin section was greater in arrester and divider geometries than in short transverse geometries.

58 GEOSCIENCES↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2024

Investigating In-Situ Fracture Behaviors of Polymer Pipeline Materials in Hydrogen and Hydrogen-Methane Blended Gas Environments

To reduce carbon emissions, the US natural gas infrastructure is seen as a primary solution for efficiently transporting hydrogen gas. Blending hydrogen gas with natural gas and transporting it across a national infrastructure could save significant infrastructure costs. To properly operate the infrastructure under the new gas system, it is critical to understand material compatibility with hydrogen under various conditions. The Blended Gas CRADA, a Hyblend project, is established to determine the material compatibility of existing natural gas pipes with hydrogen gas. In this study, we investigate the in-plane fracture behaviors of MDPEMarlex and HDPEGDB exposed to hydrogen and hydrogen-methane blended gas. Single-edge notch bending geometry is used. All tests are executed in-situ with the gas environment. The experimental results show a significant effect of the gas environment on HDPEGDB specimens, reducing 5% (H2) to 42% (Blended gas) of specific fracture energy compared to non-aged specimens. For the MDPEMarlex, the effects of the gas environment have increased the specific fracture energy by 10% (H2) to 15% (Blended gas). Fracture surfaces of the tested samples are observed using an electronic microscope. The in-plane fracture surface of HDPEGDB shows a pronounced dimple fracture pattern after exposure to hydrogen and blended gas. The expanded fracture pattern contributes to lower the specific fracture energy. These observations provide critical information for validating polymer pipeline materials when interact with hydrogen and hydrogen-blend gas.

Ko, Seunghyun↗
DOE OSTI2021

Investigating Fracture Network Deformation Using Noble Gas Release

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.

Gardner, W. Payton↗
DOE OSTI2021

Characterization of an electrically conductive proppant for fracture diagnostics

Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.

Geochemistry & Geophysics↗
DOE OSTI2020

CO2 injection with fracturing in geomechanically protected caprock: Task 6 of LLNL's Research Activities to Support DOE's Carbon Storage Program (FWP-FEW0191) (Final Technical Report)

Conventional principles of the design and operation of geologic carbon storage (GCS) require injecting CO 2 below the caprock fracturing pressure to ensure the integrity of the storage complex. In non-ideal storage reservoirs with relatively low permeability, modest injection rates can lead to pressure buildup and hydraulic fracturing of the reservoir and caprock. While the GCS community has generally viewed hydraulic fractures as a key risk to storage integrity, the actual behavior of a caprock hydraulic fracture, particularly from a geomechanical perspective, has not been thoroughly studied. The prevailing method of treating a hydraulic fracture as a high-permeability wing of the storage reservoir might have resulted in erroneous understanding of caprock hydraulic fracture behavior. Comprehensive analyses of monitoring data from the In Salah project suggested that one or more hydraulic fractures may have been created in the reservoir and lower caprock system during injection operations, but did not cause detectable leakage of CO 2 out of the storage complex. This observation coincides with a well-known phenomenon in the oil and gas industry: sedimentary rock formations have many inherent features that naturally protect the formations from unbounded vertical growth of hydraulic fractures.

58 GEOSCIENCES↗
DOE OSTI2022

Controlling Sustainability of Hydraulic Fracture Permeability in Ductile Shales

Hydraulic fracturing technology, along with horizontal drilling and associated completion technologies, dramatically increased US oil and gas production, by increasing the permeability of tight shales and allowing for cost-effective extraction of trapped hydrocarbons. For efficient and economical production, the targeted shale has to be fractured easily (i.e., good ‘frackability’), but also the permeability of induced fractures have to be sustained during production (good ‘sustainability’). For these reasons, currently, ductile shale with high clay content (>~40%) is difficult to exploit as a hydrocarbon resource, although hydrocarbons can still be found in it. Good frackability, high-TOC, low-clay-content resource shales—the ‘low-hanging fruit’—will be depleted, and we need to develop tools that allow tapping into currently underutilized, clay-rich ductile shales. This research project aims (1) to understand the behavior of fractures in clay-rich, ductile (and sometimes swelling) shales and (2) to begin to develop technologies for efficient and economical production from such shales. With these objectives, we examine time-dependent, coupled mechanical-hydrological behavior of open and proppant-filled fractures within different types of shales. In the preceding project, the focus of the research was to understand the behavior of the fractures and to predict their permeability changes due to fracture closure from shale matrix deformation and proppant embedment. In the current project, we investigate the possibility of chemically manipulating the sustainability of hydraulic fractures in ductile shales—specifically by altering the proppant-embedding behavior—via chemical means, at different stages of oil and gas production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Fracture resistance of vintage cast iron in gaseous hydrogen

In an effort to decarbonize legacy energy systems, several projects around the world are exploring alternatives to natural gas. Gaseous hydrogen is proposed as a carbon-free fuel to displace natural gas in existing legacy natural gas distribution systems, some of which continue to operate after 100 years (or more) in service. These systems, particularly in older industrial centers, contain cast iron pipe. However, the fracture resistance of most metals is degraded in gaseous hydrogen environments. This study evaluated the fracture resistance of several legacy cast iron pipe materials while exposed to gaseous hydrogen. Measurements were performed in three environments: air, a gas mixture with hydrogen partial pressure of 1 bar and pure hydrogen with a partial pressure of 34 bar. Although cast iron is generally considered a low ductility metal, elastic-plastic fracture methods are needed to assess the fracture resistance of the relatively small specimens that can be extracted from legacy pipe. Hydrogen reduced the fracture resistance of these cast iron materials by 10-40%. In air, the fracture resistance was determined to be as high as 21 MPa m 1/2 , whereas in gaseous hydrogen at pressure of about 1 bar the fracture resistance was as low as 13 MPa m 1/2 . Additional modest degradation of the fracture resistance was assessed at higher partial pressure (as low as 12 MPa m 1/2 ).

08 HYDROGEN↗
DOE OSTI2026

Development, Monitoring, and Control of Fracture Thermal Energy Storage (FTES) in Crystalline Rock Formations (DEMO-FTES) (CRADA Final Report)

The DEMO-FTES project sought to demonstrate the thermal efficiency of fracture thermal energy storage (FTES) through numerical simulations, laboratory and meso-scale field tests. A detailed dimensional and scaling analysis was performed to identify key parameters and how they can be most effectively scaled to the laboratory and decameter scale. Numerical models were developed and used for three purposes: 1. Before field experiments, numerical modelling can be used to estimate fracture properties based on previous data from the EGS Collab experiment and then predict thermal hydrological behaviors of the fracture system with hot water injection/withdrawal, therefore, help to design the experiments (e.g., to decide the duration of the cycles based on the flow rate the pump can provide, and the estimated fracture properties); 2. After the field experiment, to estimate the system properties during the experiment (as the size and shape of a fracture could change over time), and understand why system performance is different than what has been predicted, i.e., to help understand the meso-scale test; and, 3. To model the lab experiments and estimate fracture properties and storage efficiency. Ultimately, the experiment and numerical models could shed light on the processes and uncertainty happening during fracture activation and help understand the scaling between lab and field test, and finally, the design and optimization of potential fracture thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Fracture Mechanics—Theory, Modeling and Applications

The field of fracture mechanics was developed during the throes of World War II, and since then, it has been a very active area of research. This is a very energetic and vibrant field because fracture processes are relevant for a vast range of engineering applications spanning different temporal and length scales. Depending on the application, material fracture and fragmentation may be either a desirable feature (for example, for the purpose of hydraulic fracturing operations), or something that is to be avoided at all costs (for example, in key components of mechanical systems). In some other instances, fractures are already present and must be accommodated through novel design (for example, in the case of earthquake ruptures occurring in the Earth’s crust). Because of this wide spectra of motivation, there is a continuous need for the fracture mechanics research community to gain more insight on how fractures are created, evolve, and interact with themselves and with material microstructures. In addition, there are still many remaining questions about how fracture propagation, growth, and interaction impact overall material behavior, and how to develop accurate models to predict this behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Insights into the hydro-mechanical behavior of a decimeter-scale fracture using the mini-SIMFIP probe

ABSTRACT: Understanding hydro-mechanical couplings in fractured rocks is essential for predicting the rock mass response during high-pressure fluid injection, including the stimulation of enhanced geothermal systems. However, fluid-driven fracture dislocations are challenging to measure due to the need for local displacement data at high fluid pressures. In this study, fluid-driven displacement across a decimeter-scale laboratory fracture was investigated using the mini-SIMFIP (step rate injection method for fracture in-situ properties) probe, which is a smaller version of the SIMFIP tool (Guglielmi et al., 2014). The mini-SIMFIP probe is able to resolve the full 3D displacement vector of a fracture. The probe was installed in one of two boreholes across a decimeter-scale saw-cut granite fracture. Two pressure step injection tests were conducted under the same isotropic stress conditions. By varying injection between the boreholes, we estimated the aperture profile across the fracture. The results show a consistent pressure-dependent opening as long as steady-state flow was maintained. At a certain pressure step, steady-state conditions were no longer achievable and the fracture opened rapidly. The pressure-opening relationship diverged between the two tests and indicated a homogenization of the aperture profile and an increase in the non-linearity of the flow regime.

Osten, J↗
DOE OSTI2020

Study of Foam Fracturing for Enhanced Geothermal Systems Using Model Material

Foam fracturing is considered a potential approach to address water concerns with hydraulic fracturing in the development of enhanced geothermal systems (EGS). In many EGS sites, water required for hydraulic fracturing is either unavailable, extremely costly or environmentally unsustainable. This paper presents work performed in the first year of a project sponsored by the U.S. DOE GTO Waterless Stimulation Initiative that investigates foam fracturing as an alternative approach to hydraulic fracturing. One element of this work explores foam formulations that are suitable for EGS conditions. A separate paper submitted to this workshop covers the relevant work on foam characterization (Thakore et al., 2020).This paper describes the development of an experimental setup for studying the applicability of foam fluids for hydraulic fracturing along with enhanced fracture efficiency through cyclic pressurization. The system is equipped with pulse rate automation, dual pulse valves for enhanced pulsation amplitudes, and other functions. It can perform static injection at a given pressurization rate, and pulsed injection at specified pulse size and rate. Both single-phase and foamed fluids can be used. Initial experimental results of foam fracturing using cement as a model material are reported. The observations from the experimental work will be presented and discussed.

Wang, Hong↗
DOE OSTI2020

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗
DOE OSTI2021

Study on Foam Fracturing of Granite for the Development of Enhanced Geothermal Systems

Foam fracturing has been considered as a potential technical approach to address water concerns with hydraulic stimulation in the development of enhanced geothermal systems (EGS). This paper reports the recent progress of a project sponsored by the U.S. DOE GTO Waterless Stimulation Initiative. The project investigates the foam fracturing as an alternative approach to the hydraulic fracturing. The paper describes the development of foam fracturing testing system at Oak Ridge National Laboratory (ORNL). The system can be used to perform the hydraulic fracturing of rock-like materials with both water and foamed fluids under pressure to 6,000 psi (41.4 MPa). The system features the sections of foam generation and foam injection and has the capability of testing materials in both monotonic and cyclic injection modes.Experimental results of foam fracturing on Charcoal Black granite will be reported. Cylindrical granite specimens with a blind hole were studied along with water and nitrogen-gas-in-water foam as the fracturing fluids. The effects of fracturing fluids and injection modes on the breakdown pressure and failure responses of the material were examined. The observations from the experimental work and the implications to the EGS application will be presented and discussed.

Wang, Hong↗
DOE OSTI2023

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗