Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “brine”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

SALSA_python (SALSython) v1

SALSA_python (Semi-Analytical Leakage Solutions for Aquifers) is a software that computed semi-analytical solutions for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. It can simulate brine leakage into aquifers in a multi-aquifer-aquitard system with multiple injection and leaky wells. This situation is encountered in underground CO2 storage wherein brine leakage from pressurized reservoirs into aquifers is of concern. SALSA_python calls the original SALSA[1] subroutines in python by using the salsa2.so library. This enables the incorporation and coupling of SALSA computations into existing python-based codes and tools. SALSA_python is available for Linux and Mac operating systems. [1] Cihan, A., Oldenburg, C. M., & Birkholzer, J. T. (2022). Leakage from coexisting geologic forcing and injection-induced pressurization: A semi-analytical solution for multilayered aquifers with multiple wells. Water Resources Research, 58, e2022WR032343.

Bhuvankar, Pramod↗
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 OSTI2020

Computational and Experimental Characterization of Intermediate Amorphous Phases in Geological Materials

In the subsurface, MgO engineered barriers are employed at the Waste Isolation Pilot Plant (WIPP), a transuranic waste repository near Carlsbad, NM. During service, the MgO will be exposed to high concentration brine environments and may form stable intermediate phases that can alter the barriers effectiveness. Here, MgO was aged in water and three different brine solutions. X-ray diffraction (XRD) and 1 H nuclear magnetic resonance (NMR) analysis were performed to identify the formation of secondary phases. After aging, ~4% of the MgO was hydrated and fine-grained powders resulted in greater loss of crystallinity than hard granular grains. 1 H magic angle spinning (MAS) NMR spectra resolved minor phases not visible in XRD, indicating that diverse 1 H environments are present along with Mg(OH) 2 . Density functional theory (DFT) simulations for several proposed Mg-O-H, Mg-CI-O-H, and Na-O-H containing phases were performed to index peaks in the experimental 1 H MAS NMR spectra. While proposed intermediate crystal structures exhibited overlapping 1 H NMR peaks, Mg-O-H intermediates were attributed to the growth of the 1.0-0.0ppm peak while the Mg-CI-O-H structures contributed to the 2.5- 5.0ppm peak in the chloride containing brines. Overall, NMR analysis of aged MgO indicates the formation of a range of possible intermediate structures that cannot be resolved with XRD analysis alone.

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

NRAP-Open-IAM: Open Wellbore Component (V.2.0)

The Open Wellbore component of NRAP-Open-IAM is applicable to the calculation of Area of Review at a carbon storage site and fast estimation of dense gas dispersion from multiple continuous CO 2 surface leakage sources for risk assessment. The Open Wellbore Component of NRAP-Open-IAM is a lookup table model based on the drift-flux approach. The National Research Assessment Partnership’s second-generation integrated assessment model, NRAP-Open-IAM, is open-source software written in Python for use in performance and quantitative risk assessment of geologic carbon storage (GCS) systems. The look-up table for open well leakage is an updated version the previous model which provides finer resolution of reservoir depth and more accurate fluid properties than the previous open well look-up table, especially in calculating the gas (CO 2 -rich) phase properties. The component model input parameters include the reservoir transmissivity, aquifer transmissivity, brine salinity, well radius, and well top and bottom depths. The reservoir transmissivity and brine salinity were set to values consistent with the reservoir component, and the well top and bottom depths varied with location. The possible outputs from the Open Wellbore component are leakage rates of CO 2 and brine to either an aquifer or atmosphere, depending on the depth of the well top. Several python test scripts are presented that demonstrate that the Open Wellbore component works correctly. Use of the open wellbore model in NRAP-Open-IAM may result in a large CO 2 leakage rates, comparable to the leakage rates of CO 2 blow out.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Flue-Gas Desulfurization Effluent Management using an Innovative Low-Energy Biosorpotion Treatment System to Remove Key Contaminants

Among the most critical water contaminants of concern affecting wide geographical regions and a number of industries and natural systems is selenium. Selenium found in surface, ground and wastewater in originates from natural sources, as well as industrial sources such as petroleum refineries, electronics manufacturing, pesticides, and coal power plants and mining also contribute to selenium contamination in water in the US. At high concentrations, selenium is toxic to human and wildlife. There are a number of technologies that have been used to treat selenium and other similar contaminants in water. Biological treatment of selenium has been used in the past to reduce soluble SeVI and/or SeIV to insoluble Se0, which is then filtered in the same vessel. The insoluble selenium (Se0) is then backwashed from the system and solids are separated for subsequent disposal, if they meet the leaching and water content criteria. In order to promote biological reduction to insoluble Se0, heating of bioreactor is needed in some applications, and excess food source (electron donor) is added so that all selenium can be filtered. An additional disadvantage of these systems is the significant amount of water lost due to extensive and frequent backwash and rinse cycles. When comparing the advantages and energy requirements of the various treatment technologies, RO membrane filtration immediately stands out due to the excessive energy expenditure needed to pump water across the membrane although RO is an effective way to remove selenium. In addition, RO requires extensive pretreatment, such as MF membrane, and frequent maintenance, rendering it an expensive option that may be out of reach for certain applications. In fact, although the performance was good during the pilot testing by the NSMP Working Group for treatment naturally occurring selenium in the surface water, the high electricity requirements and significant reject water stream made it an infeasible alternative. While conventional ion exchange maybe an effective treatment option, it requires frequent regeneration of the resin when applied to highly contaminated water, which leads to several tons of contaminant-laden, high-salinity brine that needs to be disposed off-site each day. One of the water systems in the west coast currently uses ion exchange for selenium treatment and has been trucking selenium laden hazardous brine waste weekly in the last several years. Pneumatic pumping and rinse water pumping required for ion exchange also increase the energy usage. In comparison, adsorption process is a passive treatment system where contaminated water comes in contact with an adsorption media in a vessel. Typically, there is no mixing, backwash, or recycle pumping required, thus significantly reducing the energy usage. A passive single-use adsorption system does not require backwash, thereby generating small amount of process waste, and producing the highest water yield among the alternatives. The energy and water efficiencies, and applicability for SeVI and SeIV are summarized in Table 1. Despite these benefits though, adsorption typically does not work for the most oxidized form of selenium (SeVI). The innovative biosorption process integrates both process to increase the treatment efficiency while minimizing energy, chemical, and time required to treat both SeVI and SeIV. Additional advantages include simple partial biological reduction with reduced on-site waste generation, which lead to water and electricity savings, and less operational need compared to biological treatment alone. This makes biosorption especially suitable for remote areas, where liquid backwash and brine disposal may be cost prohibitive or infeasible.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Recovery of Rare Earths, Precious Metals and other Critical Materials from Geothermal Waters with Advanced Sorbent Structures - CRADA 355 (Abstract)

The ability to recover valuable trace level minerals from geothermal brines using high-performance solid-phase sorbents will be explored and developed. A compressive range of sorbent materials will be screened for application to metal extraction from geothermal brines. Preferred sorbents from extraction of trace levels of rare earths (REs), precious metals (PMs), and other critical/strategically valuable materials (CMs) such as Zn, Mn, Te, Sc, Se and U from geothermal brines will be identified. For the preferred sorbents PNNL will determine the volumes they are capable of providing efficient extraction from. The thermal and chemical limits (including, acid, sulfur, salt) for performance of the preferred sorbent materials will be determined; with a target of at least 125°C and perhaps as high as 400° C. Sorbent form factors (including, packed bed, membrane, mats) that can function efficiently and be installed cost effectively in geothermal energy plants will be assessed for chemical and economic viability. Material regeneration and cyclic utilization will be demonstrated, targeting hundreds to thousands of cycles. Options for recovery and purification (including, selective separation of heavy REs) of collected materials will be explored. A techno-economic analysis (TEA) will be performed to assess the best approach to provide a value-added extraction process for geothermal energy systems. The sorbent materials and engineering analysis will be applicable to other industrial processes in which secondary recovery of valuable materials could provide economic benefit.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Novel Zwitterionic Chromatography to Separate Lithium from Unconventional Resources

Since lithium (Li) is a key element for clean energy technologies, the global lithium demand is anticipated to increase rapidly. As a result, efficient lithium extraction technologies that allow the exploitation of unconventional lithium resources, such as geothermal brines and inland brine streams, are essentially needed to meet the Li demand and keep supply chains stable. However, because of low Li concentration, low ratio of Li/Mg or Li/Ca, and complex feed compositions in these feedstocks, using conventional hydrometallurgy or current Direct Lithium Extraction (DLE) processes are often chemical-intensive, requiring substantial amounts of reagent chemicals to recover Li at scale, resulting in large chemical footprints. Therefore, we developed a novel Zwitterionic Chromatography (ZIC) process to separate Li from these resources. Because Li can be partitioned from divalent salts on ZIC under water elution without using chemicals, exploiting the unconventional resources via a ZIC approach can be environmentally benign and sustainable. In this work, we will present the development of the ZIC process for DLE applications. Simulation studies for the Li separation mechanism in ZIC and the demonstration of a continuous ZIC process to separate Li from domestic inland brine and mining wastewater will be presented. Finally, the chemical and energy footprints of ZIC in comparison to other DLE processes will be discussed.

09 BIOMASS FUELS↗
DOE OSTI2021

Multi-fluid renewable geo-energy systems and methods

A geo-energy production system and method extracts thermal energy from a reservoir formation, and stores either thermal waste heat or excess heat in a storage zone of the reservoir formation. A compressed fluid injection injects an unheated, compressed working fluid into the storage zone. A fluid injection well injects a working fluid laden with thermal waste heat or excess heat into the storage zone. The storage zone is located below a caprock layer and above a native brine zone of the reservoir formation and is partially circumscribed by a hot brine storage zone. The compressed working fluid assists with a withdrawal of pressurized brine residing below and/or to the sides of the storage zone. A compressed CO2, N2, or air production well helps to remove compressed working fluid from the storage zone for use in power production.

Buscheck, Thomas A.↗
DOE OSTI2021

A simulation study of carbon storage with active reservoir management

As part of the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project led by Battelle Memorial Institute, a study was conducted to determine the feasibility of storing carbon dioxide (CO 2 ) in the stacked saline rock formations of the Sleepy Hollow Field (SHF), located in Red Willow County, southern Nebraska. A series of CO 2 injection simulation scenarios, with and without active reservoir management (ARM; brine extraction), were evaluated to investigate the feasibility of storing 50+ million tonnes (Mt) of CO 2 . The results indicated CO 2 injection combined with ARM may enable permanent storage of 50+ Mt of CO 2 . The area of review (AOR), the area in which underground sources of drinking water (USDWs) might be endangered during CO 2 injection, was assessed for the simulation scenarios. In comparison to a case without ARM, brine extraction resulted in a much smaller AOR, covering an area of 42 square miles (108.8 km 2 ), roughly one-fourth the size of an AOR resulting from CO 2 injection without ARM. Our findings, presented in this paper, indicate that CO 2 injection with ARM can improve the CO 2 storage capacity of a geologic storage complex up to 100% while also reduce the rate of pressure buildup in the subsurface, resulting in a 75% reduction in AOR. This may help in lowering carbon capture and storage project costs, risks, and effort needed to meet monitoring requirements for a storage project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Determination of transition metal ions in fossil fuel associated wastewaters using chelation ion chromatography

Here, this study outlines the development and subsequent validation of a method using chelation ion chromatography (CIC) pretreatment followed by traditional ion chromatography (IC) and post column UV/vis detection to measure transition metals in fossil fuel wastewaters, such as oil & gas (O&G) brines and coal mine drainage (CMD) waters. Measurement of transition metals is often an important characterization step in the research of environmental and energy systems. IC represents one way to measure these metals with the advantages of being versatile, simple and relatively low cost compared to other analytical methods. However, high concentrations of alkali and alkaline earth metals present in fossil fuel wastewaters will decrease IC detectability of transition metals in these waters. In this study, a CIC method was developed for the analysis of transition metal ions (Fe 3+ , Cu 2+ , Ni 2+ , Zn 2+ , Co 2+ , Mn 2+ , and Fe 2+ ) in fossil fuel associated wastewaters such as Appalachian CMD and O&G wastewaters from the Permian and Bakken shale basins in the United States. CIC system incorporated an on-line chelator column (e.g., the MetPac CC-1) with high selectivity for transition metals over alkali and alkaline earth metals for salt matrix removal prior to transition metal separation and detection. Additional method developments also included acidifying all samples to 2% v/v HCl and using gradient elution rather than isocratic. The recoverability of transition metals in simple salt solutions commonly found in CMD and brine samples (e.g. NaCl, Na 2 SO 4 , CaCl 2 ) using CIC was evaluated and compared to that using traditional IC. Our results found that the CIC system significantly improved transition metal recoveries for samples in 10,000 mg/L CaCl 2 matrix, reaching 87%-108% recovery for all analytes, as opposed to 2-323% recovery in traditional IC. The limits of detection in this study achieved 10.09 – 161.2 μg/L, comparable to reported values in similar IC studies. The developed method was also verified with certified water samples, resulting in 89% - 111% recoveries in samples with higher analyte concentrations (i.e. >4x the LoDs). The developed method achieved 87% -112% recoveries for most analytes in CMD samples and 72%-138% recoveries for Bakken shale samples, relative to ICP-MS values. Overall, the current IC method can be a very good screening tool for fast and cheap analysis for transition metals at mg/L level, to facilitate selection of samples for more detailed ICP-MS analysis.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Comparative techno-economic assessment of osmotically-assisted reverse osmosis and batch-operated vacuum-air-gap membrane distillation for high-salinity water desalination

New developments in pressure- and thermally driven membrane desalination technologies offer the potential to cost-effectively treat high-salinity waters, especially when powered by low-cost solar electricity and thermal energy. This paper presents a comparative techno-economic assessment of the state-of-the-art most promising pressure- and thermally driven membrane technologies for high recovery desalination, namely, osmotically-assisted reverse osmosis (OARO) and batch-operated vacuum-air-gap membrane distillation (batch V-AGMD), to produce potable water while concentrating brine within the range of 140–290 g/L TDS for minimum-liquid-discharge (MLD) and zero-liquid-discharge (ZLD) applications. It is shown that both OARO and batch V-AGMD can treat feedwater and brines with TDS in the range of 30–125 g/L with corresponding fresh water recovery rates of 85–25%. When low cost solar electricity and thermal energy are used, the resulting levelized cost of water (LCOW) from OARO is in the range of 0.70–6.28 $/m 3 , and that from batch-V-AGMD is in the range of 1.74–2.77 $/m 3 . OARO is more cost-effective than batch V-AGMD when feedwater salinity is below 70 g/L and recovery below 75%, whereas batch V-AGMD is more cost-effective at higher recovery rates and salinity levels. Finally, the sensitivity of this comparison on energy prices and module costs is discussed.

42 ENGINEERING↗
DOE OSTI2025

Laboratory evaluation of cyclic underground hydrogen storage in the Temblor sandstone of the San Joaquin Basin, California

Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.

08 HYDROGEN↗
DOE OSTI2020

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗
DOE OSTI2022

Analysis of Brackish Water Desalination for Municipal Uses: Case Studies on Challenges and Opportunities

Brackish water is a widely available, nontraditional water resource that can augment limited freshwater supplies. Although brackish water desalination has been continuously implemented in the United States and worldwide, it is necessary to reduce further its energy consumption, costs, and environmental impacts. This study conducted technoeconomic analyses to evaluate the current desalination and brine management technologies, focusing on the key factors and opportunities for sustainable brackish water desalination for municipal uses. Three case studies were selected as baseline representative of different geographic and operational conditions, including water quality, plant size, pretreatment, desalination, and concentrate management. Here, the technoeconomic analyses and model simulations identified challenges, opportunities, and research priorities to achieve specific pathways for enhanced brackish water desalination regarding levelized costs of water, electricity intensity, water recovery, zero liquid discharge, and brine valorization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Surfactant-Aided Low-Salinity Waterflooding for Low-Temperature Carbonate Reservoirs

Carbonate reservoirs tend to be oil-wet/mixed-wet and heterogeneous because of mineralogy and diagenesis. The objective of this study is to improve oil recovery in low-temperature dolomite reservoirs using low-salinity and surfactant-aided spontaneous imbibition. The low-salinity brine composition was optimized using ζ-potential measurements, contact-angle (CA) experiments, and a novel wettability-alteration measure. Significant wettability alteration was observed on dolomite rocks at a salinity of 2,500 ppm. We evaluated 37 surfactants by performing CA, interfacial-tension (IFT), and spontaneous-imbibition experiments. Here, three (quaternary ammonium) cationic and one (sulfonate) anionic surfactants showed significant wettability alteration and produced 43–63% of original oil in place (OOIP) by spontaneous imbibition. At a low temperature (35°C), oil recovery by low-salinity effect is small compared with that by wettability-altering surfactants. Coreflood tests were performed with a selected low-salinity cationic surfactant solution. A novel coreflood was proposed that modeled heterogeneity and dynamic imbibition into low-permeability regions. The results of the “heterogeneous” coreflood were consistent with that of spontaneous-imbibition tests. These experiments demonstrated that a combination of low-salinity brine and surfactants can make originally oil-wet dolomite rocks more water-wet and improve oil recovery from regions bypassed by waterflood at a low temperature of 35°C.

02 PETROLEUM↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗