Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relative permeability”

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.

395 records · Page 22

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

Effects of ingested nanocellulose on intestinal microbiota and homeostasis in Wistar Han rats

Micron scale crystalline and fibrillar cellulose materials are “generally regarded as safe” (GRAS) as binders and thickeners in food products. Yet, partly due to a lack of relevant toxicological data, nanocellulose (NC) materials, which have unique properties that can be exploited to improve food quality and safety, have yet to receive FDA approval as food ingredients. Results of in vitro and in vivo toxicological studies of ingested NC, detailed in a recent companion report, revealed minimal acute in vitro cytotoxicity, and no toxicity in a subacute rat gavage model, suggesting that NC materials are non-hazardous. However, ingested materials may also modulate gut microbial populations, or alter aspects of intestinal function not evaluated in standard totoxicity testing, which could have important health implications. Here, we report the results of studies of the effects of ingested cellulose nanofibrils (CNF) on the fecal microbiome and metabolome, intestinal (ileal) epithelial cell mRNA expression of cell junction genes, and ileal cytokine production. Feces, plasma, and ileum samples were collected from Wistar Han rats before and after five weeks of biweekly gavages of water or cream, with or without 1% w/w CNF. Analysis of fecal microbial populations revealed that CNF caused changes in both genus and species diversity, with specific effects on species that produce short chain fatty acids, and that have been associated with increased IgA production and elevation of insulin levels. Fecal metabolomic analysis revealed relatively few effects of CNF, with significant changes in the levels of only ten metabolites out of 366 measured. Exposure to CNF also caused differential regulation of mRNA expression of several genes involved in epithelial cell junctions, and increased production of cytokines that both promote and inhibit proliferation of CD8 T cells. These perturbations in intestinal homeostasis contrast somewhat with the absence of in vitro and in vivo toxicity observed previously but would appear to represent minor effects. Additionally, the gut microbiome is highly sensitive to ingested substances, and such disturbances of the intestinal microbial ecosystem do not necessarily represent or predict meaningful pathology. Further studies, including chronic feeding studies, are needed to assess the real health implications, if any, of these changes.

59 BASIC BIOLOGICAL SCIENCES↗

Establishing a silica gel zone in well annulus and evaluating its performance in blocking vertical water flow

Wells are often constructed for monitoring purposes with relatively long screen lengths (e.g., >10 m). Vertical water flows can occur within the artificial or natural filterpack annulus that surrounds the screened interval, bypassing packer assemblies installed inside the wellbore. Attempts to isolate discrete vertical zones during groundwater sampling are unsuccessful when annular flow occurs and lead to remedy decisions based on biased or incorrect interpretations. Blocking vertical annular water flow and contaminant transport will help obtain more accurate concentrations of contaminants from sampling in targeted depth intervals. The application of silica gels formed from the injected colloidal silica CS suspensions is a novel approach to minimize or prevent movement of vertical movement of groundwater in the surrounding filterpack annulus. In this work, we tested the feasibility of injecting CS suspensions to target locations and developed a modified CS formulation that is injectable and prevents gravity sinking. We studied the distribution and penetration of silica gel at laboratory scale in mock well annulus with surrounding formations. We evaluated the performance of the silica gel in blocking vertical water flow in the annulus and in minimizing chemical transport through the gel zone. CS suspension formulations have been defined that are ready for injection, stay in target locations, and form gel within desired time frames. Injection of CS suspensions achieved uniform distribution in a well annulus filter pack, fully occupied the annulus pore space, and penetrated the formation surrounding the filter packer with a sufficient distance to create a hydraulic annular seal when the injection was applied at a sufficient rate. The depth of penetration into the formation was dependent on the permeability contrast between the filter pack and the surrounding formation. Silica gel that formed in the annulus blocked vertical water flow and stopped the chemical transport through the gel zone. In conclusion, this research reveals that using CS suspension injection and sequential gelation (CS-GEL) is a promising technology for blocking vertical water flow and chemical transport through the filter pack in targeted zones within the annulus of long-screened well systems.

Colloidal silica suspension↗

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps

In this project, we developed a techno-economic analysis (TEA) model that can be used to evaluate the viability of a proposed Geological Thermal Energy Storage (GeoTES) design. This MATLAB-based model integrates distinct subsystem models for the reservoir, wells, power cycle, and solar field to capture their distinct characteristics. It applies this approach in simulating GeoTES storage and dispatch operations for durations ranging from hourly to seasonal. Using cases studies based on GeoTES designs provided by industry partners - Premier Resource Management (PRM) and EarthBridge Energy - we validated the TEA model estimations of system performance and costs (such as thermal and electrical power/energy inflow and outflow, capital costs, and levelized costs of energy and storage) for both concentrating solar thermal (CST) and Carnot Battery (CB) pairings with GeoTES (CST-GeoTES and CB-GeoTES). For the CST-GeoTES case, the model was validated against the proposed system designed by PRM. It showed good agreement with PRM's estimations when well and pump costs derived from PRM's estimations were used. When GETEM-based costs were used, there was a slight overprediction due to GETEM's project/site agnostic assumption of these costs. From a sensitivity analysis perspective, the levelized cost of electricity (LCOE) of the CST-GeoTES case was most sensitive to well flow rate and the charging temperature. An optimal design scenario resulted in an LCOE of 0.11 $\$$/kWhe. CST-GeoTES can also provide a source of heat to meet seasonal demands. With 12-hour and 24-hour levelized cost of heat (LCOH) of 0.018 $\$$/kWhth and 0.022 $\$$/kWhth, respectively, CST-GeoTES could be competitive in the California market with an average industrial price of natural gas in California between 0.041-0.047 $\$$/kWhth. The levelized cost of storage (LCOS) for CST-GeoTES depends on the energy storage duration. Although the LCOS is relatively higher for shorter durations (e.g., ~0.50 $\$$/kWhe for 1 hour of storage), it is an order of magnitude lower (0.06 $\$$/kWhe) for longer storage durations and competitive with lithium-ion batteries (beyond 12 hours of storage) and molten-salt thermal energy storage (beyond 32 hours). Energy. Three options were explored and applied to the EarthBridge case study: (1) A Carnot Battery design using R125 working fluid with both hot and cold storage; (2) A Carnot Battery design using R125 working fluid with only hot storage; (3) A Carnot Battery using a commercially available heat pump with carbon dioxide (CO2) working fluid and hot storage only. The CB-GeoTES with cold storage only had a slight (round-trip) efficiency advantage over the system without (43.4% vs. 42.8%). This is because the cold storage is limited by the freezing point of water, so the cold storage is not much colder than the environment. The system using commercially available technologies was the least efficient - partly because different cycles were used in the heat pump (CO2) and heat engine (binary cycle) which leads to some inefficiencies. Using the commercially available design, the levelized cost of energy (LCOS) from the model (0.10 $\$$/kWhe) was higher than that estimated by EarthBridge (0.068 $\$$/kWhe). This is because of the low round-trip (38.7%) efficiency of the commercially available design. Sensitivity analysis reveals that the model is most sensitive to electricity price. Including electricity price in the TEA for CB-GeoTES leads to an increase in LCOS from the base value to 0.25 $\$$/kWhe. To determine storage sites suitable for GeoTES, we gathered and analyzed geological, petrophysical, and geophysical data of oil and gas reservoir and aquifers in California and Texas. We down-selected possible sites based on cut-off values for site characteristics (e.g., reservoir temperature, formation thickness, permeability, porosity, depth, and brine salinity) and preliminary costs. Using this approach, the Carrizo-Wilcox, Yegua-Jackson, and Dockum brackish aquifers in Texas were identified as having the highest suitability. Similarly, in the central California region, the White Wolf, Belridge South Tulare, and Belridge South Reef Ridge were the most suitable. Going further, we assessed the storage potential in the selected sites. To do this we developed distributions of reservoir characteristic data and applied a Monte Carlo-based analysis to account for intrinsic uncertainty in the acquired data. The analysis revealed that the Carrizo-Wilcox aquifer had the highest storage potential with a mean capacity of 554 TWhth (i.e., 63 TWhe). The estimated capacity serves as an upper limit of storage potential given that not all fields in the basin will be developed. We participated in multiple outreach activities including conference presentations, panel session discussions, and the facilitation of a GeoTES workshop at the NREL Golden campus.

15 GEOTHERMAL ENERGY↗

Highly Insulated Wall Systems with Exterior Insulation of Polyisocyanurate under Different Facer Materials: Material Characterization and Long-Term Hygrothermal Performance Assessment

The application of exterior insulation in both new construction and retrofits is a common practice to enhance the energy efficiency of buildings. In addition to increased thermal performance, the rigid insulation can serve to keep the sheathing board warm and serve as a water-resistive barrier to keep moisture-related problems due to condensation and wind-driven rain. Polyisocyanurate (PIR) rigid boards have a higher thermal resistance in comparison to other commonly used exterior insulation boards. However, because of its perceived lower permeance, its use as exterior insulation is not very common. In this study, the hygrothermal property of PIR boards with different facer types and thicknesses is characterized. The material data obtained through experimental test and extrapolation is used in a long term hygrothermal performance assessment of a wood frame wall with PIR boards as exterior insulation. Results show that PIR with no facer has the smallest accumulated moisture on the sheathing board in comparison to other insulation boards. Walls with a bigger thickness of exterior insulation perform better when no vapor barrier is used. The PIR exterior insulation supports the moisture control strategy well in colder climates in perfect wall scenarios, where there is no air leakage and moisture intrusion. In cases where there is trapped moisture, the sheathing board has a higher moisture content with PIR boards with both aluminum or fiberglass type facers. An innovative facer material development for PIR boards can help efforts targeting improved energy-efficient and durable wall systems.

36 MATERIALS SCIENCE↗

Associations between regional blood-brain barrier permeability, aging, and Alzheimer’s disease biomarkers in cognitively normal older adults

Background Increased blood-brain barrier permeability (BBBp) has been hypothesized as a feature of aging that may lead to the development of Alzheimer’s disease (AD). We sought to identify the brain regions most vulnerable to greater BBBp during aging and examine their regional relationship with neuroimaging biomarkers of AD. Methods We studied 31 cognitively normal older adults (OA) and 10 young adults (YA) from the Berkeley Aging Cohort Study (BACS). Both OA and YA received dynamic contrast-enhanced MRI (DCE-MRI) to quantify K trans values, as a measure of BBBp, in 37 brain regions across the cortex. The OA also received Pittsburgh compound B (PiB)-PET to create distribution volume ratios (DVR) images and flortaucipir (FTP)- PET to create partial volume corrected standardized uptake volume ratios (SUVR) images. Repeated measures ANOVA assessed the brain regions where OA showed greater BBBp than YA. In OA, K trans values were compared based on sex, Aβ positivity status, and APOE4carrier status within a composite region across the areas susceptible to aging. We used linear models and sparse canonical correlation analysis (SCCA) to examine the relationship between K trans and AD biomarkers. Results OA showed greater BBBp than YA predominately in the temporal lobe, with some involvement of parietal, occipital and frontal lobes. Within an averaged ROI of affected regions, there was no difference in K trans values based on sex or Aβ positivity, but OA who were APOE4carriers had significantly higher K trans values. There was no direct relationship between averaged K trans and global Aβ pathology, but there was a trend for an Ab status by tau interaction on K trans in this region. SCCA showed increased K trans was associated with increased PiB DVR, mainly in temporal and parietal brain regions. There was not a significant relationship between K trans and FTP SUVR. Discussion Our findings indicate that the BBB shows regional vulnerability during normal aging that overlaps considerably with the pattern of AD pathology. Greater BBBp in brain regions affected in aging is related to APOE genotype and may also be related to the pathological accumulation of Aβ.

Science & Technology - Other Topics↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide

The project Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. Sandia Technologies, LLC, and co-investigator Conrad Geoscience Corporation, examined the potential for large-scale, permanent CO2 storage in sedimentary strata within the Newark Rift Basin. The Newark Rift Basin underlies an industrialized, developed region comprising parts of New York, New Jersey, and Pennsylvania. The project characterized and investigated the suitability of Triassic age sedimentary formations for potential geologic CO2 storage. The project team drilled and cored two test wells to define the sedimentary geologic formations underlying the basin and to document or reach basement rock. With this geologic characterization phase, an integration of seismic, geologic, borehole, and formation core results provided a higher resolution assessment of CO2 storage potential. The Stockton Formation is known to be a potentially favorable geologic storage formation in the basin. In 2011, the 1-NYSTA Tandem Lot stratigraphic test well was drilled to a depth of 6,855 feet in the northern portion of the Newark Basin in southern New York State. Approximately 9 miles south-southeast on the Lamont Doherty Campus, TW-4 was drilled and cored in 2013 to a depth of 1,802 feet and contacted apparent igneous basement at a depth of 1,712 feet. Both wells penetrated the Palisades Sill ranging from 800 feet thick in the eastern well to approximately 1,800 feet in thickness at the 1-NYSTA Tandem Lot deep drill site. A diabase sill can provide an excellent seal and dense confining layer for potential CO2 storage reservoirs and flow layers that are situated beneath it within the Stockton Sandstone. The Stockton Sandstone was encountered beneath the sill in the TW-4 well on the Lamont campus, and data integration suggests that it was likely observed near total depth in the deep 1-NYSTA Tandem Lot well. The test wells confirm and define reservoirs are present beneath the sill and offer CO2 storage potential. The integration of geologic and reservoir characterization of well logs, formation cores, and formation fluids indicated Triassic age lacustrine playa lake and mudbank shales of the Upper Passaic Group can provide an effective seal for the porous and permeable underlying sandstone reservoir layers. This project acquired seismic data, drilled borehole well logs, acquired core samples, and integrated these findings to provide a better understanding of the subsurface geologic formations in the Newark Rift Basin. These findings have contributed to a higher degree of accuracy in predicting potential geologic storage opportunities, while refining geologic storage capacity estimates for the indicated reservoirs and flow units.

.las↗

Metabolic multireactor: Practical considerations for using simple oxygen sensing optodes for high-throughput batch reactor metabolism experiments

We present a system for carrying out small batch reactor oxygen consumption experiments on water and sediment samples for environmental questions. In general, it provides several advantages that can help researchers achieve impactful experiments at relatively low costs and high data quality. In particular, it allows for multiple reactors to be operated and their oxygen concentrations to be measured simultaneously, providing high throughput and high time-resolution data, which can be advantageous. Most existing literature on similar small batch-reactor metabolic studies is limited to either only a few samples, or only a few time points per sample, which can restrict the ability for researchers to learn from their experiments. The oxygen sensing system is based very directly on the work of Larsen, et al. [2011], and similar oxygen sensing technology is widely used in the literature. As such we do not delve deeply into the specifics of the fluorescent dye sensing mechanism. Instead, we focus on practical considerations. We describe the construction and operation of the calibration and experimental systems, and answer many of the questions likely to come up when other researchers choose to build and operate a similar system themselves (questions we ourselves had when we first built the system). In this way, we hope to provide an approachable and easy to use research article that can help other researchers construct and operate a similar system that can be tailored to ask their own research questions, with a minimum of confusion and missteps along the way.

3D printing↗

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Miniaturized Magnetoelastic Sensor System

This article describes the design, assembly, and implementation of a hand-held, magnetic-field-based sensor system that can be adapted for a variety of sensing applications. The miniaturized system is based on Chemical Identification by Magneto-Elastic Sensing (ChIMES) technology, which uses three concentric solenoid coils to wirelessly interrogate a sensor body comprised of a response material coupled to a magnetoelastic wire. The response material expands when it encounters a target, imposing mechanical stress on the wire and altering its magnetic permeability. The sensor bodies are passive, requiring no external power source, and they are small, measuring about 15 mm in length and 3.0 mm in diameter. Up to four sensor bodies can be configured as an evenly-spaced linear array. The sensor system operates by applying a low-frequency, current-stabilized, filtered triangle wave to a uniform-density excitation coil to switch the magnetic domains within the wire. Further, the responses from the sensors are picked up by a detection coil as stress-induced changes in the Faraday voltage, and the strong magnetic field induced by the excitation coil in the detection coil is nullified by a cancellation coil reverse-wound in series with the detection coil. The responses of the sensors in an array are separated in time by a linear gradient dc biasing coil. The sensors can be interrogated through metallic and nonmetallic barriers. The signals from the detection coil and the excitation coil are digitized by a pair of bipolar analog-to-digital converters (ADCs). A Raspberry Pi single-board computer (SBC) and associated software perform data acquisition and control all aspects of the sensor system hardware. The program allows the user to select the number of sensors in the array, the type of signal that is being collected, and the number of samples to take. The program also allows for signal processing of the sensor data, such as baseline correction. The program can differentiate sensor peaks from each other and calculate the magnitude of each sensor response with less than 1% error. The data are then displayed along with a graph of the signal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Porous Membranes for Moisture Probe Protection (SR18024)

In the presence of contaminants, moisture probes can give false moisture readings and have shorter lifetimes. Absolute humidity probes are composed of porous, hygroscopic metallic oxide thin films (Al 2 O 3 ) coupled to electrodes (e.g. Au and Al). While these materials provide absolute humidity readings at low water levels at a relatively low cost, they are sensitive to the presence of contaminants, especially ammonia, due to its high hygroscopicity and ability to bind to the probe. Without frequent recalibration, especially in the presence of contaminants, the humidity levels are inaccurately reported due to measurement drift that occurs when the pores close, changing the impedance of the probe. The goal of this project was to improve the stability and reliability of these moisture probes by incorporating a molecular sieve into the probe to allow only the water to reach the sensor. Molecular sieves contain pores that allow molecule specific permeation. Two dimensional inorganic nanosheets (e.g. BN, MoO 2 , MoS 2 ) have recently been shown to behave as molecular sieves for water. Compared to zeolites, nanoporous materials have higher flux rates through the pores, leading to faster water permeation. Graphene oxide sheets have also been successfully employed for unimpeded water permeation; however, at lower humidity levels, the structure is unstable and the pores shrink, preventing the sieving of molecules. Compared to graphene oxide, sheets based off of inorganic materials, such as molybdenum, have been found to be more stable at lower humidity levels and have faster permeation rates. By creating a barrier in before the Al 2 O 3 sensor chip to prevent the permeation of contaminants into the probe, the lifetime of the probe will be extended. While the molybdenum sheets are a promising material for this application, their performance against ammonia and tritiated compounds have not yet been evaluated. This technology can be applied to other applications, including gas and liquid purification for environmental remediation, single molecule sensors, field effect transistors, and catalysts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS↗

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗