SEARCH · Search NASA
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.
This is the conference paper accompanying an oral presentation at the 58th US Rock Mechanics/Geomechanics Symposium (ARMA) June 23 to 26 in Golden, CO. The paper highlights using machine learning to quantify fracture network for carbon storage and is part of the SMART Initiative project applications.
This is an invited guest speaker's presentation. The present covers three use cases by applying machine learning techniques. The use cases include fracture analysis for CCUS: IBDP study, multiple level of fracture network analysis and tool: HFTS1 study, Frac-Hit with Middleland Basin datasets from collaborations with Company A.
This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).
Fracture mechanics calculations using the finite element method can be computationally expensive, which makes them challenging to use in engineering evaluations of flaws in pressure vessels. The weight function (WF) technique is a reduced order fracture modeling method that is widely used to greatly reduce these computational costs. Its computational efficiency is essential for use in probabilistic fracture mechanics evaluations of embrittled nuclear reactor pressure vessels (RPVs), due to the large number of sampled flaws that must be evaluated. Although the WF technique is general, it is typically used only for axis-aligned flaws. Recent discoveries of off-axis flaws in operating nuclear reactors have necessitated detailed simulations of such flaws and the interactions between neighboring flaws. This paper presents a generalized WF approach applicable for analyzing arbitrary flaw geometries in thick-walled cylindrical vessels, including surface-breaking and subsurface flaws, which can either axis-aligned or off-axis, and can account for interactions with other flaws. Herein, this approach is demonstrated on representative simulations of multiple flaw geometries in a RPV subjected to transient loading conditions. In all cases, the WF technique gives good comparison with benchmark results from direct simulation with greatly reduced computational effort.
A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.
To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.
Presentation at the 58th US Rock Mechanics/Geomechanics Symposium (ARMA) June 23 to 26 in Golden, CO. The presentation highlights using machine learning to quantify fracture network for carbon storage and is part of the SMART Initiative project applications.
Not provided.
This is an invited talk for Alberta Energy Regulator, Canada, for machine learning applications for CCS, unconventional hydrocarbon, and frac-hit identification based on NETL recent activities (Virtual, October 1, 2024).
The High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is approaching the production phase of the US-contributed Q1 and Q3 Interaction Region Quadrupoles (MQXFA). The structures for the MQXFA prototypes were design and inspected by the US-LARP (LHC Accelerator Research Program), AUP developed criteria, which will be used for the pre-series structures. As the first two full-length prototypes with 4.2 m magnetic length, MQXFAP1 and MQXFAP2, were designed and assembled at Lawrence Berkeley National Laboratory (LBNL), and tested at Brookhaven National Laboratory (BNL). The end aluminum short shell of MQXFAP2 was fractured along the shell length during the test, and tests were stopped. Analytical and Finite Element analysis were performed in light of the graded procedure defined in the Structure Design Criteria to investigate the fracture failure for MQXFAP2. In this paper, we report the fracture analysis of the current shell design, including the elasto-plastic simulations with submodel technique, and calculations with Linear Elastic Fracture Mechanics (LEFM). Test material properties are also presented. The results of this analysis explain why the end shell of MQXFAP2 failed, and suggest fillets on the end shell notches to meet the margin specified in the Structural Design Criteria.
This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.
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.
Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.
Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.
Abstract A tensile properties testing study was conducted to understand the influence of thickness, cross‐head speed (speed of testing), gauge length (GL; specimen test length), and sample shape on important tensile properties of polyvinyl alcohol (PVA) nanofiber webs. The effects of each testing parameter on load at break, extension at break, Young's modulus, and tensile stress–strain curve of PVA nanofiber webs are analyzed. The Welch two sample t ‐tests show the significant difference among tested data. Using interaction plots, two‐way analysis of variance, and margin mean plots, the interaction effects among testing parameters have been analyzed. Of all the factors, cross‐head speed, the interaction among GL, and sample thickness (GL: Thickness) and the interaction among GL, testing speed and sample thickness (GL: Speed: Thickness) have significant influence on the tensile properties of PVA nanofiber webs. Moreover, the hypothesized model of mechanism of tensile strain–stress curve of PVA nanofiber webs has been proposed. Based on the model, the tensile strain–stress curve can be split into three stages: linear elastic, partial break up, and complete breakage. This study will provide a better understanding of tensile testing parameters' effects and their interaction effects on the tensile properties of nanowebs.
The deformation and failure mechanisms in an ultrasonic additively manufactured Al-6061 alloy are studied using in-situ x-ray computed tomography during tensile deformation as well as ex-situ characterization with scanning electron microscopy. The role of pre-existing voids between the build layers on the failure mechanism is discussed with a focus on individual foil failures and void coalescence.
Explore the source record for details and available documents.
Explore the source record for details and available documents.