Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subsurface characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Deep Vadose Zone Treatability Test for the Hanford Central Plateau: Soil Desiccation Pilot Test Results

Over decades of operation, the U.S. Department of Energy (DOE) and its predecessors have released nearly 2 trillion L (450 billion gal) of liquid into the vadose zone at the Hanford Site. Much of this liquid waste discharge into the vadose zone occurred in the Central Plateau, a 200 km2 (75 mi2 ) area that includes approximately 800 waste sites. Some of the inorganic and radionuclide contaminants in the deep vadose zone at the Hanford Site are at depths below the limit of direct exposure pathways, but may need to be remediated to protect groundwater (DOE 2008a; Dresel et al. 2011). The Tri-Party Agencies (DOE, U.S. Environmental Protection Agency, and Washington State Department of Ecology) established Milestone M-015-50, which directed DOE to submit a treatability test plan for remediation of Tc-99 and uranium in the deep vadose zone. These contaminants are mobile in the subsurface environment and have been detected at high concentrations deep in the vadose zone, and at some locations have reached groundwater. Testing technologies for remediating Tc-99 and uranium will also provide information relevant for remediating other contaminants in the vadose zone. The desiccation test described herein was conducted as an element of the test plan published in March 2008 to meet Milestone M-015-50 (DOE 2008a). This field test report was prepared and submitted in response to Milestone M-015-110D, which set a date for desiccation field test activities to be reported to the Tri-Party Agencies. The desiccation field test was conducted at the Hanford Site 200-BC-1 Operable Unit. This waste site contains 26 cribs and trenches that received about 110 million L (29 million gal) of liquid waste primarily in the mid-1950s. The waste contained about 410 curies of Tc-99 (Corbin et al. 2005). There is no evidence the contamination has reached groundwater, located about 100 m (330 ft) below ground surface (bgs) in this area. Initial characterization efforts indicated the Tc-99 inventory is located mostly at a depth in the vadose zone of between about 30 and 70 m (98 and 230 ft) bgs. However, transport model predictions have indicated the potential for this contamination to adversely impact groundwater in the future (Ward et al. 2004). The test was conducted in two steps to provide information about desiccation that is intended for use in subsequent feasibility studies for waste sites with inorganic and radionuclide contaminants in the deep vadose zone. The first step was field-scale test site characterization conducted as described in a characterization work plan (DOE 2008b). Results of the characterization effort have been previously reported in DOE (2010a) and Um et al. (2009). A field test plan (DOE 2010b) was prepared and used to guide the desiccation field testing effort. Laboratory and numerical modeling efforts (Truex et al. 2011; Ward et al. 2008; Oostrom et al. 2009, 2011, 2012, In Press) preceded and accompanied the field test and are incorporated herein as their results pertain to assessment of desiccation for future feasibility studies.

desiccation↗
DOE OSTI2015

Metaproteomics Identifies the Protein Machinery Involved in Metal and Radionuclide Reduction in Subsurface Microbiomes and Elucidates Mechanisms and U(VI) Reduction Immobilization

The overall goal for this funded project was to develop and exploit environmental metaproteomics tools to identify biomarkers for monitoring microbial activity affecting U speciation at U-contaminated sites, correlate metaproteomics profiles with geochemical parameters and U(VI) reduction activity (or lack thereof), elucidate mechanisms contributing to U(VI) reduction, and provide remediation project managers with additional information to make science-based site management decisions for achieving cleanup goals more efficiently. Although significant progress has been made in elucidating the microbiology contribution to metal and radionuclide reduction, the cellular components, pathway(s), and mechanisms involved in U trans-formation remain poorly understood. Recent advances in (meta)proteomics technology enable detailed studies of complex samples, including environmental samples, which differ between sites and even show considerable variability within the same site (e.g., the Oak Ridge IFRC site). Additionally, site-specific geochemical conditions affect microbial activity and function, suggesting generalized assessment and interpretations may not suffice. This research effort integrated current understanding of the microbiology and biochemistry of U(VI) reduction and capitalize on advances in proteomics technology made over the past few years. Field-related analyses used Oak Ridge IFRC field ground water samples from locations where slow-release substrate biostimulation has been implemented to accelerate in situ U(VI) reduction rates. Our overarching hypothesis was that the metabolic signature in environmental samples, as deciphered by the metaproteome measurements, would show a relationship with U(VI) reduction activity. Since metaproteomic and metagenomic characterizations were computationally challenging and time-consuming, we used a tiered approach that combines database mining, controlled laboratory studies, U(VI) reduction activity measurements, phylogenetic analyses, and gene expression studies to support the metaproteomics characterizations. Growth experiments of target microorganisms (Anaeromyxobacter, Shewanella, Geobacter) revealed tremendous respiratory versatility, as evidenced by the ability to utilize a range of electron donors (e.g. acetate, hydrogen, pyruvate, lactate, succinate, formate) and electron acceptors (e.g. nitrate, fumarate, halogenated phenols, ferric iron, nitrous oxide, etc.). In particular, the dissimilatory metabolic reduction of metals, including radionuclides, by target microorganisms spurred interest for in situ bioremediation of contaminated soils and sediments. Distinct c-type cytochrome expression patterns were observed in target microorganisms grown with the different electron acceptors. For each target microorganism, the core proteome covered almost all metabolic pathways represented by their corresponding pan-proteomes. Unique proteins were detected for each target microorganism, and their expression and possible functionalities were linked to specific growth conditions through proteomics measurements. Optimization of the proteomic tools included in-depth comprehensive metagenomic and metaproteomic analyses on a limited number of samples. The optimized metaproteomic analyses were then applied to Oak Ridge IFRC field samples from the slow-release substrate biostimulation. Metaproteomic analysis and pathway mapping results demonstrated the distinct effects of metal and non-metal growth conditions on the proteome expression. With these metaproteomic tools, we identified which previously hypothetical metabolic pathways were active during the analyzed time points of the slow release substrate biostimulation. Thus, we demonstrated the utility of these tools for site assessment, efficient implementation of bioremediation and long-term monitoring. This research of detailed protein analysis linked with metal reduction activity did (1) show that c-type cytochrome isoforms, previously associated with radionuclide reduction activity, are suitable biomarkers, (2) identify new biomarker targets for site assessment and bioremediation monitoring, and (3) provide new information about specific proteins and mechanisms involved in U(VI) reduction and immobilization. This expanded metagenomic and metaproteomic toolbox contributed to implementing science-driven site management with broad benefits to the DOE mission.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computed Tomography Scanning and Geophysical Measurements of Appalachian Basin Core from the Jones and Laughlin #1 Well, Beaver County, PA

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia, were used to characterize Appalachian Basin core from Beaver County, Pennsylvania. The primary impetus of this work is a collaboration between the U.S. Department of Energy (DOE) and the Pennsylvania Geological Survey to characterize and make publicly available core information from the Onondaga-Huntersville formations of the Appalachian Basin. This stratigraphic well and the core data produced in this report will aid in understanding the structural complexities of the Onondaga-Huntersville formations. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/jonesandlaughlin1well. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with imaging performed with the NETL medical CT scanner over the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for further study. Targeted higher resolution CT scanning of select sections was performed with NETL’s micro-CT scanner. The combination of methods used provides a multiscale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

58 GEOSCIENCES↗
DOE OSTI2020

Pump-and-Treat Impacts on Aquifer Conditions and the Development of System-Wide Performance Strategies - 20377

The pump-and-treat (P and T) system on the Hanford Site's Central Plateau currently treats multiple contaminants within the 200-ZP-1 operable unit, including nitrate through its biological treatment processes. Although it was initially designed to support the 200-ZP-1 Remedial Action Objectives, the P and T system has evolved in use to address contaminants from other regions of the aquifer as well. The P and T influents are biogeochemically distinct, treating groundwaters with elevated concentrations of cyanide, uranium, technetium-99, and nitrate. Given the varying geochemical signatures from current and potential future sources, operational conditions, and amendments considered for treatment processes, there is a need to identify the impacts that the P and T effluent has on the aquifer, and how the aquifer affects P and T performance. To this end, laboratory efforts have focused on characterizing P and T effluent and its biogeochemical effects on the sediments through a series repacked column experiments. The P and T effluent is rich in particulate matter. A systematic inspection of this particulate matter in the P and T effluent showed significant levels of manganese/iron oxides, calcite crystals, relatively large organic matter molecules, and biomass. Therefore, the column experiments were designed and conducted for four different treatment systems: artificial groundwater, filtered P and T effluent, non-filtered P and T effluent, and chlorinated P and T effluent. Different levels of microbiologically- and geochemically-induced changes in the columns occurred and modified the flow systems. For example, physical filtration occurred in the sediments where the particulate material had accumulated in the top portion of the column. Reactive transport modeling of the laboratory experiments is being carried out so that the reaction network can be extrapolated to the field scale. These simulations will provide predictive evaluations of system performance based on P and T operational changes and aquifer mass loading and capacity estimates. This information can be used to address combined subsurface and operational challenges, such as the need for frequent well redevelopment to maintain injection capacity. This information also provides the technical basis for decisions related to the P and T system operations that support remedy optimization efforts and short- and long-term remedy decisions related to multiple operable units on the Hanford Site Central Plateau. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

A conceptual geochemical model of the geothermal system at Surprise Valley, CA

Characterizing the geothermal system at Surprise Valley (SV), northeastern California, is important for determining the sustainability of the energy resource, and mitigating hazards associated with hydrothermal eruptions that last occurred in 1951. Previous geochemical studies of the area attempted to reconcile different hot spring compositions on the western and eastern sides of the valley using scenarios of dilution, equilibration at low temperatures, surface evaporation, and differences in rock type along flow paths. These models were primarily supported using classical geothermometry methods, and generally assumed that fluids in the Lake City mud volcano area on the western side of the valley best reflect the composition of a deep geothermal fluid. In this contribution, we address controls on hot spring compositions using a different suite of geochemical tools, including optimized multicomponent geochemistry (GeoT) models, hot spring fluid major and trace element measurements, mineralogical observations, and stable isotope measurements of hot spring fluids and precipitated carbonates. We synthesize the results into a conceptual geochemical model of the Surprise Valley geothermal system, and show that high-temperature (quartz, Na/K, Na/K/Ca) classical geothermometers fail to predict maximum subsurface temperatures because fluids re-equilibrated at progressively lower temperatures during outflow, including in the Lake City area. We propose a model where hot spring fluids originate as a mixture between a deep thermal brine and modern meteoric fluids, with a seasonally variable mixing ratio. The deep brine has deuterium values at least 3 to 4‰ lighter than any known groundwater or high-elevation snow previously measured in and adjacent to SV, suggesting it was recharged during the Pleistocene when meteoric fluids had lower deuterium values. The deuterium values and compositional characteristics of the deep brine have only been identified in thermal springs and groundwater samples collected in proximity to structures that transmit thermal fluids, suggesting the brine may be thermal in nature. On the western side of the valley at the Lake City mud volcano, the deep brine-meteoric water mixture subsequently boils in the shallow subsurface, precipitates calcite, and re-equilibrates at about 130 °C. On the eastern side of the valley, meteoric fluid mixes to a greater extent with the deep brine, cools conductively without boiling, and the composition is modified as dissolved elements are sequestered by secondary minerals that form along the cooling and outflow path at temperatures <130 °C. Re-equilibration of geothermal fluids at lower temperatures during outflow explains why subsurface temperature estimates based on classical geothermometry methods are highly variable, and fail to agree with temperature estimates based on dissolved sulfate-oxygen isotopes and results of classical and multicomponent geothermometry applied to reconstructed deep well fluids. Furthermore, the proposed model is compatible with the idea suggested by others that thermal fluids on the western and eastern side of the valley have a common source, and supports the hypothesis that low temperature re-equilibration during west to east flow is the major control on hot spring fluid compositions, rather than dilution, evaporation, or differences in rock type.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

Groundwater Remediation in a Flood plain Aquifer at Shiprock, New Mexico - 16097

The Shiprock, New Mexico, Disposal Site is the location of a former uranium and vanadium ore-processing mill that operated from 1954 to 1968. During milling years, tailings leachate and other mill-related liquids migrated into underlying sediments, resulting in the contamination of groundwater beneath the river terrace containing the mill and a nearby flood plain alluvial aquifer adjacent to the San Juan River. The US DOE Office of Legacy Management manages the site and is implementing the current groundwater remedy for the flood plain area, consisting of natural flushing supplemented by extraction of alluvial groundwater to accelerate reduction in contaminant concentrations. The alluvial aquifer contains elevated concentrations of sulfate, uranium, and nitrate derived largely from the historical eastward migration of groundwater through terrace alluvium and underlying bedrock to the flood plain. The bedrock formation beneath the entire Shiprock site is Mancos shale, which forms an escarpment that separates the terrace from the flood plain. Contamination in the flood plain alluvial aquifer is influenced by a dynamic flow system that changes seasonally and from year to year. Simulations with a flood plain flow model based on monitored groundwater levels and coincident river-surface elevations show that background groundwater flow, unaffected by remediation pumping, is characterized by distinct hydraulic conditions for three separate periods each year. During months of low river flow in winter and early spring, river losses to the aquifer create a large hypo-rheic zone in the south half of the flood plain. Simultaneously, recharge of continually flowing surface water emptying onto the north half of the flood plain from Bob Lee Wash creates divergent, radial flow from the wash outlet and diverts northwestward-migrating contaminant plumes to the northeast where they discharge to the river. Subsequent bank storage processes in May and June from high snow melt runoff change background flow directions, increasing the spread of contamination in the alluvium. From mid-summer through fall, river flows are typically similar in value to those observed prior to the onset of snow melt, but groundwater elevations in the aquifer are noticeably lower than water levels observed in winter and early spring and may induce river losses to the aquifer. Several different contaminant sources have impacted the alluvial aquifer. Much of the contamination originated as tailings leachate and raffinate wastewater that infiltrated the terrace subsurface and eventually reached the flood plain at the escarpment. Additional sources of contamination in the aquifer included discharge of mill effluent to Bob Lee Wash and a pond at the base of the escarpment, and aqueous mobilization of solid-phase contamination in the form of windblown deposits and contaminants adsorbed to alluvial sediments. The degree to which remnant, secondary-source processes, such as leaching of solid-phase deposits and discharge of terrace groundwater across the escarpment, contribute to existing contaminant plumes is uncertain. The flood plain remediation system consists of two groundwater extraction wells, pumping from two horizontal wells placed in trenches near the base of the escarpment, a sump collecting discharges from seeps in the escarpment wall, and discharge of collected water to an evaporation pond on the terrace. Flow-model simulations indicate that most contaminated areas are captured by remediation pumping, resulting in several zones of freshwater between the river and the escarpment. Despite the highly transient nature of groundwater flow beneath the flood plain and the varied contaminant sources, initiation of remediation pumping in 2003 and coordinated operation of the system's various components have successfully removed large amounts of contaminant mass and helped prevent discharge of contaminated groundwater to the river. The percentage reductions of dissolved sulfate, uranium, and nitrate mass in the aquifer estimated for the first 9 years of remediation range from 40% to 72%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

H-2-SALT: Storing Fossil Energy as Hydrogen in Salt Caverns

Hydrogen storage in subsurface caverns is known as perhaps largest and longest-term energy storage system. While few such systems are currently in service (e.g., Texas; Teesside, UK), they are envisioned as being a pillar of the current energy transition that will allow intermittent power production from renewable sources (e.g., wind, solar) to be balanced with demand at grid-scale. In addition, hydrogen energy storage can allow existing fossil power plants to run more economically either by minimizing startups/shutdowns or by allowing them to take advantage of arbitrage. Finally, stored hydrogen can feed a variety of non-power users such as pipeline gas, heating, transportation, and manufacturing. The H-2-SALT paper study assessed the feasibility of a power-to-hydrogen system utilizing salt cavern storage of hydrogen in bedded salt in central Kansas, where over 750 such caverns have been constructed to date, of which approximately 350 are still in service. Legacy well and salt cavern data were collected to develop a regional geological database for two sites in Central Kansas and adjacent areas for use in more detailed geological and geomechanical site characterization efforts. This study found commercial viability and competitiveness of large-scale, electrolytic H 2 production and storage ($1.78 per kg H 2 ) that can be used for both electrical power supply to the grid during high-priced electricity periods (including co-firing with NG and fuel cell power production) and sale of H 2 for various industries, such as petrochemicals or transportation. The commercial potential of the H-2-SALT system is also relatively high because each of its components operate commercially today. It is recommended that a stratigraphic test well be drilled to collect log data and core samples that can reduce uncertainty surrounding geological and geomechanical properties of the salt beds that will be required for further cavern and energy storage system design studies.

03 NATURAL GAS↗
DOE OSTI2026

Development Fiber Optic Distributed System for Direct Detection of Subsurface Gases Leakages

Carbon, natural gas, and hydrogen gas storage is an emerging solution to safeguard us against pollution, support goals of negative carbon emission, and protect sources of renewable energy. Properly constructed storage wells provide a virtually impervious barrier to any unintended subsurface transmission. The ability to ensure the long-term integrity of such wells is vital to the success of any storage operation and be successful in the public eyes. Therefore, robust monitoring of any gas migration into the subsurface is highly sought. A fiber-optic distributed chemical sensor (DCS) enables monitoring of long-term well integrity along its depth, ensuring the success of any storage operation and bolsters public acceptance of the safety of the reservoir via leak early detection. The same technique can be applied to gas monitoring in pipeline networks and nuclear stockpile monitoring applications. Fiber based Raman spectroscopy enables DCS, as optical fibers can be deployed in virtually any environment and relay spectroscopic information over long distances back to the user. Hollow core fibers (HCF) make excellent DCSs as the air core of the fiber allows gas from the environment to diffuse into the core, which interacts with the laser signal that is carried in the air core. This work builds upon the previous LDRD project, Fiber Optic System for Direct Detection of Carbon Dioxide Leakage in Carbon Storage Wells (21-FS-003), in which the feasibility of Raman spectroscopy detection of Carbon Dioxide (CO2) in HCF detection was demonstrated. We mitigated the risk of this DCS technology by establishing and completing five objectives. The first objective was to model and optically characterize HCF uptake of CO2, establishing the relationship between HCF length, gas diffusion time, detectable gas concentration, and measured Raman intensity. In objective two, we developed a fiber core drilling recipe to enable additional diffusion ports in the fiber core and established a method for maintaining fiber strength and integrity post drilling. Objective three characterized the drilled fibers against the undrilled fibers, establishing the differences in the gas mechanics and optical properties and provided parameters to iterate the drilling process. In objective four, a fusion splicing technique was developed to join the HCF to conventional single-mode fibers, localizing the gas detection point at the drilled HCF hole, emulating a DCS. Lastly, objective five was the testing of the sensor in Edgar Mines at Colorado School of Mines on a CO2 pipeline with a simulated leak, to showcase the ability to detect CO2 leaks. This capstone result showed CO2 leak detection in < 10 minutes, raising the technology readiness level of HCF segments as deployable DCS.

organic↗
DOE OSTI2018

Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO 2 and N 2

Rivers in climatic zones characterized by dry and wet seasons often experience periodic transitions between losing and gaining conditions across the river-aquifer continuum. Infiltration shifts can stimulate hyporheic microbial biomass growth and cycling of riverine carbon and nitrogen leading to major exports of biogenic CO 2 and N 2 to rivers. In this study, we develop and test a numerical model that simulates biological-physical feedback in the hyporheic zone. We used the model to explore different initial conditions in terms of dissolved organic carbon availability, sediment characteristics, and stochastic variability in aerobic and anaerobic conditions from water table fluctuations. Our results show that while highly losing rivers have greater hyporheic CO 2 and N 2 production, gaining rivers allowed the greatest fraction of CO 2 and N 2 production to return to the river. Hyporheic aerobic respiration and denitrification contributed 0.1–2 g/m 2 /d of CO 2 and 0.01–0.2 g/m 2 /d of N 2 ; however, the suite of potential microbial behaviors varied greatly among sediment characteristics. We found that losing rivers that consistently lacked an exit pathway can store up to 100% of the entering C/N as subsurface biomass and dissolved gas. Finally, our results demonstrate the importance of subsurface feedbacks whereby microbes and hydrology jointly control fate of C and N and are strongly linked to wet-season control of initial sediment conditions and hydrologic control of seepage direction. These results provide a new understanding of hydrobiological and sediment-based controls on hyporheic zone respiration, including a new explanation for the occurrence of anoxic microzones and large denitrification rates in gravelly riverbeds.

58 GEOSCIENCES↗
DOE OSTI2015

Dynamics of Magnesite Formation at Low-Temperature and High pCO2 in Aqueous Solution

Like many metal carbonate minerals, despite conditions of supersaturation, precipitation of magnesite from aqueous solution is kinetically hindered at low temperatures, for reasons that remain poorly understood. The present study examines precipitation products from reaction of Mg(OH)2 in aqueous solutions saturated with supercritical CO2 at high pressures (90 atm and 110 atm) and low temperatures (35 °C and 50 °C). Traditional bulk characterization (X-ray diffraction) of the initial solid formed indicated the presence of hydrated magnesium carbonates (hydromagnesite and nesquehonite), thermodynamically metastable phases that were found to slowly react during ageing to the more stable anhydrous form, magnesite, at temperatures as low as 35 °C (135-140 days) and at a faster rate at 50 °C (56 days). Undetected by bulk measurements, detailed examination of the precipitates by scanning electron microscopy (SEM) showed that magnesite is present as a minor component at relatively early reaction times (7 days) at 50 °C. In addition to magnesite dominating the solid phases over time, we find that mangesite nucleation and growth occurs more quickly with increasing partial pressure of CO2, and in electrolyte solutions with high bicarbonate content. Furthermore, formation of magnesite was found to be enhanced in sulfate-rich solutions, compared to chloride-rich solutions. We speculate that much of this behavior is possibly due to sulfate serving as sink of protons generated during carbonation reactions. These results support the importance of integrating magnesite as an equilibrium phase in reactive transport calculations of the effects of carbon dioxide sequestration on subsurface formations at long time scales.

Qafoku, Odeta↗
DOE OSTI2020

Temporal flow variations interact with spatial physical heterogeneity to impact solute transport in managed river corridors

The interactions between surface water and groundwater in river corridors lead to temporal fluctuations in subsurface water fluxes which have a critical role on solute transport dynamics. Here, we develop a framework to analyze the relative impacts of different temporal frequencies of the flow field in a spatially heterogeneous aquifer on solute transport. Our analysis indicates that the advection-dispersion equation behaves as a low-pass filter by wiping out the effect of high-frequency velocity fluctuations on the first two spatial moments of the solute plume, namely its center of mass and spreading. The concepts discussed in the theoretical analysis are then applied to understand solute transport dynamics at the 300 Area of the Hanford site (USA) adjacent to the Columbia River. We examine the temporal behavior of the solute plume's spatial moments for different temporal frequencies utilizing geostatistical parameters estimated in the 300 Area. Due to the proximity to the Columbia river, groundwater fluxes at the Hanford site are highly dynamic resulting in a large range of characteristic temporal frequencies. Nonetheless, similar to the theoretical analysis, our results show that the effect of high-frequency fluctuations is filtered, with most of the solute transport dynamics being controlled by fluctuations characterized by a large characteristic period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Area V: Waste Management and Subsurface Science (FY21-23 Final Report)

The objective of the Waste Management and Subsurface Science Topic Area V project is to evaluate the feasibility and safety of an intermediate-depth borehole for disposal of Israel Atomic Energy Commission (IAEC) nuclear waste. The IAEC nuclear waste is currently stored in above ground and/or in near-surface facilities that provide short-term radiological worker safety but do not provide long-term, full security and safeguards protection for the nuclear material. Disposal of the IAEC nuclear waste several hundred meters below the surface in a sealed borehole would provide long-term security and safeguards due to its depth underground, as well as long-term radiological safety, due to its isolation in a stable geological formation. This project supports the nuclear security and safeguards initiatives of the Science and Technology Working Group of the NNSA/IAEC Memorandum of Understanding (MOU) and addresses the Topic Area V objective and associated goals. The goals for Topic Area V are: (1) Assess the suitability of the vadose zone in the Yamin Plain of the northern Negev Desert, Israel as a host location for intermediate-depth borehole disposal of radioactive waste, (2) Define and refine the necessary characterization needs to better understand the regional climate, structure, and stratigraphy, geomechanics, geochemistry, hydrogeology, and seismicity of the stratigraphic sequence from the surface to the water table in the vicinity of the proposed host location for the disposal borehole(s), and (3) Develop a preliminary plan and schedule, including necessary resources and costs, for a borehole disposal concept and robust and transparent pre-closure and post-closure safety assessments. The safety assessments will require iteration with process level models at both site and regional scales. Upon completion, the information collected by this effort will allow the IAEC to make decisions regarding the viability and applicability of an intermediate depth (<500m) borehole for disposal in the vadose zone of the Yamin Plain. This FY21-23 final report is provided in lieu of a workshop report that was originally planned for delivery in January, 2024. Due to the conflict in Israel and associated delays, the workshop has been delayed to April 2024. The workshop will be held at LLNL but offsite and will switch to hybrid if travel restrictions remain in place. The associated workshop report will be provided as part of a new FY24-26 LCP.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Evaluating Impacts of CO2 Gas Intrusion Into a Confined Sandstone aquifer: Experimental Results

Deep subsurface storage and sequestration of CO2 has been identified as a potential mitigation technique for rising atmospheric CO2 concentrations. Sequestered CO2 represents a potential risk to overlying aquifers if the CO2 leaks from the deep storage reservoir. Experimental and modeling work is required to evaluate potential risks to groundwater quality and develop a systematic understanding of how CO2 leakage may cause important changes in aquifer chemistry and mineralogy by promoting dissolution/precipitation, adsorption/desorption, and redox reactions. Sediments from the High Plains aquifer in Kansas, United States, were used in this investigation, which is part of the National Risk Assessment Partnership Program sponsored by the US Department of Energy. This aquifer was selected to be representative of consolidated sand and gravel/sandstone aquifers overlying potential CO2 sequestration repositories within the continental US. In this paper, we present results from batch experiments conducted at room temperature and atmospheric pressure with four High Plains aquifer sediments. Batch experiments simulate sudden, fast, and short-lived releases of the CO2 gas as would occur in the case of well failure during injection. Time-dependent release of major, minor, and trace elements were determined by analyzing the contacting solutions. Characterization studies demonstrated that the High Plains aquifer sediments were abundant in quartz and feldspars, and contained about 15 to 20 wt% montmorillonite and up to 5 wt% micas. Some of the High Plains aquifer sediments contained no calcite, while others had up to about 7 wt% calcite. The strong acid extraction tests confirmed that in addition to the usual elements present in most soils, rocks, and sediments, the High Plains aquifer sediments had appreciable amounts of As, Cd, Pb, Cu, and occasionally Zn, which potentially may be mobilized from the solid to the aqueous phase during or after exposure to CO2. However, the results from the batch experiments showed that the High Plains sediments mobilized only low concentrations of trace elements (potential contaminants), which were detected occasionally in the aqueous phase during these experiments. Importantly, these occurrences were more frequent in the calcite-free sediment. Results from these investigations provide useful information to support site selection, risk assessment, and public education efforts associated with geological CO2 storage and sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Transmission electron microscopy characterization of the interfacial structure of a galvanized dual-phase steel

Site-specific studies were carried out to characterize the interface of a galvanized dual-phase (DP) steel. Focused ion beam (FIB) was used to prepare specimens in the interface region (~ 100 nm thick) between the coating and the substrate. Transmission electron microscopy (TEM), scanning TEM (STEM), and high resolution TEM (HRTEM) were performed to resolve the phases and the structures at the interface between the zinc (Zn) coating and the steel substrate. The STEM and TEM results showed that a continuous manganese oxide (MnO) film with a thickness of ~ 20 nm was present on the surface of the substrate while no silicon (Si) oxides were resolved. Internal oxide particles were observed as well in the sub-surface region. Despite the presence of the continuous oxide film, a well-developed inhibition layer was observed right on top of the oxide film. The inhibition layer has a thickness of ~ 100 nm. Possible mechanisms for the growth of the inhibition layer were discussed. - Highlights: •Site-specific examinations were performed on the Zn/steel interface. •Continuous external MnO oxides (20 nm) were observed at the interface. •No Si oxides were observed at the interface. •Internal oxide particles were distributed in the subsurface. •A continuous inhibition layer grew on top of the external oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling Frequency–Independent Q Viscoacoustic Wave Propagation in Heterogeneous Media

Quantifying the attenuation of seismic waves propagating in the Earth interior is critical to study the subsurface structure. Previous studies have proposed fractional anelastic wave equations to model the frequency–independent Q seismic wave propagation. Such wave equations involve fractional derivatives that pose computational challenges for the numerical schemes in terms of accuracy and efficiency when dealing with heterogeneous Earth media. To tackle these challenges, here we derive a new viscoacoustic wave equation, where the power terms of the fractional Laplacian operators are spatially independent, thus accurate and efficient methods (e.g., the Fourier pseudospectral method) can be adopted. Our derivation enables the resultant equation to capture both amplitude and phase signatures of the anelastic wave propagation by matching the complex wave numbers for all the frequencies of interest. We verify the derivation by comparing the dispersion curves of both the attenuation factor and the phase velocity produced by the new wave equation with their theoretical values as well as the Pierre Shale in situ measurements. Following that, we use a synthetic attenuating gas chimney model to demonstrate the attenuation effects on seismic waveforms and then construct the Q–compensated reverse time migration to undo these effects for seismic image enhancement. Finally, we find that our forward modeling results can characterize the spatiotemporal attenuation effects revealed in the Frio–II CO2 injection time–lapse seismic monitoring data. Here, we expect this proposed equation to be useful to quantify the attenuation in seismic data to push the resolution limits of seismic imaging and inversion.

25 ENERGY STORAGE↗
DOE OSTI2018

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development

This DOE GTO funded project was a collaborative effort between Array Information Technology and the University of California at Berkeley. The project goal is to improve technology to assess in-situ stress magnitude and orientation, kinematic fracture parameters, rupture size, rupture orientation as well as temporal changes and volumetric distribution of the injected fluid during Enhanced Geothermal System (EGS) resource development. The project leveraged high-frequency seismic data recorded by the LBNL 34-station permanent geophone network and seismic broadband data recorded by a temporary 33-station seismometer network that operated in The Geysers during the injection phase of the DOE GTO funded EGS demonstration project at the well Prati-32. The operation of the broadband network with high station density and good azimuthal coverage to monitor the injection phase produced an unprecedented dataset that is typically unavailable for EGS operations. The two datasets offer the advantages of studying the parameters under investigation over a broader frequency band and the higher sensitivity of the broadband seismometers resulted in a richer dataset with a lower magnitude of completeness. Moment tensor (MT) inversion have been employed to develop a catalog of earthquake source solutions including moment tensor and finite source rupture area. This resulted in a 168-event waveform-based seismic moment tensor catalog for events 0.7≤Mw≤3.9 based on a semi-automated moment tensor analysis technique. The MT solutions were used for in-situ stress estimation during the EGS injection phases. We found an approximate 15-degree counterclockwise rotation of the least compressive stress sigma-3, and a rotation of the maximum compressive stress sigma-1 toward the vertical as the injected volume of water increased. We developed a rupture-area magnitude scaling relationship of The Geysers earthquakes obtained from finite-source inversion for fault slip that reveals a very high correlation to published results. Based on the rupture area-magnitude relationship we mapped Mw from the LBNL earthquake catalog to estimate the activated fracture area, its location and orientation within the injection volume of Prati-32. The double-difference Wadati (DDW) technique was applied to high-accuracy P- and S-wave differential travel times derived from waveform cross correlation to image fluid saturation in the subsurface based on Vp/Vs-ratio. The recently developed technique has been applied to estimate temporal changes of fluid distribution in fault zones and in volcanic regions but not yet in engineered geo-reservoirs. The improvement of the technology includes high-resolution waveform cross correlation applied to the high-quality datasets with lower magnitude of completeness resulting in a higher number of recorded events. The improvement allowed to estimate both, the spatial and temporal distribution of injected fluids and the state of the injected fluid during reservoir stimulation. Results of the DDW analysis revealed the development of an injection-derived steam plume during the early injection phase, particularly in the northern part of the Prati-32 study area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A dynamic Bayesian optimized active recommender system for curiosity-driven partially Human-in-the-loop automated experiments

Optimization of experimental materials synthesis and characterization through active learning methods has been growing over the last decade, with examples ranging from measurements of diffraction on combinatorial alloys at synchrotrons, to searches through chemical space with automated synthesis robots for perovskites. In virtually all cases, the target property of interest for optimization is defined a priori with the ability to shift the trajectory of the optimization based on human-identified findings during the experiment is lacking. Thus, to highlight the best of both human operators and AI-driven experiments, here we present the development of a human–AI collaborated experimental workflow, via a Bayesian optimized active recommender system (BOARS), to shape targets on the fly with human real-time feedback. Here, the human guidance overpowers AI at early iteration when prior knowledge (uncertainty) is minimal (higher), while the AI overpowers the human during later iterations to accelerate the process with the human-assessed goal. We showcase examples of this framework applied to pre-acquired piezoresponse force spectroscopy of a ferroelectric thin film, and in real-time on an atomic force microscope, with human assessment to find symmetric hysteresis loops. It is found that such features appear more affected by subsurface defects than the local domain structure. This work shows the utility of human–AI approaches for curiosity driven exploration of systems across experimental domains.

36 MATERIALS SCIENCE↗