SEARCH · Engineering Papers

Results for “FRACTURES”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Relationships between lineal fracture intensity and chemical composition in the Marcellus Shale, Appalachian Basin

Within mudrock reservoirs, brittle zones undergo failure during hydraulic stimulation, creating numerous artificial fractures which enable hydrocarbons to be liberated from the reservoir. Natural fractures in mudrock reduce the tensile strength of the host rock, creating planes of weaknesses that are hypothesized to be reactivated during hydraulic stimulation. Combined, brittleness and natural fractures contribute to creating more abundant and complex fracture networks during hydraulic stimulation. Research efforts toward quantifying rock brittleness have resulted in numerous mineral-/compositional-based indices, which are used during petrophysical analysis to predict zones most conducive to hydraulic stimulation. In contrast, investigations on the relationship between chemical composition and core-scale natural fractures are limited. For this study, we collected high-resolution energy-dispersive X-ray fluorescence (XRF) data, calibrated with a wave-dispersive XRF, from a Marcellus Shale core. Additionally, we characterized corescale natural fractures in terms of length, width, in-filling material or lack thereof, and orientation. Following the characterization, we transformed the natural fracture data into a continuous P10 (lineal fracture intensity) curve, expressed as the number of fractures per a one-half foot window. Using these data sets, we investigated the relationship between rock composition and natural fracture intensity. Regression analyses recorded positive relationships between natural fracture intensity and calcium, silicon/aluminum, and total organic carbon (TOC), and negative relationships with silicon and aluminum. Aluminum recorded the strongest (negative) relationship ([Formula: see text]) with natural fracture intensity. To access the degree to which natural fractures can be predicted based on chemical composition, we applied a partial least-squares analysis, a multivariate method, and recorded an [Formula: see text]. Our study illustrates that although numerous factors are responsible for natural fracture genesis, such fractures predictively concentrate in areas of similar chemical composition, largely in zones with low aluminum concentrations.

Geochemistry & Geophysics↗

DOE OSTI2021

Fracture toughness evaluations for spent nuclear fuel dry storage canister welds and spent nuclear fuel clad-pellet structures

Integrity of spent nuclear fuel (SNF) interim storage canisters is very important to the safety of the back-end nuclear fuel cycle. Stress corrosion cracking (SCC) potential of interim storage canister has been considered as a high priority. Because no post-weld heat treatment was required for forming these canisters, the high tensile residual stress existed within these canister welds. This can change the fracture resistance capacity significantly as well as increase SCC potential. Due to relative thin shell thickness of a canister weldment, the spiral notch torsion test (SNTT) method was used to estimate the canister weldment fracture toughness. SNTT was developed to measure the intrinsic fracture toughness of structural materials using small specimens. The SNTT method has been applied to a wide variety of structural materials, such as low-alloy steels, stainless steel, aluminum alloy, ceramics, concrete, and composites. The SNTT system operates by applying pure torsion to cylindrical specimens with a notch line that spirals around the specimen at a 45° pitch. In order to carry out pure torsion load mode, biaxial tension/torsion tester was developed accordingly to perform SNTT protocol. Moreover, applying fracture mechanics approach to 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. Furthermore, due to thin wall and small dimension of clad tubing structure, the SNTT method was used to estimate the clad tubing structure fracture toughness. Fracture testing were performed on the received stainless steel canister weldment, most SNTT weld samples fracture initiation sites are at heat-affected zone (HAZ) regions. The estimated fracture toughness J Q ’ for the baseline SS304 steel is at 283 kJ/m². The estimated JQ’ for the SS304/308 weld and baseline metals are 148 kJ/m 2 and 459 kJ/m 2 , respectively. Out of cell fracture testing for spent fuel structure were carried out on the surrogate rods made of Zr-4 clad and alumina inserts, the estimated fracture toughness values for baseline Zr-4 cladding with alumina-pellet inserts are: (1) For SNTT samples with a short or medium crack length, between 5.4-mm and 8-mm, the estimate J IQ upon fracture initiation for the baseline Zr-4 cladding is at 50 kJ/m 2 with 2-sigma uncertainty of 3.26 kJ/m 2 , and the associated K IQ is at 67.46 MPa$\sqrt{m}$; and (2) 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 J MQ is at 18.9 kJ/m 2 , the associated K MQ is at 41.4 MPa$\sqrt{m}$.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗

DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗

DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗

DOE OSTI2015

Stimuli-Responsive/Rheoreversible Hydraulic Fracturing Fluids as a Greener Alternative to Support Geothermal and Fossil Energy Production

Cost-effective yet safe creation of high-permeability reservoirs within deep bedrock is the primary challenge for the viability of enhanced geothermal systems (EGS) and unconventional oil/gas recovery. Although fracturing fluids are commonly used for oil/gas, standard fracturing methods are not developed or proven for EGS temperatures and pressures. Furthermore, the environmental impacts of currently used fracturing methods are only recently being determined. Widespread concerns about the environmental contamination have resulted in a number of regulations for fracturing fluids advocating for greener fracturing processes. To enable EGS feasibility and lessen environmental impact of reservoir stimulation, an environmentally benign, CO2-activated, rheoreversible fracturing fluid that enhances permeability through fracturing (at significantly lower effective stress than standard fracturing fluids) due to in situ volume expansion and gel formation is investigated herein. The chemical mechanism, stability, phase-change behavior, and rheology for a novel polyallylamine (PAA)-CO2 fracturing fluid was characterized at EGS temperatures and pressures. Hydrogel is formed upon reaction with CO2 and this process is reversible (via CO2 depressurization or solubilizing with a mild acid) allowing removal from the formation and recycling, decreasing environmental impact. Rock obtained from the Coso geothermal field was fractured in laboratory experiments under various EGS temperatures and pressures with comparison to standard fracturing fluids, and the fractures were characterized with imaging, permeability measurement, and flow modeling. This novel fracturing fluid and process may vastly reduce water usage and the environmental impact of fracturing practices and effectively make EGS production and unconventional oil/gas exploitation cost-effective and cleaner.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2014

Fractured rock stress‐permeability relationships from in situ data and effects of temperature and chemical‐mechanical couplings

Abstract The purpose of this paper is to (i) review field data on stress‐induced permeability changes in fractured rock; (ii) describe estimation of fractured rock stress‐permeability relationships through model calibration against such field data; and (iii) discuss observations of temperature and chemically mediated fracture closure and its effect on fractured rock permeability. The field data that are reviewed include in situ block experiments, excavation‐induced changes in permeability around tunnels, borehole injection experiments, depth (and stress) dependent permeability, and permeability changes associated with a large‐scale rock‐mass heating experiment. Data show how the stress‐permeability relationship of fractured rock very much depends on local in situ conditions, such as fracture shear offset and fracture infilling by mineral precipitation. Field and laboratory experiments involving temperature have shown significant temperature‐driven fracture closure even under constant stress. Such temperature‐driven fracture closure has been described as thermal overclosure and relates to better fitting of opposing fracture surfaces at high temperatures, or is attributed to chemically mediated fracture closure related to pressure solution (and compaction) of stressed fracture surface asperities. Back‐calculated stress‐permeability relationships from field data may implicitly account for such effects, but the relative contribution of purely thermal‐mechanical and chemically mediated changes is difficult to isolate. Therefore, it is concluded that further laboratory and in situ experiments are needed to increase the knowledge of the true mechanisms behind thermally driven fracture closure, and to further assess the importance of chemical‐mechanical coupling for the long‐term evolution of fractured rock permeability.

58 GEOSCIENCES↗

DOE OSTI2013

Evaluation of known-boundary and resistivity constraints for improving cross-borehole DC electrical resistivity imaging of discrete fractures

There is a need to better characterize discrete fractures in contaminated hard rock aquifers to determine the fate of remediation injections away from boreholes and also to evaluate hydraulic fracturing performance. A synthetic cross borehole electrical resistivity study was conducted assuming a discrete fracture model of an existing contaminated site with known fracture locations. Four boreholes and two discrete fracture zones, assumed to be the dominant electrical and hydraulically conductive pathways, were explicitly modeled within an unstructured tetrahedral mesh. We first evaluated different regularization constraints starting with an uninformed smoothness-constrained inversion, to which a priori information was incrementally added. We found major improvements when (1) smoothness regularization constraints were relaxed (or disconnected) along boreholes and fractures, (2) a homogeneous conductivity was assumed along boreholes, and (3) borehole conductivity constraints that could be determined from a specific conductance log were applied. We also evaluated the effect of including borehole packers on fracture zone model recovery. We found the fracture zone conductivities with the inclusion of packers were comparable to similar trials excluding the use of packers regardless of electrical potential changes. The misplacement of fracture regularization disconnects (FRDs) can easily be misinterpreted as actual fracture locations. Conductivities within these misplaced disconnects were near the starting model value and removing smoothing between boreholes and assumed fracture locations helped in identifying incorrectly located FRDs. Analysis confirmed improved model sensitivity when regularization disconnects were (1) applied along the boreholes and fracture zones, and (2) FRDs were placed where actual fractures existed. Our findings demonstrate that structural constraints used after careful evaluation of a priori information are critical to improve imaging of fracture electrical conductivities, locations and orientations.

geophysics, inversion, resistivity, fracture, cont↗

DOE OSTI2022

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

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

02 PETROLEUM↗

DOE OSTI2020

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

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

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

DOE OSTI2012

Permeability Calculation in a Fracture Network - 12197

Laminar flow of a viscous fluid in the pore space of a saturated fractured rock medium is considered to calculate the effective permeability of the medium. The effective permeability is determined from the flow field which is calculated numerically by using the finite element method. The computation of permeability components is carried out with a few different discretizations for a number of fracture arrangements. Various features such as flow field in the fracture channels, the convergence of permeability, and the variation of permeability among different fracture networks are discussed. The longitudinal permeability in general appears greater than the transverse ones. The former shows minor variations with fracture arrangement whereas the latter appears to be more sensitive to the arrangement. From the calculations of the permeability in a rock medium with a fracture network (two parallel fractures aligned in the direction of 45-deg counterclockwise from the horizontal and two connecting fractures(narrowing, parallel and widening) the following conclusions are drawn. 1. The permeability of fractured medium not only depends on the primary orientation of the main fractures but also is noticeably influenced by the connecting fractures in the medium. 2. The transverse permeability (the permeability in the direction normal to the direction of the externally imposed macro-scale pressure gradient) is only a fraction of the longitudinal one, but is sensitive to the arrangement of the connecting fractures. 3. It is important to figure out the pattern of the fractures that connect (or cross) the main fractures for reliable calculation of the transverse permeability. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Computational Modeling of Hydraulic Properties of a Sheared Single Rock Fracture

Flow through a single mechanically sheared Marcellus shale fracture was investigated computationally and experimentally. To provide a better understanding of the variation of hydraulic and geometrical characteristics of a fracture subjected to shearing, coupled shear flow tests on the fracture for four shearing displacement steps under constant normal stress were performed. At the end of each shearing step, computed tomography (CT) scans with resolution 26.8 μm were obtained and the corresponding fracture geometries were evaluated. The CT images were used to generate full aperture maps of the fracture configuration. In addition, average aperture maps were also created by averaging the full-resolution data over 10 × 10 pixels, smoothing out fine structural details. Computational modeling of water flow through the fractures at different shearing steps was performed using a modified local cubic law approach and the 3D full Navier–Stokes equations with the use of the ANSYS-Fluent software. Both the average aperture maps and full maps were used in these simulations. The experimental pressure drops of the fracture at shearing step 1, which has very small apertures, poorly matched the numerical results, quite likely because the fracture structure was inadequately captured by the scanning resolution. Shearing typically increased the aperture height of the fracture, whose features were then better captured by the CT scan. Good agreement between the experimental data and the numerical results of the full map for shearing step 2 was observed. The simulations were performed for both full and average aperture maps, and the effects of scan resolution and surface roughness on the accuracy of the results were studied. The modified local cubic law and full Navier–Stokes simulations of the averaged map fracture were found to be in good agreement. It was conjectured that this was because the nonlinear losses were insignificant for the smoothed out averaged map fracture. Similar comparisons with those of the full map showed agreement in trends, but there were some quantitative differences. The averaged fracture map simulations also predicted lower pressure drops compared to the full map, particularly for high flow rates. These differences were due to the fine-scale geometrical complexity (surface roughness) of fracture geometry that affects the fluid flow in the fracture. In conclusion, an improved cubic law model was also proposed, and its accuracy was verified by comparing its predictions with those of the Navier–Stokes simulations.

58 GEOSCIENCES↗

DOE OSTI2023

Using in-situ strain measurements to evaluate the accuracy of stress estimation procedures from fracture injection/shut-in tests

Fracture injection/shut-in tests are commonly used to measure the state of stress. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of Hayashi and Haimson. In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. In conclusion, the calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

58 GEOSCIENCES↗

DOE OSTI2018

Injection‐Induced Shear Slip and Permeability Enhancement in Granite Fractures

Abstract Shear stimulation, so‐called hydroshearing , is believed to reactivate critical or near‐critical preexisting fractures by pressurized fluid injection, causing them to slip and dilate for economic production from geothermal and unconventional petroleum reservoirs. However, little or no experiments have been conducted to directly probe the coupled hydromechanical responses and permeability evolution during fracture shearing. In this paper, we present the results of novel injection‐induced shear tests on cylindrical granite samples each containing a single tensile or saw cut fracture. The shear test on rough fractures demonstrates that retainable permeability enhancement can be achieved through dilatant shear slip. The characteristic stick‐slip behavior of the fracture during shearing is also observed. Fracture aperture controlled by effective normal stress dominates the fluid flow in stick state, while in slip state, the major contributor to production increase is the irreversible normal dilation caused by fracture shear slip. The experimental data show that depending on the fracture roughness, the fracture slip process includes two quasi‐static slip intervals with a slip velocity of ~10 −9 to ~10 −6 m/s and a dynamic slip interval with a slip velocity of ~10 −7 to ~10 −5 m/s. The fracture slip correlates well with the associated stress relaxation: a faster fracture slip induces a quicker stress relaxation, and vice versa. The influence of surface roughness on the hydromechanical responses during fracture shearing and the shear‐induced asperity degradation are also studied. The various experimental observations on rough and smooth fractures clearly indicate that the key component to shear stimulation is fracture self‐propping by asperities.

Ye, Zhi↗

DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in theMidwest United States

This report presents final technical results for the project Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in the Midwest United States (DE-FE0023330). The project was a three-year effort consisting of seven technical tasks focused on defining geomechanical factors for CO 2 storage applications in deep saline rock formations in Ohio and the Midwest United States, because geomechancial issues have been identified as a significant risk factor for large-scale CO 2 storage applications. A basin-scale stress-strain analysis was completed to describe the geomechanical setting for rock formations of Ordovician-Cambrian age in Ohio and adjacent areas of the Midwest United States in relation to geologic CO 2 storage applications. The tectonic setting, stress orientation-magnitude, and geomechanical and petrophysical parameters for CO 2 storage zones and caprocks in the region were cataloged. Ten geophysical image logs were analyzed for natural fractures, borehole breakouts, and drilling-induced fractures. The logs indicated mostly less than 10 fractures per 100 vertical feet in the borehole, with mostly N65E principal stress orientation through the section. Geophysical image logs and other logs were obtained for three wells located near the sites where specific models were developed for geomechanical simulations: Arches site in Boone County, Kentucky; Northern Appalachian Basin site in Chautauqua County, New York; and E-Central Appalachian Basin site in Tuscarawas County, Ohio. For these three wells, 9,700 feet of image logs were processed and interpreted to provide a systematic review of the distribution within each well of natural fractures, wellbore breakouts, faults, and drilling induced fractures. There were many borehole breakouts and drilling-induced tensile fractures but few natural fractures. Concentrated fractures were present at the Rome-basal sandstone and basal sandstone-Precambrian contacts at the Arches and East-Central Appalachian Basin sites. Geophysical logs were utilized to develop local-scale geologic models by determining geomechanical and petrophysical parameters within the geologic formations. These data were ported to coupled fluid-flow and reservoir geomechanics multi-phase CO 2 injection simulations. The models were developed to emphasize the geomechanical layers within the CO 2 storage zones and caprocks. A series of simulations were completed for each site to evaluate whether commercial-scale CO 2 could be safely injected into each site, given site-specific geologic and geomechanical controls. This involved analyzing the simulation results for the integrity of the caprock, intermediate, and reservoir zones, as well quantifying the areal uplift at the surface. Simulation results were also examined to ensure that the stress-stress perturbations were isolated within the subsurface, and that there was only limited upward migration of the CO 2 . Simulations showed capacity to inject more than 10 million metric tons of CO 2 in a single well at the Arches and East Central Appalachian Basin sites without excessive geomechanical risks. Low-permeability rock layers at the Northern Appalachian Basin study area well resulted in very low CO 2 injection capacity. Fracture models developed for the sites suggests that the sites have sparse fracture network in the deeper Cambrian rocks. However, there were indicators in image logs of a moderate fracture matrix in the Rose Run Sandstone at the Northern Appalachian Basin site. Dual permeability fracture matrix simulations suggest the much higher injection rates may be feasible in the fractured interval. Guidance was developed for geomechanical site characterization in the areas of geophysical logging, rock core testing, well testing, and site monitoring. The guidance demonstrates that there is a suitable array of options for addressing geomechanical issues at CO 2 storage sites. Finally, a review of Marcellus and Utica-Point Pleasant shale gas wells and CO 2 storage intervals indicates that these items are vertically separated, except for the Oriskany sandstone and Marcellus wells in southwest Pennsylvania and northern West Virginia. Together, project results present a more realistic portrayal of geomechanical risk factors related to CO 2 storage for existing and future coal-fired power plants in Ohio.

58 GEOSCIENCES↗

DOE OSTI2017

Report on Developing a US-China Joint Project on CO2-Based Fracturing Techniques

There are several types of hydraulic fracturing methods that have been investigated in the past decades that use little or no water. The criteria for selecting an optimal hydraulic fracturing fluids are many, and no single fluid can meet all the desired characteristics. Developing a set of technologies by using different fracturing fluids can better meet individual geologic, economic, and environmental needs. CO2 as an alternative of water for hydraulic fracturing has its unique place among all other fracturing fluids, and can meet both economic and environmental objectives under certain circumstances. There are three main categories of CO2-based fracturing fluids, simple liquid CO2 (LCO2), directly thickened LCO2, and CO2-foams/emulsions. Although LCO2 has many advantages, critical disadvantages arise due to its high operating pressure, low viscosity, poor proppant-carrying capacity, and generation of shallower fracture networks. To date, techniques for direct thickening of LCO2 to overcome these limitations have proven ineffective. To increase the apparent viscosity of CO2, we propose to further develop CO2-foam (and emulsion) fracturing techniques. While CO2-foams have long been used in EOR, the required properties for hydraulic fracturing are different, and relatively few peer-reviewed publications and field test reports are available regarding their applications in fracturing, and the surfactants and additives used are often proprietary This report describes a brief state-of-the-art analysis of CO2-based hydraulic fracturing methods and then discusses a joint plan for a long-term R&D collaboration to advance the most promising CO2-based hydraulic fracturing techniques as part of the U.S.-China Clean Energy Research Center, Water Energy Technology (CERC-WET). The objective of this U.S.-China joint research project is to design and laboratory test a CO2-foam fracturing fluid system, and then conduct pilot field-tests on them in China’s Jilin Oilfield. This joint research project will benefit both U.S. and China for advancing their shared goal of developing better technologies for CO2-based fracturing. The Chinese CERC-WET has two teams joining the project, one from RIPED, responsible for aboveground infrastructure design and construction. Another team from the Jilin Oilfield will be responsible for the field test design and conducting the test. The U.S.CERC-WET team at LBNL will be responsible for providing the CO2-foam formulations and their rheology parameters to be tested in the field (a challenging task and a key to success); participating in the equipment design and field test design; assisting in analyzing field data; and reporting to U.S. DOE and CERC-WET.

58 GEOSCIENCES↗

NASA NTRS2017

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗