SEARCH · Engineering Papers

Results for “fracturing fluid”

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.

607 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗

NASA NTRS2003

Mineralogy and Petrology of Amoeboid Olivine Inclusions in CO3 Chondrites: Relationship to Parent-Body Aqueous Alteration

Petrographic and mineralogic studies of amoeboid olivine inclusions (AOIs) in CO3 carbonaceous chondrites reveal that they are sensitive indicators of parent-body aqueous and thermal alteration. As the petrologic subtype increases from 3.0 to 3.8, forsteritic olivine (Fa(sub 0-1)) is systematically converted into ferroan olivine (Fa(sub 60-75)). We infer that the Fe, Si and O entered the assemblage along grain boundaries, forming ferroan olivine that filled fractures and voids. As temperatures increased, Fe(+2) from the new olivine exchanged with Mg(+2) from the original AOI to form diffusive haloes around low-FeO cores. Cations of Mn(+2), Ca(+2) and Cr(+3) were also mobilized. The systematic changes in AOI textures and olivine compositional distributions can be used to refine the classification of CO3 chondrites into subtypes. In subtype 3.0, olivine occurs as small forsterite grains (Fa(sub 0-1)), free of ferroan olivine. In petrologic subtype 3.2, narrow veins of FeO-rich olivine have formed at forsterite grain boundaries. With increasing alteration, these veins thicken to form zones of ferroan olivine at the outside AOI margin and within the AOI interior. By subtype 3.7, there is a fairly broad olivine compositional distribution in the range Fa(sub 63-70), and by subtype 3.8, no forsterite remains and the high-Fa peak has narrowed, Fa(sub 64-67). Even at this stage, there is incomplete equilibration in the chondrite as a whole (e.g., data for coarse olivine grains in Isna (CO3.8) chondrules and lithic clasts show a peak at Fa(sub39)). We infer that the mineral changes in A01 identified in the low petrologic types required aqueous or hydrothermal fluids whereas those in subtypes greater than or equal to 3.3 largely reflect diffusive exchange within and between mineral grains without the aid of fluids.

Chizmadia, Lysa J.↗

DOE OSTI2020

Blunt Impact Brain Injury using Cellular Injury Criterion

The Advanced Combat Helmet (\ACH") military specification (\mil-spec") requires a helmeted magnesium (\Mg") Department of Transportation (\DOT") headform be dropped vertically, with an impact speed of 3.1 m/s (10 ft/s), onto a steel hemispherical target. The pass/fail criteria are based on translational acceleration (150 G) alone, absent of any rotational component. Without a rotational component, the specification's injury risk application is limited to skull fracture and peripheral hematomas (subdural, subarachnoid), since this translational acceleration injury risk assessment is based on the Wayne State Tolerance Curve (\WSTC"). To provide a more comprehensive view of injury for the entire brain, an alternative approach is needed. To meet this need, we worked with a larger group called PANTHER, a collaboration between national laboratories, industry, and academia. Collaborations specific to research and results presented here come from efforts led by Mr. Ron Szalkowski and Mr. Sushant Malave, Ms. Alice Fawzi, and Dr. Christian Franck. We have developed a prototypical injury risk criterion based on the neuronal response to abrupt changes in general motion (translation, rotation, or both). The cellular-based mild traumatic brain injury (\cbmTBI") criterion utilizes both the strain and strain rate of brain tissue to account for the stretch and rate of stretch that occurs throughout the brain as a result of blunt impact to the head. We conducted physical experiments of an ACH-fitted magnesium headform, which produced repeatable headform peak accelerations. Then, we developed a simulation of the experiment, and validated the simulation output with the experimental data. We then substituted the magnesium headform with a human headform, consisting of skin, muscle, bone, gray matter, white matter, cerebral-spinal fluid, membranes, vasculature, intravertebral discs, airway and sinus. We quantified brain injury risk using the cbmTBI criterion, using the current mil-spec test and a modified test. The modified mil-spec test used an inclined anvil target that was located posterior to the crown of the helmet in the axial plane. While the current mil-spec test produced brain deformation from head translation alone, the modified test produced brain deformation from head translation and rotation, which is closer to most real world and combat theater impacts (e.g., such as occur in tertiary blast exposure). Compared to the current mil-spec test, the modified test produced elevated strains in the human digital twin. These data, mapped to the cbmTBI criterion, suggest increased injury risk for blunt impacts that cause rotation and translation, rather than just translation alone. Moreover, these data may lead to a rotational performance metric, which is rooted in the biology and pathology of the brain's response to impact and blast, and which should be used to improve next-generation helmet designs.

42 ENGINEERING↗

NASA NTRS1993

Double Cantilever Beam and End Notched Flexure Fracture Toughness Testing of Two Composite Materials

Two different unidirectional composite materials were provided by NASA Langley Research Center and tested by the Composite Materials Research Group within the Department of Mechanical Engineering at the University of Wyoming. Double cantilever beam and end notched flexure tests were performed to measure the mode I (crack opening) and mode II (sliding or shear) interlaminar fracture toughness of the two materials. The two composites consisted of IM7 carbon fiber combined with either RP46 resin toughened with special formulation of LaRC IA resin, known as JJS1356; or PES chain extended thermoplastic resin known as JJS1361. Double Cantilever Beam Specimen Configuration and Test Methods As received from NASA, the test specimens were nominally 0.5 inch wide, 6 inches long, and 0.2 inch thick. A 1 inch long Kapton insert at the midplane of one end of the specimen (placed during laminate fabrication) facilitated crack initiation and extension. It was noted that the specimens provided were smaller than the nominal 1.5 inch wide, 9.0 inch long configuration specified. Similarly, the Kapton inserts were of greater length than those in the present specimens. Hence, the data below should not be compared directly to those generated with the referenced methods. No preconditioning was performed on the specimens prior to testing. In general, the methodology was used for the present work. Crack opening loads were introduced to the specimens via piano hinges attached to the main specimen faces at a single end of each specimen. Hinges were bolted to the specimens using the technique presented. The cracks were extended a small distance from the end of the Kapton insert prior to testing. Just before precracking, the sides of the specimens were coated with water-soluble typewriter correction fluid to aid in crack visualization. Scribe marks were then made in the coating at half-inch intervals.

Kessler, Jeff A.↗

DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗

DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Fitness For Service Assessment of a Corroded Heat Exchanger

Within the Fermi National Accelerator complex, there exist various water systems that support accelerator operations. One of these systems is extremely vital to the operation of the machine; that is the cooling system. The cooling system consists of nine relatively large heat exchangers that take untreated pond water and use it to cool the process fluid that further cools machine components. Over the 30 years these heat exchangers have been in operation, they have undergone significant material loss on the channels. This material loss, due to various forms of corrosion such as galvanic and microbiologically influenced corrosion (MIC) and possibly others, has deteriorated more than 80% of the nominal wall thickness of some of the exchangers and placed them in a questionable state. ASME FFS-1 (API 579) has been applied to address the condition of the heat exchangers due to their noncompliance with the governing code, BPVC Sec. VIII Div. 1. The assessments encompassed ASME FFS-1 parts 4: General Metal Loss and 9: Crack Like Flaw using level 1, 2, and 3 analysis techniques based on inspection data obtained by API 510 inspections. Level 1 and 2 assessments were deemed unfit for the corroded regions due to their location relative to a major structural discontinuity (channel to tube-sheet joint), so a level 3 analysis was conducted according to ASME Sec. VIII Div. 2 (design by analysis) rules for pressure vessels. Supplemental information included pond water tests to determine an accurate future corrosion allowance due to lacking inspection history. A leak before break (LBB) route was chosen to evaluate the possibility of leaking prior to the onset of failure. The analysis of one heat exchanger shows that the possibility the channel will develop a pinhole leak over 2.5 more years of operation should not be overlooked, but burst was unlikely from operation. The use of fracture mechanics show, that if a through-wall crack were to develop, it would not propagate further than the channel geometry and cause a leak not greater than 35 GPM. Using ASME Section XI Code Case N-705-1, allowing us to operate with a leak until the next outage given certain operating conditions and developing a leak mitigation procedure, this heat exchanger is deemed fit-for-service.

Humenik, Alex [Fermilab]↗

DOE OSTI2022

Characterization of Acoustic Emissions From Analogue Rocks Using Sparse Regression‐DMDc

Abstract Moisture loss in rock is known to generate acoustic emissions (AE). Phenomena that result in AE during drying are related to the movement of fluids through the pores and induced‐cracks that arise from differential mineral shrinkage, especially in clay‐bearing rock. AE from the movement of fluids occurs from the reconfiguration of fluid interfaces during drying, while AE from mineral shrinkage involves the debonding within or between minerals. Here, analogue rock samples were used to examine the differences in the AE signatures when one or both AE source‐types are present. An unsupervised sparse regression model, Dynamic Mode Decomposition with control, that extends Dynamic Mode Decomposition is used to characterize the AE signals recorded during the drying of porous analogue rock samples fabricated with ordinary Portland cement, with and without clay. This method can effectively and accurately reconstruct acoustic signals emitted from samples that only experience moisture loss without cracking. However, the method struggles to reconstruct signals from samples with intricate crack networks that formed during drying because AE generating mechanisms can emit contemporaneously, and the resulting waves propagate through drying‐induced cracks that can lead to multiple internal reflections. Thus, the differential reconstruction accuracy of time series generated by different underlying physical processes provides a robust filter for reducing large data catalogs. In general, both dynamics and sparse initiating events are learned directly from data and this method exposes a data hierarchy based on the complexity of the intrinsic dynamics.

58 GEOSCIENCES↗

DOE OSTI2022

Diatom influence on the production characteristics of hydrate-bearing sediments: Examples from Ulleung Basin, offshore South Korea

The Ulleung Basin Gas Hydrate field expeditions in 2007 (UBGH1) and 2010 (UBGH2) sought to assess the Basin's gas hydrate resource potential. Coring operations in both expeditions recovered evidence of gas hydrate, primarily as fracture-filling (or vein type) morphologies in mainly silt-sized, fine-grained sediment, but also as pore-occupying hydrate in the coarser-grained layers of interbedded sand and fine-grained systems. A commonality across many of these occurrences is the presence of diatoms in the fine-grained sediment. Here we tested fine-grained sediment (median grain size <12.5 μm) associated with hydrate occurrences at four UBGH2 sites (UBGH2-2-2, UBGH2-3, UBGH2-6 and UBGH2-11) to investigate potential impacts of diatoms on efforts to extract methane from hydrate, or to tap hydrocarbon reservoirs beneath hydrate-bearing sediment. Two key considerations are: the extent to which diatoms control sediment mechanical properties, and the extent to which pore-water freshening, which occurs as gas hydrate breaks down during resource extraction, alters the diatom control on sediment mechanical properties. We conducted experiments to measure sediment index properties, sedimentation behavior and compressibility to address these considerations. We relied on scanning electron microscope (SEM) imagery and X-ray powder diffraction (XRD) to characterize the sediment mineralogy. Our high-level findings are that at the ~20–45% (by volume) diatom concentrations observed at these UBGH2 sites, sediment compressibility increases with diatom content, but diatoms only appear to increase porosity and permeability at the highest diatom concentration (~45%). Our measurements suggest in situ compression indices of 0.35–0.55 and permeabilities on the order of 0.01milliDarcies (1 × 10 -17 m 2 ) can be anticipated at these sites. Importantly, these properties are not expected to vary significantly upon pore water freshening that accompanies gas hydrate dissociation during production.

58 GEOSCIENCES↗

DOE OSTI2025

Transportability of exogenous microbial community correlates with interwell connectivity in deep aquifers

Subsurface resource engineering operations often utilize continuous injection of externally-sourced water into geological reservoirs for formation pressure maintenance, resource recovery or energy/waste storage. Such injected water generally contains naturally occurring microbes. Little is known, however, about how the injectate microbes transport through geological media as a community, how such transportability is affected by injector-producer connectivity, and whether such knowledge can be utilized for flowpath characterization. In this study, we analyzed daily-to-weekly timeseries microbial community data from the injected- and produced-fluids of a ten-month flow test at a deep, well-characterized engineered aquifer. We found that the injectate microbial community was distinct from the indigenous community at the amplicon sequence variant (ASV) level, and that the transportability of injectate community towards a given producer, quantified by an “nASV-Overlap” metric we propose, had strong and significant positive correlation with known injector-producer connectivities at our site. This suggests that the better the connectivity, the higher the probability for more injectate species to flow through the interwell region and arrive at a producer. Because interwell connectivity is an important yet usually unknown parameter in subsurface resource engineering, such correlation in turn points to nASV-Overlap as a useful indicator of interwell connectivity for aquifer characterization and long-term monitoring. Based on our findings, an nASV-Overlap-based microbial tracing approach was developed for characterizing and monitoring the relative connectivities across multiple producers with a given injector. A side-by-side comparison between the new nASV-Overlap approach and traditional artificial tracer methods is presented, and their respective strengths and limitations are discussed.

Deep biosphere↗

DOE OSTI2022

EGS Collab Experiment 2: Continuous Broadband Seismic Waveform Data

Two broadband seismometers were installed on the 4100 level and recorded for the duration of EGS Collab Experiment #2. Inspired by published data from similar instruments installed in the Aspo Hard Rock Lab, these long-period instruments aimed to measure the tilting of the drift in response to the injection of fluid into the testbed. One instrument was installed underneath the wellheads in Site A (aka the "battery" alcove) and the other was installed along the east wall of the drift, south of Site B. Due to the feet of gravel (ballast) laid along the floor of the drift, we were unable to anchor the sensors directly to the rock. As a result, the coupling of the sensors to the experiment rock volume is likely poor. In addition, there are a number of noise sources that complicate the interpretation of the data. For example, sensor BBB is installed adjacent (within 3 ft) to the rail line that runs towards the Ross shaft. Trains (motors) run along this line almost daily and produce a large signal in these data. Careful extraction of periods of interest, as well as filtering for specific signals, is necessary. The sensors are Nanometrics Trillium Compact Posthole seismometers, sensitive down to 120 seconds period. They were installed as close to the drift wall and as deep as we could manually excavate (only about 1 ft or so). The holes were leveled with sand and the sensors were placed on a paver before backfilling with sand. The hole was then covered by a bucket filled with insulation to improve the sensor's isolation from daily temperature variations, which are minor but present due to drift ventilation from the surface. Data were recorded on Nanometrics Centaur digitizers at 100 Hz. The full response information is available in the StationXML file provided here, or by querying the sensors through the IRIS DMC (see links below). These instruments were provided free of charge through the IRIS PASSCAL instrument center. The network code is XP and the station codes are BBA and BBB. The waveform data can be queried through the IRIS FDSN server using any method the user likes. One convenient option is to use the Obspy python package: https://docs.obspy.org/packages/obspy.clients.fdsn.html

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗