Search NASA⌕ Search

SEARCH · Search NASA

Results for “depth migration”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Modeling magnetic fields from a DC power cable buried beneath San Francisco Bay based on empirical measurements

Here, the Trans Bay Cable (TBC) is a ±200-kilovolt (kV), 400 MW 85-km long High Voltage Direct Current (DC) buried transmission line linking Pittsburg, CA with San Francisco, CA (SF) beneath the San Francisco Estuary. The TBC runs parallel to the migratory route of various marine species, including green sturgeon, Chinook salmon, and steelhead trout. In July and August 2014, an extensive series of magnetic field measurements were taken using a pair of submerged Geometrics magnetometers towed behind a survey vessel in four locations in the San Francisco estuary along profiles that cross the cable’s path; these included the San Francisco-Oakland Bay Bridge (BB), the Richmond-San Rafael Bridge (RSR), the Benicia- Martinez Bridge (Ben) and an area in San Pablo Bay (SP) in which a bridge is not present. In this paper, we apply basic formulas that ideally describe the magnetic field from a DC cable summed vectorially with the background geomagnetic field (in the absence of other sources that would perturb the ambient field) to derive characteristics of the cable that are otherwise not immediately observable. Magnetic field profiles from measurements taken along 170 survey lines were inspected visually for evidence of a distinct pattern representing the presence of the cable. Many profiles were dominated by field distortions unrelated to the cable caused by bridge structures or other submerged objects, and the cable’s contribution to the field was not detectable. BB, with 40 of the survey lines, did not yield usable data for these reasons. The unrelated anomalies could be up to 100 times greater than those from the cable. In total, discernible magnetic field profiles measured from 76 survey lines were regressed against the equations, representing eight days of measurement. The modeled field anomalies due to the cable (the difference between the maximum and minimum field along the survey line at the cable crossing) were virtually identical to the measured values. The modeling yielded a pooled cable depth below the bay floor of 2.06 m (±1.46 std dev), and estimated the angle to the horizontal of the imaginary line connecting the crosssectional center of the cable’s two conductors (0.1143 m apart) as 178.9° ±61.9° (std dev) for Ben, 78.6°±37.0° (std dev) for RSR, and 139.9°±27.4° (std dev) for SP. The mean of the eight daily average currents derived from the regressions was 986 ±185 amperes (A) (std dev), as compared to 722 ±95 A (std dev) provided by Trans Bay Cable LLC. Overall, the regressions based on fundamental principles (Biot Savart law) and the vectorial summation of cable and geomagnetic fields provide estimates of cable characteristics consistent with plausible expectations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deformation band corrosion as a result of applied stress and cold-work in 304 stainless steel in simulated pressurized water reactor primary water conditions

Stainless steels like 304SS are commonly used in pressurized water reactor primary water environments due to their excellent resistance to stress corrosion cracking (SCC). However, SCC susceptibility has been observed in cold-worked 304SS. Cold-rolling introduces a high density of defects, such as deformation bands or dislocations, into the microstructure; applied stress, even if macroscopically below the yield stress, can result in regions with high localized stresses. Cold-worked and non-cold-worked samples were exposed to 360 °C simulated PWR primary water for 1500 h with and without active loading. Corrosion susceptibility of dislocations and deformation bands induced by cold-rolling and/or applied stress was studied using high-resolution (analytical) transmission electron microscopy. Deformation bands were generated in all specimens examined in the study. Penetrative oxidation was observed in some deformation bands in each specimen. The depth of penetrative deformation band oxidation frequently exceeded the depth of grain boundary oxidation in all specimens. However, deformation band oxides generally were narrowly confined to the deformation band plane. Highly Ni-rich regions were observed beyond all oxidation fronts. Surface oxides were thicker in the specimens exposed under active loading. In conclusion, the implications of these findings at surfaces for deformation bands intersecting stress corrosion crack paths are discussed.

EFTEM↗
NASA NTRS2002

Distinguishing Indigenous from Contaminating Microorganisms in Rock Samples from a Deep Au Mine in South Africa

The concentration and distribution of microbial biomass within deep subsurface rock strata is not well known To date, most analyses are from water samples and a few cores. Hand samples, block samples and cores from an actively mined Carbon Leader ore zone at 3.2 kilometers depth were collected for microbial analyses. The Carbon Leader was comprised of quartz, S-bearing aromatic hydrocarbons, Fe(III) oxyhydroxides, sulfides, uraninite, Au and minor amounts of sulfate. The porosity of the ore was 1% and the maximum pore throat diameter was less than 0.1 microns; whereas, the porosity of the adjacent quartzite was .02 to .9% with a maximum pore throat diameter of 0.9 microns. Rhodamine dye, fluorescent microspheres, microbial enrichments, autoradiography, phospholipid fatty acid (PLEA) and 16S rDNA analyses were performed on these rock samples and the mining water. The date indicate that the levels of solute contamination less than 0.01% for pared rock samples. Despite this low level of contamination, PLEA, microbial enrichment, DNA and tracer analyses and calculations indicate that most of the viable microorganisms in the Carbon Leader represent gram negative aerobic heterotrophs and ammonia oxidizers that are phylogenetically identical or closely related to service water microorganisms. These microbial contaminants probably infiltrated the low permeability rock through mining-induced microfractures. Geochemical data also detected drilling water in a fault zone approx. 1 meter behind the rock face encountered during coring. The mining induced macrofractures that are common at these great depths act as pathways for the drilling water borne microorganisms into the lower temperature zone that extends several meters into rock strata from the rock face. Combined PLEA and T- RFLP analyses of the service water and Carbon Leader samples indicate that the concentration of indigenous microorganisms was less than 10(exp 2) cells/gram. Such a low concentrations result from the submicron pore throat diameters. PLFA. SO4-35 autoradiography and tracer analyses indicate that the bounding quartzite contains thermophilic sulfate reducing bacteria at 10(exp 3) cells/gram that are not attributable to drilling water contamination. The microorganisms may be surviving on sulfate generated by oxidation of sulfide by radiolytic reactions resulting from the high U concentration in the ore zone. The presence of up to 8,000 ppm of Fe(III) oxyhydroxides in the host rock will also act to recycle sulfide generated by the sulfate reducing bacteria into sulfate. The activity of these sulfate-reducing bacteria may be enhanced by mining induced fracturing which can propagate up to 40 meters into virgin rock where the temperatures are ca. 50 C, and decrepitate of sulfate rich fluid inclusions. In ultra deep mines, judicious application of tracers and multiple microbial characterization techniques can distinguish microbial contamination caused by the near field fracturing and drilling water migration from the indigenous microbial communities in rock strata. The importance of far field fracturing on indigenous microbial communities, however, remains unknown.

Onstott, T. C.↗
DOE OSTI2020

Safety Functions and Features, Events and Processes for the E-Area Performance Assessment

The DOE Technical Standard, “Disposal Authorization Statement and Tank Closure Documentation,” (DOE 2017) recommends the use of safety functions and features, events and processes (FEPs) to support development of conceptual models and identification of scenarios to be considered in a performance assessment (PA). The FEP process provides a means to describe how a PA considers and addresses the factors that could influence the performance of key barriers. Understanding the roles of barriers in terms of limiting migration helps to focus on how changes in the system could lead to a situation where those roles cannot be fulfilled and there is the potential for compromised performance. The FEPs screening and review process was used to identify FEPs that are relevant for the EArea Low-Level Waste Facility (LLWF) and specifically those FEPs that could have a detrimental impact on the effectiveness of a given safety function. For this PA, a default list of FEPs developed at the International Atomic Energy Agency (IAEA 2004) and an approach implemented for PAs at the Hanford and Idaho sites (Mehta et al. 2016, DOE-ID 2019) are used to identify processes and events that could influence the effectiveness of a given safety function for the E-Area LLWF (e.g., subsidence can impact the safety function of the cover system and lead to increased infiltration). The Hanford and Idaho PAs represent two of the most recent applications of this approach. The PA evaluates the potential impacts of changes in performance of different features of the system and demonstrates that the safety functions represent multiple and redundant barriers. Barrier analyses, assuming a safety function is not present, also test the robustness of the system in the event of the loss of one or more safety functions. Such evaluations also support a qualitative illustration of the concept of defense in depth. The safety concept for closure of the E-Area LLWF (generically referred to as “E-Area”) encompasses a variety of different features (i.e., administrative controls, natural site features, and engineered barriers) that reduce the potential impacts on human health and the environment from the residual waste that will remain after closure. These features can be represented as a collection of safety functions acting independently and as a system to provide for overall safety. In some applications, there have been attempts to assign numerical expectations to specific safety functions, but that is not the intent in this case. The concept of safety functions is used more qualitatively in two ways for this PA: 1. To illustrate the robustness of the E-Area design, operational practices and closure approach by documenting features that are and are not credited in different modeling cases. 2. To identify the roles of the different features and potential processes and events that could compromise the performance of safety features and need to be considered when developing the modeling approach. This report addresses both safety functions and FEPs for the E-Area PA.

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

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗
DOE OSTI2016

Order-of-magnitude differences in retention of low-energy Ar implanted in Si and SiO{sub 2}

The retention of 1 and 5 keV Ar implanted at 45° in Si and 4.3 nm SiO{sub 2} on Si was studied at fluences between 3 × 10{sup 14} and 1.5 × 10{sup 16} cm{sup −2}. X-ray photoelectron spectroscopy (XPS) served to monitor the accumulation of Ar as well as the removal of SiO{sub 2}. Bombardment induced changes in oxygen chemistry caused the O 1s peak position to move toward lower binding energies by as much as 2.2 eV. Plotted versus depth of erosion, the fluence dependent changes in oxygen content, and peak position were similar at 1 and 5 keV. The Ar content of Si increased with increasing exposure, saturating at fluences of ∼2 × 10{sup 15} cm{sup −2} (1 keV) and ∼6 × 10{sup 15} cm{sup −2} (5 keV). Much less Ar was retained in the SiO{sub 2}/Si sample, notably at 1 keV, in which case the low-fluence Ar signal amounted to only 8% of the Si reference. The results imply that essentially no Ar was trapped in undamaged SiO{sub 2}, i.e., the Ar atoms initially observed by XPS were located underneath the oxide. At the lowest fluence of 5 keV Ar, the retention ratio was much higher (43%) because the oxide was already highly damaged, with an associated loss of oxygen. The interpretation was assisted by TRIM(SRIM) calculations of damage production. Partial maloperation of the ion beam raster unit, identified only at a late stage of this work, enforced a study on the uniformity of bombardment. The desired information could be obtained by determining x,y line scan profiles of O 1s across partially eroded SiO{sub 2}/Si samples. Fluence dependent Ar retention in Si was described using an extended version of the rapid relocation model which takes into account that insoluble implanted rare-gas atoms tend to migrate to the surface readily under ongoing bombardment. The range parameters required for the modeling were determined using TRIM(SRIM); sputtering yields were derived from the literature. The other three parameters determining the Ar signal, i.e., (1) the thickness w of the near-surface Si region devoid of Ar, (2) the relocation efficiency Ψ{sub rlc}, and (3) the effective attenuation length L in XPS analysis were varied within reasonable limits until the calculated retention curves for 1 and 5 keV Ar in Si agreed with experimental data to better than 8%, using the same XPS sensitivity factor throughout. Results: w = 1.4 ± 0.1 nm, Ψ{sub rlc} = 6.6 ± 0.5, and L = 2.7 ± 0.2 nm. Combining experimental and calculated data, it was found that the Ar trapping efficiency of the damaged oxide is intimately correlated with the loss of oxygen. The calculated stationary areal densities of all retained Ar are compared with results obtained by high-resolution medium-energy ion scattering spectrometry. Attractive areas of future research in rare gas retention and nanobubble formation are sketched briefly.

42 ENGINEERING↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2017

Foamed Cement Interactions with CO 2

Geologic carbon storage (GCS) is a potentially viable strategy to reduce greenhouse emissions. Understanding the risks to engineered and geologic structures associated with GCS is an important first step towards developing practices for safe and effective storage. The widespread utilization of foamed cement in wells may mean that carbon dioxide (CO 2 )/brine/foamed cement reactions may occur within these GCS sites. Characterizing the difference in alteration rates as well as the physical and mechanical impact of CO 2 /brine/foamed cement is an important preliminary step to ensuring offshore and onshore GCS is a prudent anthropogenic CO 2 mitigation choice. In a typical oil and gas well, cement is placed in the annulus between the steel casing and formation rock for both zonal isolation and casing support. The cement must have sufficient strength to secure the casing in the hole and withstand the stress of drilling, perforating, and fracturing (e.g. API, 1997, 2010 Worldwide Cementing Practices). As such, measuring the mechanical and properties of cement is an important step in predicting cement behavior under applied downhole stresses (Nelson, 2006). Zonal isolation is the prevention of fluids migrating to different zones outside of the casing and is strongly impacted by the permeability of the wellbore cement (Nelson, 2006). Zonal isolation depends on both the mechanical behavior and permeability (a physical property) of the cement (Mueller and Eid, 2006; Nelson, 2006). Long-term integrity of cement depends on the mechanical properties of the cement sheath, such as Young’s Modulus (Griffith et al., 2004). The cement sheath’s ability to withstand the stresses from changes in pressure and temperature is predominantly determined by the mechanical properties, including Young’s modulus, Poisson’s ratio, and tensile strength. Any geochemical alteration may impact both the mechanical and physical properties of the cement, thus ultimately impacting the structural integrity of the wellbore. In this study, atmospheric foamed cements were generated using a neat cement and three foam qualities (volume of entrained gas in the cement) - 10%, 20%, and 30 % gas volume. The samples were immersed in a 0.25 M NaCl brine followed by the injection of supercritical CO 2 at 28.9 MPa and 50°C. Petrophysical properties were examined for representative samples using computed tomography (CT) and scanning electron microscopy (SEM). CT scanning of representative samples across the range of reacted cements revealed macroscopic changes in structure due to brine/CO 2 /cement interactions. The high foam quality samples resulted in more CO 2 -saturated brine infiltrating radially deeper into the cement and thus were more susceptible to alteration. After 56 days of exposure, the 30% foam quality sample had the most reaction resulting in an alteration depth of 8.35 ± 0.13 mm with a calculated 34.6 ± 0.2% reacted area and 5.76 ± 0.2% reacted pore space area. The neat sample on the other hand, had a reaction depth of 0.31 ± 0.13 mm with a calculated 0.15 ± 0.08% reacted area and 0.57 ± 0.05% reacted pore area. Physical measurements of the exposed samples were consistent with this degree of alteration having 47.02% porosity and the highest permeability of 0.041 mD. These results indicate that the greater surface area provided by the increase of pore space in the higher quality foam coupled with carbonate diffusion reactions enabled greater alteration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Path to Subsurface Closure: When Real and Predicted Data Don't Match, Central Nevada Test Area (CNTA) - 19560

The Central Nevada Test Area, Nevada, was the site of an underground nuclear test in 1968, which resulted in residual contamination at and near the detonation depth of 975 meters (3200 feet). The environmental corrective action strategy for the site was conducted in accordance with the Federal Facility Agreement and Consent Order and all applicable Nevada Division of Environmental Protection (NDEP) policies and regulations. The corrective action strategy initiated in 1999 included developing a conceptual site model, preparing a numerical flow and transport model, calculating contaminant boundaries, negotiating compliance boundaries with NDEP, performing model validation, and monitoring groundwater. This strategy was executed through the Corrective Action Decision Document/Corrective Action Plan phase, when it was determined that the flow and transport model could not be validated because it failed to adequately predict groundwater elevations in the monitoring/validation (MV) wells (MV-1, MV-2, and MV-3) drilled subsequent to the modeling effort. This led to a change in the corrective action strategy that is now designed to validate the compliance boundary through monitoring and institutional controls, rather than relying predominantly on a numerical flow and transport model. This new approach included enhancements to the monitoring well network and additional monitoring to confirm that data were sufficient to proceed to closure. The corrective action activities were completed in 2015, and the Closure Report, which established the long-term postclosure monitoring program, was approved in 2016. The postclosure monitoring program is designed to (1) assess the effectiveness of the compliance boundary by monitoring for the radionuclides of interest and (2) evaluate the effectiveness of monitoring locations within the groundwater flow system by monitoring water elevations to ensure that monitoring wells are located along the potential migration pathways. (authors)

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

Wallula Basalt Pilot Demonstration Project: Post-injection Results and Conclusions

Deep underground geologic formations are emerging as a reasonable option for long-term storage of CO 2 , including large continental flood basalt formations. At the GHGT-11 and GHGT-12 conferences, progress was reported on the initial phases for Wallula Basalt Pilot demonstration test (located in Eastern Washington state), where nearly 1,000 metric tons of CO 2 were injected over a 3-week period during July/August 2013. The target CO 2 injection intervals were two permeable basalt interflow reservoir zones with a combined thickness of ~20 m that occur within a layered basalt sequence between a depth of 830-890 m below ground surface. During the two-year post-injection period, downhole fluid samples were periodically collected during this post-injection monitoring phase, coupled with limited wireline borehole logging surveys that provided indirect evidence of on-going chemical geochemical reactions/alterations and CO 2 disposition. A final detailed post-closure field characterization program that included downhole fluid sampling, and performance of hydrologic tests and wireline geophysical surveys. Included as part of the final wireline characterization activities was the retrieval of side-wall cores from within the targeted injection zones. These cores were examined for evidence of in-situ mineral carbonization. Visual observations of the core material identified small globular nodules, translucent to yellow in color, residing within vugs and small cavities of the recovered basalt side-wall cores, which were not evident in pre-injection side-wall cores obtained from the native basalt formation. Characterization by x-ray diffraction identified these nodular precipitates as ankerite, a commonly occurring iron and calcium rich carbonate. Isotopic characterization (δ 13 C, δ 18 O) conducted on the ankerite nodules indicate a distinct isotopic signature that is closely aligned with that of the injected CO 2 . Both the secondary mineral nodules and injected CO 2 are measurably different from the isotopic content of basalt, injection zone groundwater and for naturally occurring calcite. Final post-injection wireline geophysical logging results also indicate the presence of free-phase CO 2 at the top of the two injection interflow zones, with no vertical migration of CO 2 above the injection horizons. Furthermore, these findings are significant and demonstrate the feasibility of sequestering CO 2 in a basalt formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Investigating the Effect of Sorbed Humic Acid on the Mobility of Uranium - 19461

The Savannah River Site (SRS) was one of the major nuclear facilities for the U.S Department of Energy. During the Cold War, these facilities produced large amounts of radioactive and hazardous waste. On site, the F-Area consisted of three unlined seepage basins that collected approximately 6.8 billion liters of low-level waste that contained nitric acid, radionuclides, and dissolved metals due to plutonium separation operations. Sufficient quantities of uranium isotopes, I-129, Tc-99, and tritium migrated into the groundwater, creating an acidic plume with a pH between 3 and 5.5. Although the basins functioned as designed, the acidic nature of the waste created a source of contamination for the groundwater and soil below the basins. With radionuclides above their maximum contaminant level (MCL) remediation efforts took place. Humic acid as a remediation technology has shown potential for controlling radionuclide mobility. Humic acid is a fraction of humic substance that are ubiquitous in the environment. Humic substances (HS) are complex heterogeneous mixtures of polydispersed materials formed by biochemical and chemical reactions during the decay and transformation of plant and microbial remains. Humic substances account for 50-80% of the organic carbon in the soil or sediment and are known for their excellent binding capacity for metals. Humic acid can be extracted with a potassium hydroxide base solution and then drying it. The high pH solubilizes the humic acid molecules and generates a dark-brown highly-concentrated solution, rich in humic acid, which can be diluted for use. Humic acid's strong sorption capacity and its ability to develop a bond with uranium at slightly acid pH constitute it as remediation technology to reduce uranium (VI) mobility. Column experiments were performed to identify the sorption and desorption of HS in flow-through porous media representative of the acidic aquifer at the F-Area Seepage Basins. Two commercially available HA were selected for comparison, Huma-K and KW15 modified humics (modified humic acid, mod-HA). By conducting a nonreactive rhenium tracer test, on columns filled with soil obtained from SRS's FAW-1 at a depth interval of 18 m -21 m, column parameters such as pore volume (PV), variance, dispersion and peclet number were determined. Columns were then used to estimate the influence of HA on uranium mobility. Data obtained from this study shows that using HA as an amendment for uranium remediation improves uranium removal significantly. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2007

Liquid Water Oceans in Ice Giants

Aptly named, ice giants such as Uranus and Neptune contain significant amounts of water. While this water cannot be present near the cloud tops, it must be abundant in the deep interior. We investigate the likelihood of a liquid water ocean existing in the hydrogen-rich region between the cloud tops and deep interior. Starting from an assumed temperature at a given upper tropospheric pressure (the photosphere), we follow a moist adiabat downward. The mixing ratio of water to hydrogen in the gas phase is small in the photosphere and increases with depth. The mixing ratio in the condensed phase is near unity in the photosphere and decreases with depth; this gives two possible outcomes. If at some pressure level the mixing ratio of water in the gas phase is equal to that in the deep interior, then that level is the cloud base. The gas below the cloud base has constant mixing ratio. Alternately, if the mixing ratio of water in the condensed phase reaches that in the deep interior, then the surface of a liquid ocean will occur. Below this ocean surface, the mixing ratio of water will be constant. A cloud base occurs when the photospheric temperature is high. For a family of ice giants with different photospheric temperatures, the cooler ice giants will have warmer cloud bases. For an ice giant with a cool enough photospheric temperature, the cloud base will exist at the critical temperature. For still cooler ice giants, ocean surfaces will result. A high mixing ratio of water in the deep interior favors a liquid ocean. We find that Neptune is both too warm (photospheric temperature too high) and too dry (mixing ratio of water in the deep interior too low) for liquid oceans to exist at present. To have a liquid ocean, Neptune s deep interior water to gas ratio would have to be higher than current models allow, and the density at 19 kbar would have to be approx. equal to 0.8 g/cu cm. Such a high density is inconsistent with gravitational data obtained during the Voyager flyby. In our model, Neptune s water cloud base occurs around 660 K and 11 kbar, and the density there is consistent with Voyager gravitational data. As Neptune cools, the probability of a liquid ocean increases. Extrasolar "hot Neptunes," which presumably migrate inward toward their parent stars, cannot harbor liquid water oceans unless they have lost almost all of the hydrogen and helium from their deep interiors.

Wiktorowicz, Sloane J.↗
DOE OSTI2018

Final Report: Optimal Model Complexity in Geological Carbon Sequestration: A Response Surface Uncertainty Analysis

The critical component of a risk assessment study in evaluating GCS is an analysis of uncertainty in CO2 modeling. In such analyses, direct numerical simulation of CO2 flow and leakage requires many time-consuming model runs. Alternatively, analytical methods have been developed which allow fast and efficient estimation of CO2 storage and leakage, although restrictive assumptions on formation rock and fluid properties are employed. In this study, an intermediate approach is proposed based on the Design of Experiment and Response Surface methodology, which consists of using a limited number of numerical simulations to estimate a prediction outcome as a combination of the most influential uncertain site properties. The methodology can be implemented within a Monte Carlo framework to efficiently assess parameter and prediction uncertainty while honoring the accuracy of numerical simulations. The choice of the uncertain properties is flexible and can include geologic parameters that influence reservoir heterogeneity, engineering parameters that influence gas trapping and migration, and reactive parameters that influence the extent of fluid/rock reactions. The method was tested and verified on modeling long-term CO2 flow, non-isothermal heat transport, and CO2 dissolution storage by coupling two-phase flow with explicit miscibility calculation using an accurate equation of state that gives rise to convective mixing of formation brine variably saturated with CO2. All simulations were performed using three-dimensional high-resolution models including a target deep saline aquifer, overlying caprock, and a shallow aquifer. To evaluate the uncertainty in representing reservoir permeability, sediment hierarchy of a heterogeneous digital stratigraphy was mapped to create multiple irregularly shape stratigraphic models of decreasing geologic resolutions: heterogeneous (reference), lithofacies, depositional environment, and a (homogeneous) geologic formation. To ensure model equivalency, all the stratigraphic models were successfully upscaled from the reference heterogeneous model for bulk flow and transport predictions (Zhang & Zhang, 2015). GCS simulation was then simulated with all models, yielding insights into the level of parameterization complexity that is needed for the accurate simulation of reservoir pore pressure, CO2 storage, leakage, footprint, and dissolution over both short (i.e., injection) and longer (monitoring) time scales. Important uncertainty parameters that impact these key performance metrics were identified for the stratigraphic models as well as for the heterogeneous model, leading to the development of reduced/simplified models at lower characterization cost that can be used for the reservoir uncertainty analysis. All the CO2 modeling was conducted using PFLOTRAN – a massively parallel, multiphase, multi-component, and reactive transport simulator developed by a multi-laboratory DOE/SciDAC (Scientific Discovery through Advanced Computing) project (Zhang et al., 2017, in review). Within the uncertainty analysis framework, increasing reservoir depth were investigated to explore its effect on the uncertainty outcomes and the potential for developing gravity-stable injection with increased storage security (Dai et al., 20126; Dai et al., 2017, in review). Finally, to accurately model CO2 fluid-rock reactions and resulting long-term storage as secondary carbonate minerals, a modified kinetic rate law for general mineral dissolution and precipitation was proposed and verified that is invariant to a scale transformation of the mineral formula weight. This new formulation will lead to more accurate assessment of mineral storage over geologic time scales (Lichtner, 2016).

58 GEOSCIENCES↗
DOE OSTI2024

Results from an Aeromagnetic Survey to Detect Steel-Cased Wells at a Marcellus Shale Well Site in Washington County, Pennsylvania

Pennsylvania has a 150-year history of oil and gas production—the longest of any state—and this enduring activity has resulted in the drilling of more than 300,000 recorded wells. However, unknown wells likely exist because innumerable wells were drilled during Pennsylvania’s intense early oil and gas history when incomplete records were kept of well locations. There is concern that early wells are likely to be ineffectively sealed because there were no laws that required plugging when the wells were abandoned. Today, many undocumented and unplugged wells are thought to be in areas of emerging shale gas and shale oil development where open wellbores can provide a pathway for undesired upward migration of fluids and gas from hydraulically fractured reservoirs. Due to this concern, Pennsylvania regulators have asked operators to locate orphaned and abandoned wells within a 1,000-ft buffer of proposed new wells. The objective of this report is to demonstrate that high-resolution aeromagnetic surveys, historic air photos, and Light Detection and Ranging (LiDAR) imagery can be rapid and effective methods to reconnoiter large, forested areas of moderate terrain for the presence of abandoned wells. These well-finding methods were evaluated at a proposed Marcellus Shale gas drilling site in Washington County, Pennsylvania, where the methods collectively located 18 confirmed wells: 15 wells were identified from aeromagnetic surveys, two wells were identified from inspection of historical air photos, and one well was identified by evaluation of state-wide LiDAR imagery. Only six wells were previously known, and their locations, as recorded in Pennsylvania’s statewide oil and gas wells database (PA/IRIS/WIS), were often too inaccurate for the wells to be found in the dense underbrush. Twelve wells identified in this study were abandoned, unmarked, and undocumented. Aeromagnetic surveys locate wells by detecting the unique magnetic signature of vertical, steel well casing, which is depicted on magnetic maps as a “bull’s eye” type anomaly that is centered directly over the well. However, when wells were drilled and found to be sub-economic, their casing was sometimes pulled and salvaged for reuse. Such wellbores provide no magnetic response and go undetected if all casing was removed. Oftentimes attempts to retrieve well casing were not 100% successful. For example, historical records for one well in the study area indicate that the well was completed in 1902 as a dry hole and that, to the extent possible, the casing was pulled for reuse. However, a section of 10-in. diameter steel casing was not recovered and remains at an unknown depth in the wellbore. This well was easily detected by the aeromagnetic survey although only deep casing remained in the well. To mitigate for the likelihood that wellbores exist where most or all casing has been removed, this study augmented aeromagnetic data with historic air photos and digital terrain models generated from LiDAR datasets—both databases are publicly available at no cost for areas within Pennsylvania. These complementary methods located three wells where the aeromagnetic anomaly, although present, was subtle and overlooked. Together, these methods determined accurate locations for six known wells within the study area and located 12 previously unknown wells. Although it is not certain that these methods successfully located all wells in the study area, the application of these methods does represent a significant improvement over relying on existing databases for well locations. For the Appendix to the report, see: https://www.netl.doe.gov/energy-analysis/details?id=b46c417a-7c9e-4d25-b810-e6248b0217f4</p>

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2019

Path Forward for Groundwater Compliance Strategy at the Tuba City, Arizona, Disposal Site - 19581

The Tuba City, Arizona, Disposal Site is located in the Navajo Nation and near the Hopi Reservation. A uranium mill operated at the site from 1956 to 1966. During the operational period, mill tailings were slurried to an unlined impoundment. Excess slurry water was transferred to unlined evaporation ponds. Groundwater contamination resulting from the infiltration of fluids in the mill tailings impoundment and evaporation ponds was characterized in the 1990's. The site's Groundwater Compliance Action Plan was prepared in 1999, and a groundwater extraction and treatment system was constructed in 2002. Issues with aquifer productivity and treatment plant reliability led to the suspension of operations in 2014. The limitations of the remediation system prompted a realization that the site conceptual model should be reevaluated. An integrated study to better understand the geochemical and hydrologic factors affecting contaminant transport, to support the development of a new groundwater compliance strategy, was initiated in 2017 and is ongoing. In addition to the groundwater science components, the study included the evaluation of enhanced evaporation methods. The integrated evaluation will lead to refinement of the site groundwater flow model and development of a contaminant transport model and subsequent development of remedial alternatives. Remedial alternatives will be evaluated to determine a new groundwater compliance strategy. The objectives of geochemical testing were to confirm and quantify geochemical reactions controlling uranium fate and transport. These reactions are important because they control the distribution of uranium on the solid phase; potential for release from the solid phase; mobility of uranium in the water phase; and potential for uranium redistribution onto the currently uncontaminated solid phase. A key outcome is the change in understanding of the mass of uranium in the subsurface that is potentially soluble. These results contradict the previous basis for remediation: that uranium was almost exclusively present in soluble form, within the boundaries of the contaminant plume. Testing showed that uranium sorption downgradient from the plume was measurable but relatively minimal. The movement of other mill-related constituents (nitrate, sulfate, alkalinity) results in the dissolution of naturally occurring uranium and a higher than background concentration of uranium in groundwater. These constituents are not subject to sorption, and their downgradient concentrations are reduced only by diffusion and dispersion. The geochemical testing provides conceptual and quantitative information on potential plume migration. While the testing revealed the potential availability of uranium contaminant sources, the aquifer buffering capacity, pH, and presence of oxygen in groundwater may, individually or in combination, be controlling factors for plume movement. Evaluation of evaporative treatment alternatives was undertaken to determine one or more sustainable and efficient methods capable of treating 151 liters per minute (L/min) (40 gallons per minute), the annual average flow rate. In addition to establishing the baseline capacity for evaporation using the site's existing pond, prospective methods of enhanced evaporation included preheating water before inflow to the evaporation pond, vertical circulation in the pond, dye addition, and Wind-Aided Intensified Evaporation. Consideration of a 151 L/min evaporative treatment arose from previous analysis of extraction well production rates and contaminant mass removal, which indicated that about 80% contaminant removal, relative to full-scale operation of the site's groundwater treatment plant, could be achieved with a smaller capacity system. The lower-flow evaporative option also represented a potentially more sustainable treatment system in comparison to restarting or reconfiguring the groundwater treatment plant. Aquifer pumping tests were designed to develop or improve estimates of hydraulic properties such as transmissivity, horizontal and vertical hydraulic conductivity, storativity, and storage coefficient. Testing involved pumping from one well at a constant flow rate for up to 5 days while monitoring groundwater drawdown in the pumping well and five surrounding wells and then shutting down the pumping well to observe groundwater recovery. The aquifer pumping test results provide insight regarding the hydraulics of the saturated zone. The saturated zone is in a low-conductivity formation with preferential flow in the horizontal direction. The aquifer exhibits low vertical conductance. The use of AQTESOLV software indicates that the saturated zone is a layered system due to the cross-bedded sand-grain orientation within the Navajo Sandstone. These indicators of limited vertical flow provide a potential justification for the observed depth of the contaminant plume of 15-46 meters (50-150 feet) below ground surface). Given the relatively low vertical conductance, groundwater extraction should be focused on the upper portions of the saturated zone. Aquifer response testing provided an improved understanding of aquifer characteristics. However, the primary objective, to identify optimal groundwater extraction strategies, has not been fully addressed. A list of path forward activities for hydrologic analyses was developed, with an overarching objective of providing input parameters for groundwater flow and contaminant transport modeling. Defensible predictions of contaminant transport, through modeling, are critical to the development of a new site groundwater compliance strategy. (authors)

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

Generic Post-Closure Safety Assessment for Deep Borehole Disposal of Heat-Generating Waste - 20451

The deep borehole disposal (DBD) concept for high-activity waste consists of drilling a large-diameter borehole into crystalline basement rock to a depth of about 5,000 m, placing waste packages in the lower portion of the borehole - the waste emplacement zone - and sealing and plugging the upper portion of the borehole with a combination of bentonite, cement plugs, and sand or crushed rock ballast. The waste emplacement zone is several times deeper than typical mined repositories and is well below the typical maximum depth of fresh groundwater resources. Recent work describes a generic reference case for deep borehole disposal of Cs and Sr capsules stored at the Hanford Site. Two scenarios are described: an undisturbed (nominal) scenario, and disturbed scenario, in which a waste package is stuck adjacent to a transmissive fracture zone in crystalline basement. Simulations of the processes affecting radionuclide release and transport in the two scenarios provide the basis for quantitative post-closure performance assessment (PA) of the reference case. This paper presents a probabilistic PA of the stuck waste package scenario, i.e. uncertainties affecting radionuclide mobilization and transport are propagated through 200 realizations of a predictive model to quantify uncertainties in predicted radionuclide concentrations within the borehole and the fracture zone. Simulations predict no radionuclide transport in the borehole to distances greater than 100 m above the stuck waste package. An assumed lateral pressure gradient results in predictions of limited radionuclide transport in the fracture zone at a lateral distance of 200 m from the stuck waste package. These results complement previous results for the nominal scenario and contribute to the development of a generic safety case for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2019

Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements

The Aerosol Robotic Network (AERONET) has provided highly accurate, ground-truth measurements of the aerosol optical depth (AOD) using Cimel Electronique Sun–sky radiometers for more than 25 years. In Version 2 (V2) of the AERONET database, the near-real-time AOD was semiautomatically quality controlled utilizing mainly cloud-screening methodology, while additional AOD data contaminated by clouds or affected by instrument anomalies were removed manually before attaining quality-assured status (Level 2.0). The large growth in the number of AERONET sites over the past 25 years resulted in significant burden to the manual quality control of millions of measurements in a consistent manner. The AERONET Version 3 (V3) algorithm provides fully automatic cloud screening and instrument anomaly quality controls. All of these new algorithm updates apply to near-real-time data as well as post-field-deployment processed data, and AERONET reprocessed the database in 2018. A full algorithm redevelopment provided the opportunity to improve data inputs and corrections such as unique filter-specific temperature characterizations for all visible and near-infrared wavelengths, updated gaseous and water vapor absorption coefficients, and ancillary data sets. The Level 2.0 AOD quality-assured data set is now available within a month after post-field calibration, reducing the lag time from up to several months. Near-real-time estimated uncertainty is determined using data qualified as V3 Level 2.0 AOD and considering the difference between the AOD computed with the pre-field calibration and AOD computed with pre-field and post-field calibration. This assessment provides a near-real-time uncertainty estimate for which average differences of AOD suggest a +0.02 bias and one sigma uncertainty of 0.02, spectrally, but the bias and uncertainty can be significantly larger for specific instrument deployments. Long-term monthly averages analyzed for the entire V3 and V2 databases produced average differences (V3–V2) of +0.002 with a ±0.02 SD (standard deviation), yet monthly averages calculated using time-matched observations in both databases were analyzed to compute an average difference of −0.002 with a ±0.004 SD. The high statistical agreement in multiyear monthly averaged AOD validates the advanced automatic data quality control algorithms and suggests that migrating research to the V3 database will corroborate most V2 research conclusions and likely lead to more accurate results in some cases.

David M. Giles↗
DOE OSTI2022

Thermocouple Testing in Support of the AGR-5/6/7 Experiment

This report documents thermocouple testing performed in IRC Lab C-15 over a period of seven years. This testing supported selection and characterization of the thermocouple set used in the AGR-5/6/7 experiment. The following summary was taken directly from the report. Temperature measurement is a challenging aspect of very high temperature irradiation experiments because commonly used high-temperature commercial thermocouples such as platinum-rhodium (Types S, R, and B) and tungsten-rhenium (Type C), suffer dramatic drift because of neutron-induced transmutation. As a result, these types of thermocouples, which are used routinely for industrial temperature measurements outside of reactors, are used only in very special circumstances for reactor experiments. Conversely, because of their low neutron cross-sections, Type N thermocouples are affected to only a limited extent by neutron irradiation. However, the use of these nickel-based thermocouples is limited when the temperature exceeds 1050°C due to drift arising from minor alloying elements migrating from the thermocouple's metal sheath to the thermoelements. This change in the composition of the thermo-elements results in significant decalibration of the signal. The issues described above were recognized during the early planning stages of the final AGR experiment (designated AGR-5/6/7), and a thermocouple furnace testing program was performed over a seven-year period (2014-2019, 2021) to first select and then characterize the best thermocouple set for the high temperature regions of the AGR-5/6/7 irradiation experiment. The calculated temperature range of the AGR-5/6/7 experiment was 600–1500°C. For temperatures below 1000°C standard Type N thermocouples were deemed adequate. The furnace testing campaign identified two thermocouple types suitable for measuring temperatures above 1000°C, a Mo/Nb thermocouple developed at INL called HTIR-TC, and a Type N thermocouple developed by Cambridge University (called herein Cambridge Type N), which featured a custom high nickel alloy sheath. One of the original goals of the furnace testing program was to identify a thermocouple capable of low drift operation near the peak temperature expected in AGR-5/6/7, i.e., about 1400°C. The HTIR-TC design appeared promising in this regard, however a manufacturing difficulty proved to be a barrier and instead the furnace testing focused on drift performance at 1250°C. The manufacturing difficulty was that the Nb sheaths of the HTIR-TCs experienced extreme embrittlement when heat treated at 1600°C or greater. Heat treatment is needed to stabilize the emf output of this TC type, and the higher the heat treatment temperature the higher the peak temperature of stable operation. Because of the sheath embrittlement the heat treatment temperature had to be lowered to 1450°C resulting in a stable operating temperature of about 1250°C. One of the successes of the furnace testing program was identification of a shortcoming in the heat treatment procedure that had been traditionally used in the production of HTIR-TCs. The shortcoming was that the entire heated length of the HTIR-TC sensor was not being heat treated, but rather only the part of the sensor expected to experience temperatures above 1000°C. The problem manifested itself when the thermocouples were removed from the heat treat furnace and placed in another furnace with a different geometry, their indicated temperatures would be widely scattered, but mostly in the negative direction. The solution was to heat treat the entire heated length of the sensor. Since the deepest immersion depth in the AGR-5/6/7 experiment was about 40 inches, a heat treatment length of 48 inches was used. After this change was implemented, thermocouples which were moved into a new environment with a different temperature profile (i.e., a different furnace), produced accurate temperature measurements. Although assembly of the AGR-5/6/7 experiment was completed in September of 2017 (and irradiation begun in 2018), furnace testing of thermocouples continued in 2018 and 2019. The main purpose of this testing was to establish very long-term drift characteristics of the HTIR and Cambridge Type N thermocouples installed in the experiment. Representative thermocouples from the same lots as those installed in the AGR-5/6/7 experiment were used. Additionally, thermocouples of different designs, (particularly variations on the HTIR-TC design) were "piggy-backed" on this testing program to provide insights for instrumenting future very high temperature irradiation experiments. This two-year testing program demonstrated that HTIR-TCs and Cambridge Type N TCs could operate at 1250°C for up to 10,000 hrs (and in some cases longer) while experiencing negative drifts on the order of 2-4°C/1000 hrs. This performance was considered acceptable given the extreme operating environment the sensors faced.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗