Search NASA⌕ Search

SEARCH · Search NASA

Results for “TDS”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

How oxygen passivates polycrystalline nickel surfaces

The passivation of polycrystalline nickel surfaces against hydrogen uptake by oxygen is investigated experimentally with low energy ion scattering (LEIS), direct recoil spectroscopy (DRS), and thermal desorption spectroscopy (TDS). These techniques are highly sensitive to surface hydrogen, allowing the change in hydrogen adsorption in response to varying amounts of oxygen exposure to be measured. The chemical composition of a nickel surface during a mixed oxygen and hydrogen exposure was characterized with LEIS and DRS, while the uptake and activation energies of hydrogen on a nickel surface with preadsorbed oxygen were quantified with TDS. By and large, these measurements of how the oxygen and hydrogen surface coverage varied in response to oxygen exposure were found to be consistent with predictions of a simple site-blocking model. This finding suggests that, despite the complexities that arise due to polycrystallinity, the oxygen-induced passivation of a polycrystalline nickel surface against hydrogen uptake can be approximated by a simple site-blocking model.

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

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Initial study on thermal stability of cold spray tantalum coating irradiated with deuterium for fusion applications

Abstract Removal of neutral hydrogen atoms in the plasma edge reduces the number of charge exchange events and thus, the net energy losses in the plasma, significantly improving performance of fusion devices. Effective control of the residual pressure of hydrogen isotopes (HIs) in the plasma edge may be achieved by utilizing a hydrogen absorbing first wall interface capable of withstanding the harsh fusion environment. In this study, we have investigated tantalum (Ta) coating deposited by cold spray technology on 316L stainless steel substrate as a potential plasma-facing material surface. High fluence low energy deuterium plasma irradiation experiments and subsequent thermal annealing cycles associated with thermal desorption spectrometry (TDS) demonstrated superior structural stability of the Ta coating. TDS experiments revealed the outgassing of deuterium (as measure of its retention) for cold spray Ta coatings to be three times higher than bulk Ta and two orders of magnitude greater than bulk polycrystalline W. X-ray photoelectron spectroscopy revealed evolution of oxidation states upon deuterium irradiation and a partial recovery of the metallic signature of Ta after the thermal treatment at 1100 K.

Physics↗
DOE OSTI2022

Reduced defect recovery in self-ion damaged W due to simultaneous deuterium exposure during annealing

Deuterium (D) plasma exposure during annealing of self-ion damaged tungsten (W) is shown to exhibit reduced defect recovery when compared to annealing without D plasma exposure. In these experiments, samples were first damaged with 20 MeV W ions. Next, samples were annealed either with or without simultaneous D 2 plasma exposure. The simultaneous annealed samples were first decorated by D 2 plasma at 383 K prior to ramping up to an annealing temperature of 473, 573, 673, or 773 K and held for 1 h with concurrent plasma exposure. The vacuum annealed samples each had a corresponding temperature history but without D 2 plasma treatment. Finally, all samples were exposed to D 2 plasma at 383 K to decorate any remaining defects. Nuclear reaction analysis and thermal desorption spectroscopy (TDS) shows that the simultaneous plasma-exposed and annealed samples exhibited virtually no defect recovery at annealing temperatures of up to 673 K, and had higher D retention than found in the vacuum annealed samples. TDS results indicate that only the lowest detrapping energy defects recover at an 773 K anneal for the simultaneous plasma annealed samples, while the vacuum annealed samples showed defect recovery at all anneal temperatures. Furthermore, this experiment clearly demonstrates that D occupied defects can significantly reduce or eliminate defect annealing in W, and is consistent with the existence of synergistic plasma exposure/displacement damage effects in fusion-energy relevant plasma facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observation of topological Dirac fermions and surface states in superconducting BaSn3

The interplay between topological electronic structure and superconductivity has attracted tremendous research interests recently as they could induce topological superconductivity (TSCs) which may be used to realize topological qubits for quantum computation. Among various TSC candidates, superconducting BaSn 3 (T c ~ 4.4K) has been predicted to be a topological Dirac semimetal (TDS) hosting two pairs of Dirac points along the Γ-A direction. Here, by combining the use of angle-resolved photoemission spectroscopy and ab initio calculations, we identified the predicted topological Dirac fermions and confirmed the TDS nature of the compound. In addition, we observed surface states connecting the Dirac points. Our observations demonstrate BaSn 3 as a superconductor with nontrivial topological electronic structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Terahertz time-domain spectroscopy imaging of pancreatic ductal adenocarcinoma tissues

Pancreatic ductal adenocarcinoma (PDAC) ranks among the malignancies with the highest fatality and morbidity rates. This is predominantly attributable to an absence of understanding the intricate and diverse microenvironment of the tumor. We use terahertz time-domain spectroscopy (THz-TDS) imaging in transmission geometry to probe ex-vivo the heterogenous microenvironment of the genetically modified murine PDAC tissue that closely resembles the PDAC heterogeneity in human malignancy. We introduced a maximum a-posteriori probability estimation algorithm to objectively the tumor’s heterogenous microenvironment using the average values of refractive index and absorption coefficient within the useable terahertz bandwidth as imaging markers. Furthermore, direct comparison of stained histopathologic images and the refractive index and the absorption coefficient high-resolution, two dimensional maps of the same PDAC samples confirms the high potential of the THz-TDS method for tumor tissue characterization.

absorption↗
DOE OSTI2025

High-speed fiber-coupled terahertz time-domain spectroscopy system using a bias-free plasmonic source

We demonstrate a high-speed fiber-coupled terahertz time-domain spectroscopy (THz-TDS) system using a bias-free plasmonic terahertz source. The system is powered by a compact, low-power femtosecond fiber laser oscillator that drives a dual-branch fiber amplifier to pump and probe a plasmonic terahertz source and detector, respectively. A polarization-maintaining, fiber-coupled optical delay line is employed to precisely control the temporal offset between the pump and probe beams through a compact, robust platform. Operating in the self-similar amplification regime, the fiber amplifiers are engineered to mitigate nonlinear pulse distortion and dispersion while delivering robust, high-energy pulses at a 1560-nm wavelength, with an output power exceeding 326 mW. In reflection mode, the THz-TDS system achieves a peak dynamic range of 81 dB and a bandwidth greater than 3.5 THz at a scan rate of 3.75 Hz. The utility of the system for rapid sensing is further demonstrated by quantifying the porosity and density of pharmaceutical tablets through measurements of their physical thickness and effective refractive index.

Jiang, Xinghe [University of California, Los Angel↗
DOE OSTI2020

Revised Kimberlina Groundwater Simulations

Kimberlina Groundwater Qaulity Simulations (Version 3: 2015 October) by Lawrence Livermore National Labratory. Submitted to the EDX team at the National Energy Technology Laboratory for research use. This dataset contains 562 updated groundwater quality simulations based on the Kimberlina site. It is a six-layer model including the aquifers and aquitards inbetween the ground surface and the reservoir. The input files and associated leakage information for these simulations are included in the zip file "nrap_kim_well31_fluxes and inputs.tar.gz". These 562 simulations are the successful results out of 1000 monte carlo simulations, with time length greater than 5 years. Each simulation consists 6 output files, therefore, there are 3372 files (6*562) in total. One single example simulation of full time length (200 year) is also provided in the zip file "single_simulation_200yr_sim0012.tar.gz". These 6 groundwater quality parameters are: - pH (use absolute value) - TDS (mg/kg) - Pressure (Pa) - CO2 Liquid Saturation (fraction range 0.0-1.0) - CO2 Gas Concentration (mol/L) - CO2 Liquid Concentration (mol/kg) Their correponding file names, take the first simulation for example, are: ph0001.ntab (for pH) , tds0001.ntab (for TDS), usnt0001.p.ntab (for pressure), usnt0001.s_liquid.ntab (for liquid saturation), usnt0001.c_co2_gas.ntab (for CO2 gas concentration) and usnt0001.c_co2_liquid.ntab (for CO2 liquid concentration). Each output file contains at least 19 columns (variables). The total number of columns depends on the simulation time length. The column data include the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the jth grid of y-axis 4. k: the kth grid of z-axis 5. element_ref: element type reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The rest (14, 15, 16...) variables are the simulated parameter values over time, take pH as an example, are: 14. ph000y: initial pH value. 15. ph001y: simulated pH value at the end of the 1st year. 16. ph002y: simulated pH value at the end of the 2nd year. ... (and so on)... The zip input files have the input fluxes of brine and CO2 (in kg/s) allocated to model layers 1,2,3 and 5 (layers 4 and 6 are assumed to be no-flow aquitards) through the following nodes: (src1 WELL 7 7 25 25 10 32) #Etchegoin (src2 WELL 7 7 25 25 33 72) #Macoma Chanac (src3 WELL 7 7 25 25 73 83) #Santa Margarita Mclure (src4 WELL 7 7 25 25 84 89) #Fruitvale Round Mountain (aquitard) (src5 WELL 7 7 25 25 90 96) #Olcese Sand (thief zone) (src6 WELL 7 7 25 25 97 99) #tremblor Freeman (aquitard) ex: brn1 is assigned to src 1 from i=(7,7), j=(25,25), k=(10,32) Edited by Ya-Mei Yang at National Energy Technology Laboratory

brine leakage↗
DOE OSTI2022

Maximizing the Proppant Carrying and Viscoelastic Properties of the Bakken Hypersaline-Produced Water with High-Viscosity Friction Reducers for Sustainable Applications

Summary The development and production of unconventional reservoirs, such as the Bakken Formation, have become a resolved mystery for operators in North America since the arrival and advancement of horizontal drilling and hydraulic fracturing technologies. As a result, unconventional reservoir assets became the central focus of the oil and gas industry at the state, national, and global levels. The produced water from these activities in the Bakken Formation have high salt contents (110,000–350,000 ppm) total dissolve solids (TDS) and can pose significant challenges to the environment if not treated. Deep injection into disposal wells is the routine method used to get rid of the Bakken produced water. However, there have been some concerns that unrestrained injections, in addition to polluting the groundwater, could potentially lead to seismic activities either at the time of injection or in the near future. To diminish the environmental impacts that may be associated with induced seismicity, including the reduction of the costs of water acquisition, the produced water can be treated and reused in the hydraulic fracturing processes. Also, the treated water could be used for irrigation purposes, for power generation, and coal mining operations. The issues of waste water and residual oil high in TDS are challenges yet to be effectively addressed despite preceding research and studies on advancing produced water technologies. The goal of this study is to explore all applicable ways by which the produced water from the Bakken Formation can serve as a replacement base fluid for use with polymers like the high-viscosity friction reducers (HVFRs) to create hydraulic fracturing fluids that can be stable with reservoir conditions and also be able to minimize environmental impacts and cost of operations. Experimental investigations using the high-salinity produced water from the Bakken Formation with HVFRs were carried out. The studies included a base case that served as a bench mark for comparing the effectiveness of the other scenarios. The results indicate that the Bakken hyper-saline produced water can withstand effect of heavy metals, salinity, hardness and remain stable through different shear rates (66–330 s−1) when treated with higher dosages [4–8 gal/1,000 gal (gpt)] of HVFRs. Filtration and dilution were the only methods used on the Bakken Formation produced wate for this research

Engineering↗
DOE OSTI2020

Intensified Flue Gas Desulfurization Water Treatment for Reuse, Solidification, and Discharge

The goal of this final project report is to summarize the work conducted on project DE-FE0031555. In accordance with the Statement of Project Objectives (SOPO), the University of Kentucky Center for Applied Energy Research (UK CAER) (Recipient) developed a hybrid electrochemical-membrane technology for treating wet flue gas desulfurization (WFGD) blowdown from coal-fired power generation plants that (1) removes government-regulated species and total suspended solids (TSS) via electrocoagulation with solid/liquid separation to avoid further treatments by a biological method specifically for selenium, and (2) reduces total dissolved solids (TDS) in a nanofiltration unit to promote water recycle and reduce freshwater intake. Additional technological configurations were pursued, including (1) nanofiltration combined with zeolite dewatering for a near-zero liquid discharge (ZLD) scenario, and (2) polishing discharge from zeolite dewatering using capacitive deionization (CDI) with 90+% salt rejection to produce water with less than 100ppm TDS for more beneficial use. Project results verified that UK CAER’s hybrid approach addressed the complexities of WFGD blowdown with simplified unit operation at relatively low capital investment by demonstrating: (1) an effective in-situ generated solid sorbent material, green-rust (GR), that specifically removes regulated dissolved species in WFGD blowdown, including selenium, arsenic, and nitrates; (2) a low-cost, high flux nanofiltration treatment that removes both monovalent and divalent species at >80% rejection to adequately provide water recycle and extend the life of FGD equipment by mitigating chloride aided pitting corrosion mechanisms; and (3) a solidification and stabilization process using pozzolanic agents that retain regulated species and satisfy the Resource Conservation and Recovery Act (RCRA) leachate requirements. UK CAER proved electrocoagulation with GR to be an intensified solution for effluent limit guidelines (ELG), and based on preliminary techno-economic analysis (TEA), in comparison to DOE’s Case 1 for Biological Wastewater Treatment that had a purchased equipment cost (PEC) of $\$$5.26 million and a Total plant cost (TPC) of $\$$30.0 million, the EC process had a PEC of $\$$2.67 million and the TPC of $\$$17.0 million. While the electricity consumption for Case 1 was 53 kW, it was 110 kW for the EC process due to dissolving iron for a three-stage reaction design. However, electricity consumption for the EC process could be significantly reduced by adopting a counter-current reactor design for at least a 6-stage chemical reaction operation.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

NRAP-Open-IAM: FutureGen2 Component Models

This report describes the development and testing of three component models for NRAP-Open-IAM, the National Risk Assessment Partnership’s open-source integrated assessment model. The FutureGen2 Lookup Table Reservoir component model is based on interpolation of data from a set of lookup tables. The lookup tables contain pressures and saturations predicted by multiphase flow simulations performed with a reservoir simulator. The FutureGen2 Above Zone Monitoring Interval (AZMI) component is a surrogate model that can be used to estimate the impact that carbon dioxide (CO 2 ) and brine leaks from the CO 2 storage reservoir at the FutureGen 2.0 site might have had on overlying aquifers or monitoring units were a leak to occur. The model estimates the size of “impact plumes” according to five metrics: pH, Total Dissolved Solids (TDS), pressure, dissolved CO 2 and temperature. The FutureGen2 Aquifer component is similar, but is limited to four metrics: pH, Total Dissolved Solids (TDS), pressure, and dissolved CO 2 . The input parameters for each model are the same, but the Aquifer component is applicable to depths between 100 m and 700 m and the AZMI component is applicable from depths between 700 m and 1050 m.

42 ENGINEERING↗
DOE OSTI2022

An Innovative Zero-Liquid Discharge Intermediate-Cold-Liquid Eutectic-Freeze Desalination System

A novel freeze desalination technology is presented for purification of water produced from various industrial processes, including oil and gas extraction. The main innovative idea is introducing a water-immiscible intermediate-cold-liquid (ICL) medium which absorbs the cold from a refrigeration cycle and transfers it to a freezing chamber where it is mixed with brine. The cold ICL enters the freezing chamber at temperatures between -10°C to -30°C. Precooled brine in injected into the freezing chamber and pure water freezes out of the brine upon mixing with the ICL. The rate of freezing inside the chamber is controlled by adjusting the flow rate of the cold ICL relative to the inlet brine. Both eutectic freeze and non-eutectic freezing are possible by controlling the temperature of the freezing chamber. The solid-liquid slurry created in the freezing chamber flows to a filter separation unit after passing through a separation column. The ice crystals are separated from the liquid by a rotary filter device. The collected ice is rinsed by a centrifugal unit. The rinsed ice is melted to produce fresh water. The developed technology resolves the issues common to current freeze-desalination systems while offering superior heat transfer performance ascribed to direct-contact between the cooling medium and the brine in addition to producing solid salt and pure water. The system operates under atmospheric pressure. The low operation temperature minimizes the corrosion issues common in systems exposed to brine. Also, due to the low operating temperature, there are no concerns about evaporation of volatile organics from the produced water during the treatment. The developed technology effectively eliminates the use of chemicals for pretreatment of brine. The system can be constructed from commercially available components. Extensive testing using a lab-scale prototype by using synthetic brine and actual produced water with TDS values from 50,000 ppm to 300,000 ppm showed that freshwater could be achieved with relatively small dependency on the quality of the input brine. The quality of the treated water was found to be mainly a function of the duration of the centrifugal process. Treated water with TDS as low as 200 ppm could be achieved in the tests. Thermo-economic modeling coupled with system level thermochemical process model were developed for prediction of the LCOW. The LCOW was predicted to be $\$$0.50 to $\$$0.90 per barrel of produced water. The potential markets for deployment of the technology are the oil and gas, in-land thermoelectric powerplants, RO desalination plants, and mining industry.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Geochemical characterization of lithium deposition in fossil energy wastes

Lithium is a critical mineral used in rechargeable batteries for electric vehicles (EVs) and future modernization of electric grids with sensitive supply chains that are subject to volatility. Therefore, methods to recover lithium from domestic unconventional sources are actively being pursued in the United States. Herein, we report on the lithium recovery potential from four fossil energy waste feedstocks: oil and gas drill cuttings (sample number n = 16); oil and gas produced waters (n > 200) from Marcellus, Bakken and Permian basins; coal byproducts (n>20), such as coal ash; and acid mine drainage treatment solids (AMD solids). Select solid samples underwent sequential extraction to explore lithium hosting phases and simple organic and inorganic acid extraction to explore lithium recovery potential. Our preliminary results found lithium concentrations up to 440 mg/kg in select AMD treatment solids and up to 300 mg/L in Marcellus Shale produced waters, demonstrating these sources have lithium yield potentials comparable to conventionally mined lithium ores and brines. Sequential extractions revealed that lithium is mostly associated with clay and silicate phases in shale drill cuttings, whereas lithium resides in Fe,Mn-oxide reducible phases in AMD solids. Further, the lithium content in produced waters has a strong linear relationship with total dissolved solid (TDS) levels in three separate oil and gas basins. At the same TDS level, Marcellus Shale produced waters contain more lithium compared to Bakken Shale and Permian Basin waters, with higher percentages of Ca and Mg, major cations that might impact lithium recovery efficiencies. Our study demonstrates the heterogeneity of lithium hosts from different fossil energy wastes. Characterization results will inform future lithium recovery from these different feedstocks.

Stuckman, Mengling↗
DOE OSTI2024

Geochemical and Microbial Dynamics of Hydrogen in a Methane Storage Reservoir

Hydrogen has been identified as a flexible energy carrier with zero or negative emission across multiple energy systems, and existing natural gas infrastructure could be leveraged if hydrogen gas (H2) was blended with methane (CH4). For example, subsurface methane storage reservoirs could be slightly modified to also store hydrogen if a methane/hydrogen blend were injected. However, the compatibility of methane storage reservoirs to include H2 injection has not been fully demonstrated, and this could lead to geochemical and microbiological reactions that alter the reservoir and stored gas content. It is essential that we understand the impact of H2 gas on the biogeochemistry of subsurface storage reservoirs before deploying large-scale H2-CH4 storage, We collected produced fluid from two separate methane storage reservoirs in the Southwestern US. First, we completed a baseline analysis of the biogeochemistry through qPCR, 16S rRNA sequencing, metagenomic sequencing, and geochemical analysis. Each reservoir was found to have unique geochemical conditions and a unique microbial community structure, with Site 1 having a higher TDS and an abundance of Shewanella and Site 2 having a lower TDS and high abundance of Eubacterium and Acetobacterium. Next, we ran a series of high pressure, high temperature reactors under hydrogen storage conditions with the biological sample from one of the storage reservoirs and a 20% H2-80% CH4 gas blend for up to 7 days. Our results show a decrease of hydrogen by 5% in reactors as early as 1-3 days. Previous hydrogen storage work has linked subsurface microorganisms with methanogenesis hydrogen sulfide production, acid production, and microbial corrosion. Our results show minimal change in the fluid chemistry, with the exception of a decrease in dissolved sulfate concentrations. Taxonomic sequencing demonstrated the presence of microorganisms capable of iron redox, acid generation, and hydrogen sulfide production throughout the reactors, suggesting microbial hydrogen consumption may occur through various metabolic pathways. This work demonstrates that site-specific geochemistry and microbiology may impact the efficiency of hydrogen storage in methane storage reservoirs.

environmental microbiology↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Proposed Potential Mitigation of Wastewater Disposal Through Treated Produced Water in Bakken Formation

The recent global expansion in the development of unconventional oil and gas assets has also resulted in a tremendous increase in the number of extended horizontal drilling and hydraulic fracturing projects in the Bakken. The United States is presently the largest global crude oil producer, and the Bakken Formation in North Dakota is one of the major contributors to this achievement. However, the wastewater produced from these increased oilfield activities are highly saline (∼170,000 to 350,000 ppm TDS) and no technology currently available can satisfactorily treat it. As a result, more than 90% of wastewater in the Bakken is disposed of by deep injection into disposal wells. However, there are growing environmental and operational concerns about the sustainability and impacts of this approach. Research has shown that cumulative wastewater injection in some areas could increase the chances of earthquakes in those areas. However, if this produced water is efficiently treated, it could be reused in hydraulic fracturing operations or to support coal mining and irrigation activities. All these applications would reduce the need for wastewater injection and reduce the demand for fresh water used in hydraulic fracturing operations across North Dakota. For this purpose, we propose an enhanced supercritical technology we call Supercritical Water Extraction – Enhance Targeted Recovery to handle the issue of high TDS of flowback and produced water in the Bakken.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2021

Southwest Pacific tropical cyclone development classification utilizing machine learning and synoptic composites

This study evaluates the ability of machine learning algorithms to classify tropical depressions (TDs) and tropical storms (TSs) in the western region of the southwest Pacific Ocean (SWPO). Decision rules are generated to predict the environment required for a depression to fully develop into a mature storm, and the most influential predictors in the classification decision are ranked. TD and TS are discriminated based on a maximum sustained wind speed threshold (≥17 ms -1 ). Various aerosol, thermodynamic, and dynamic parameters are extracted closest to the initiation point of each non-developing and developing sample. The covariates associated with each labelled sample are used to train a decision tree and random forest model. Results using a testing dataset suggest the random forest approach more accurately distinguishes between non-developing and developing samples. The classification accuracy of the decision tree and random forest are 72% and 91%, respectively. Random forest outperformed the decision tree by providing higher accuracy in test data. The most important variables for binary classification are sea salt aerosol optical depth (AOD), 1,000 mb relative humidity, and sea surface temperature. AOD is a quantitative estimate of the aerosols presents in the air through the extinction of a ray of light as it passes through the atmosphere. Mean composite maps constructed in an unsupervised manner have been created for the most important variables identified by the random forest classifier during TD and TS events to highlight the difference in geophysical and aerosol variables' climatology during the two different classifications. This work will advance the risk management strategies for northeastern Australia and other SWPO basin islands to control their tropical cyclone related losses through prioritizing forecasting variables that are the strongest predictors of the strengthening of tropical depressions into tropical cyclones.

54 ENVIRONMENTAL SCIENCES↗