Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volume visualization”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Verification of a two-dimensional infiltration model for the resin transfer molding process

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin was verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTI IM7/8HS preform and corn oil. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle.

Hammond, Vincent H.↗
NASA NTRS2010

A Goddard Multi-Scale Modeling System with Unified Physics

A multi-scale modeling system with unified physics has been developed at NASA Goddard Space Flight Center (GSFC). The system consists of an MMF, the coupled NASA Goddard finite-volume GCM (fvGCM) and Goddard Cumulus Ensemble model (GCE, a CRM); the state-of-the-art Weather Research and Forecasting model (WRF) and the stand alone GCE. These models can share the same microphysical schemes, radiation (including explicitly calculated cloud optical properties), and surface models that have been developed, improved and tested for different environments. In this talk, I will present: (1) A brief review on GCE model and its applications on the impact of the aerosol on deep precipitation processes, (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), and preliminary results (the comparison with traditional GCMs), and (3) A discussion on the Goddard WRF version (its developments and applications). We are also performing the inline tracer calculation to comprehend the physical processes (i.e., boundary layer and each quadrant in the boundary layer) related to the development and structure of hurricanes and mesoscale convective systems. In addition, high - resolution (spatial. 2km, and temporal, I minute) visualization showing the model results will be presented.

Tao, Wei-Kuo↗
NASA NTRS1998

Lessons Learned During Implementation and Early Operations of the DS1 Beacon Monitor Experiment

A new approach to mission operations will be flight validated on NASA's New Millennium Program Deep Space One (DS1) mission which launched in October 1998. The Beacon Monitor Operations Technology is aimed at decreasing the total volume of downlinked engineering telemetry by reducing the frequency of downlink and the volume of data received per pass. Cost savings are achieved by reducing the amount of routine telemetry processing and analysis performed by ground staff. The technology is required for upcoming NASA missions to Pluto, Europa, and possibly some other missions. With beacon monitoring, the spacecraft will assess its own health and will transmit one of four beacon messages each representing a unique frequency tone to inform the ground how urgent it is to track the spacecraft for telemetry. If all conditions are nominal, the tone provides periodic assurance to ground personnel that the mission is proceeding as planned without having to receive and analyze downlinked telemetry. If there is a problem, the tone will indicate that tracking is required and the resulting telemetry will contain a concise summary of what has occurred since the last telemetry pass. The primary components of the technology are a tone monitoring technology, AI-based software for onboard engineering data summarization, and a ground response system. In addition, there is a ground visualization system for telemetry summaries. This paper includes a description of the Beacon monitor concept, the trade-offs with adapting that concept as a technology experiment, the current state of the resulting implementation on DS1, and our lessons learned during the initial checkout phase of the mission. Applicability to future missions is also included.

Sherwood, Rob↗
DOE OSTI2023

Utah FORGE 2439: Well 16B(78)-32 Field-Test Data from Mini-Frac Tests

This submittal includes the field-test data collected during stress tests conducted in the Utah FORGE 16B(78)-32 wellbore to measure/characterize the stresses in the geothermal reservoir. The type of stress test performed is referred to as a mini-frac test or a micro-frac test. The test is a hydraulic fracture test that involves injecting a small volume of water into a short interval of the reservoir that is isolated by a straddle packer to create a fracture. The test provides information about the minimum and maximum horizontal stress at the test depths. A total of seven mini-frac tests were performed within the upper (vertical or nearly vertical) section of the 9-5/8-inch diameter uncased wellbore, between depths of 5,202 ft MD [5,202 ft TVD] and 5,980 ft MD [5966.18 ft TVD]. Two types of data were generated as part of the mini-frac tests. The first type of data is time-series data for each mini-frac test (e.g., injection-rate, test-interval pressure), and are contained in the "Mini-Frac Test Data.zip" file below. The second type is geophysical log data, attached below in the "Geophysical Logs.zip" file. Here, two types of image logs and acoustic logs were obtained before and after the mini-frac tests to examine each test depth for visual evidence of induced fracture(s) and to measure the azimuth (orientation) of any observed fractures. The attached Readme file includes the test depth information and contents of each data file included here. Analysis of the field-test data was conducted and is documented in the attached report below.

15 GEOTHERMAL ENERGY↗
NASA NTRS2013

NextGen Operations in a Simulated NY Area Airspace

A human-in-the-loop simulation conducted in the Airspace Operations Laboratory (AOL) at NASA Ames Research Center explored the feasibility of a Next Generation Air Transportation System (NextGen) solution to address airspace and airport capacity limitations in and around the New York metropolitan area. A week-long study explored the feasibility of a new Optimal Profile Descent (OPD) arrival into the airspace as well as a novel application of a Terminal Area Precision Scheduling and Spacing (TAPSS) enhancement to the Traffic Management Advisor (TMA) arrival scheduling tool to coordinate high volume arrival traffic to intersecting runways. In the simulation, four en route sector controllers and four terminal radar approach control (TRACON) controllers managed traffic inbound to Newark International Airport's primary runway, 22L, and its intersecting overflow runway, 11. TAPSS was used to generate independent arrival schedules for each runway and a traffic management coordinator participant adjusted the arrival schedule for each runway 11 aircraft to follow one of the 22L aircraft. TAPSS also provided controller-managed spacing tools (slot markers with speed advisories and timelines) to assist the TRACON controllers in managing the arrivals that were descending on OPDs. Results showed that the tools significantly decreased the occurrence of runway violations (potential go-arounds) when compared with a Baseline condition with no tools. Further, the combined use of the tools with the new OPD produced a peak arrival rate of over 65 aircraft per hour using instrument flight rules (IFR), exceeding the current maximum arrival rate at Newark Liberty International Airport (EWR) of 52 per hour under visual flight rules (VFR). Although the participants rated the workload as relatively low and acceptable both with and without the tools, they rated the tools as reducing their workload further. Safety and coordination were rated by most participants as acceptable in both conditions, although the TRACON Runway Coordinator (TRC) rated neither as acceptable in the Baseline condition. Regarding the role of the TRC, the two TRACON controllers handling the 11 arrivals indicated that the TRC was very much needed in the Baseline condition without tools, but not needed in the condition with tools. This indicates that the tools were providing much of the sequencing and spacing information that the TRC had supplied in the Baseline condition.

NextGen↗
NASA NTRS1994

HST photometry of the trapezium cluster

We have obtained images of 11 fields in the Trapezium cluster with the Planetary Camera (PC) of the Hubble Space Telescope (HST) in order to extend Herbig & Terndrup's (1986) study of this prototype, dense embedded cluster to fainter magnitudes than is possible from the ground. Using these images, we have identified 319 stars within an area of approximately 12 sq arc min corresponding roughly to a volume of approximately 0.065 cu pc assuming the cluster is approximately spherically symmetric. Our completeness limits for star identification in V-band and I-band images are V approximately = 20 and I(sub c) approximately = 19 respectively, corresponding to a mass limit of approximately 0.15 solar mass if the faintest stars have the same average A(sub v) as that estimated for the brighter stars in the cluster. We have compared the V versus V-I color-magnitude diagram derived from the HST photometry to new theoretical isochrones. Star formation in the Trapezium appears to be remarkably coeval, with greater than or = 80% of the stars having inferred ages less than 1 Myr. Over the somewhat limited mass range of the observations, there is no evidence for 'bimodal' star formation-the high- and low-mass stars appear to have the same ages. The sharp cores of the HST images and the small angular size of the PC pixels has allowed us to identify 35 new visual binaries in the cluster with separations from approximately 0.06 arc sec (approximately 26 AU) to approximately 1.0 arc sec (approximately 440 AU). For the range of binary separations that we are sensitive to, the observed binary frequency for the Trapezium is essentially identical to that estimated for field low-mass stars by Duquennoy & Mayor (1991). The most straightfoward inference from this result is that binaries in this separation are unlikely to be formed by a tidal capture process. We have also identified three stars which have associated compact nebulosity visible in the HST images. One of these star + nebulosity cases was previously identified by O'Dell, Wen, & Hu (1993)-these objects appear to form a class of objects whose circumstellar matter is being 'lit up', most likely by Theta(sup 1) Ori C, enabling the gas to be observable at both optical and radio wavelengths (Felli et al. 1993a, b). We provide a brief summary of the optical properties of the other radio sources which appear in our PC images.

Prosser, Charles F.↗
NASA NTRS2014

Effects of Long Duration Spaceflight on Venous and Arterial Compliance

The visual impairment and intracranial pressure syndrome (VIIP) is a newly described space flight-associated medical condition made up of a constellation of symptoms affecting at least 34% of American astronauts who have flown International Space Station (ISS) missions. VIIP is defined primarily by visual acuity deficits and anatomical changes to eye structures, and is thought to be related to elevated intracranial pressure secondary to space flightinduced cephalad fluid shifts. Loss of visual acuity could be a significant threat to crew health and performance and may be suggestive of other adaptations with implications for years post-flight. Our primary objective is to determine whether vascular compliance is altered by space flight and whether such adaptations are related to the incidence of VIIP. In particular, we will measure ocular parameters and vascular compliance in vessels of the head and neck in astronauts who have no space flight experience, in astronauts before, during, and after space flight, and in bed rest subjects with conditions similar to space flight. Additionally, we will analyze astronaut data from the Lifetime Surveillance of Astronaut Health (LSAH) archive to determine which factors might be predictive of the development of VIIP. The project will be conducted in four separate but related parts. To understand the baseline condition of astronauts without any prior space flight experience, we will study 10 astronauts who have never flown in space by performing a comprehensive evaluation of the vasculature of the head, neck and eyes. Hemodynamic data (stroke volume and blood pressure), ocular (tonometry and ocular ultrasound), venous and arterial parameters will be acquired across a range of tilt angles (20, 10, 0, -10, -20 degrees). Vessels to be studied include the temporal, jugular, and vertebral veins and the cerebral, carotid and vertebral arteries. Ophthalmic data from the annual physical will be obtained through data sharing. To examine the relation between vascular compliance in the head and neck and the development of VIIP after a long duration space flight, we will study 10 astronauts before, during, and after long-duration ISS missions. Pre- and post-flight testing will be identical to that described above. During flight, images of the same vessels of interest will be obtained for later analysis. Ophthalmic data including VIIP scores will be obtained through data sharing from medically-required tests. To investigate the effects of age and elevated sodium intake, two potential contributors to VIIP, we will study 24 men (in two age groups: 25-35 and 45-55) during a 14 day 6deg head-down bed rest, a well-accepted analog of space flight. Standard NASA bed rest conditions will be maintained except for dietary sodium. Sodium intake will be similar to that of ISS astronauts, which is higher than consumed in previous bed rest studies. Pre- and post-bed rest testing procedures will be identical to the testing protocol described above for astronauts. Ophthalmic testing (optical coherence tomography, fundoscopy, and tonometry) will be conducted on the same day that vascular compliance measures are obtained. To identify parameters that may relate to an increase in an astronaut's susceptibility to developing VIIP, we will use data mining techniques to evaluate astronaut data obtained from the LSAH. Medical history, family history, space flight history and its related exposures, and history of high performance jet aircraft exposure will be examined for their potential relationship to ocular data. We hypothesize that the cephalad fluid shift induced by space flight will result in structural and functional adaptations in head and neck vessels leading to decreased vascular compliance and related to the development of VIIP symptoms. Further, although VIIP has not been observed in previous bed rest studies, we hypothesize that an elevated sodium intake will increase the incidence of VIIP symptoms in this space flight analog. Finally, we hypothesize that data mining analyses will reveal relationships between health history, previous exposures (including space flight and high performance aircraft), and the development of VIIP in the astronaut population.

Ribeiro, L. C.↗
NASA NTRS1997

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗
NASA NTRS1992

Adaptability of the oxidative capacity of motoneurons

Previous studies have demonstrated that a chronic change in neuronal activation can produce a change in soma oxidative capacity, suggesting that: (i) these 2 variables are directly related in neurons and (ii) ion pumping is an important energy requiring activity of a neuron. Most of these studies, however, have focused on reduced activation levels of sensory systems. In the present study the effect of a chronic increase or decrease in motoneuronal activity on motoneuron oxidative capacity and soma size was studied. In addition, the effect of chronic axotomy was studied as an indicator of whether cytoplasmic volume may also be related to the oxidative capacity of motoneurons. A quantitative histochemical assay for succinate dehydrogenase activity was used as a measure of motoneuron oxidative capacity in experimental models in which chronic electromyography has been used to verify neuronal activity levels. Spinal transection reduced, and spinal isolation virtually eliminated lumbar motoneuron electrical activity. Functional overload of the plantaris by removal of its major synergists was used to chronically increase neural activity of the plantaris motor pool. No change in oxidative capacity or soma size resulted from either a chronic increase or decrease in neuronal activity level. These data indicate that the chronic modulation of ionic transport and neurotransmitter turnover associated with action potentials do not induce compensatory metabolic responses in the metabolic capacity of the soma of lumbar motoneurons. Soma oxidative capacity was reduced in the axotomized motoneurons, suggesting that a combination of axoplasmic transport, intracellular biosynthesis and perhaps neurotransmitter turnover represent the major energy demands on a motoneuron. While soma oxidative capacity may be closely related to neural activity in some neural systems, e.g. visual and auditory, lumbar motoneurons appear to be much less sensitive to modulations in chronic activity levels.

NASA Discipline Musculoskeletal↗
DOE OSTI2026

Data Fusion for the Development of a Multimodal Freight Transload Facilities Dataset in the U.S.

To withstand the growing demand of commodity volume and its strain on the transportation infrastructure, it is necessary to identify the flow of commodities by route and mode. However, a national multimodal freight routing model does not exist for the U.S. The development of such model requires multiple building blocks, such as virtual representations of roadway, railway, and waterway networks, transload facilities (TFs), and access/egress links. Most of these blocks have a robust database in the U.S., except for the TFs. Here, this paper presents the fusion of dispersed and heterogeneous representations of multimodal TFs into a single, comprehensive, geospatial freight TF dataset. The TF dataset is derived from several sources, including the U.S. Army Corps of Engineers Master Docks Plus, the National Transportation Atlas Database, the Intermodal Association of North America, industry publications, and other public information. First, individual datasets were queried and reconciled. A geocoding/reverse geocoding process was applied to get the best street address and latitude/longitude location for each terminal. Then, duplicate terminals were identified by a fuzzy match algorithm based on terminal name and location, and removed. Validation was performed by visual inspection of random facilities. The main contributions of this work are: a publicly available version of the TF dataset, including facility location and multimodal transfer capability of 9,003 facilities, and an enterprise-version with the same facilities but including commodity handling capabilities. The main purpose of developing the TF dataset is to inform multimodal routing algorithms. The proposed TF dataset allows for credibly modeling the multimodal transfer of commodities within shipment routes.

Commodity Routing↗
DOE OSTI2018

A complex approach to assessing porous structure of structured ceramics obtained by SPS technique

Highlights: • A complex approach to microstructural analysis of nanostructured SPS ceramics was presented. • A bunch of gas, intrusion and microscopic methods were applied for open/closed porosity determination in nanoceramics. • Kinetic TPR analysis for structural changes in nanoceramic was carried out. • 3D pore visualization for nanoceramics was provided via SEM/FIB tomography • Porosity calculations were performed via BET, BJH, t-plot, Washburn equation, FEI Avizo 9.1 software results. - Abstract: Spark plasma sintering (SPS) is a highly prospective technology providing rapid consolidation of nano-disperse compositions into structured ceramics with complex structural organization. Characterizing pore structure of such systems remains a challenge due to different types, size and shape of the pores as well as structural defects occurring in the material, thus requiring special methods to quantify certain porosity type. The paper presents a complex approach to studying structure of SPS hematite-based structured ceramics. The method involves a number of gas, intrusion and electron microscopy probing techniques to characterize open and closed pores in the ceramics bulk. Low temperature nitrogen sorption has been implemented to estimate open micro- and mesoporosity in the size range 0.7–50 nm using BET, BJH and t-plot methods to process experimental data. Additionally, structural peculiarities of consolidated structured ceramics have been studied using temperature-programmed reduction (TPR) with the kinetic analysis of the results being carried out. Volume of macropores sized above 50 nm has been measured with mercury intrusion. Closed pores and structural defects have been quantified with SEM and FIB/SEM tomography. By comparing experimental and estimated results we have shown that addressing all porosity types provides a highly precise structural characterization of the material.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Defining a Metric Space of Host Logs and Operational Use Cases

Host logs, in particular, Windows Event Logs, are a valuable source of information often collected by security operation centers (SOCs). The semi-structured nature of host logs inhibits automated analytics, and while manual analysis is common, the sheer volume makes manual inspection of all logs impossible. Although many powerful algorithms for analyzing time-series and sequential data exist, utilization of such algorithms for most cyber security applications is either infeasible or requires tailored, research-intensive preparations. In particular, basic mathematic and algorithmic developments for providing a generalized, meaningful similarity metric on system logs is needed to bridge the gap between many existing sequential data mining methods and this currently available but under-utilized data source. In this paper, we provide a rigorous definition of a metric product space on Windows Event Logs, providing an embedding that allows for the application of established machine learning and time-series analysis methods. We then demonstrate the utility and flexibility of this embedding with multiple use-cases on real data: (1) comparing known infected to new host log streams for attack detection and forensics, (2) collapsing similar streams of logs into semantically-meaningful groups (by user, by role), thereby reducing the quantity of data but not the content, (3) clustering logs as well as short sequences of logs to identify and visualize user behaviors and background processes over time. Overall, we provide a metric space framework for general host logs and log sequences that respects semantic similarity and facilitates a wide variety of data science analytics to these logs without data-specific preparations for each.

Verma, Miki↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
NASA NTRS1993

Preliminary design of an intermittent smoke flow visualization system

A prototype intermittent flow visualization system that was designed to study vortex flow field dynamics has been constructed and tested through its ground test phase. It produces discrete pulses of dense white smoke consisting of particles of terephthalic acid by the pulsing action of a fast-acting three-way valve. The trajectories of the smoke pulses can be tracked by a video imaging system without intruding in the flow around in flight. Two methods of pulsing the smoke were examined. The simplest and safest approach is to simply divert the smoke between the two outlet ports on the valve; this approach should be particularly effective if it were desired to inject smoke at two locations during the same test event. The second approach involves closing off one of the outlet ports to momentarily block the flow. The second approach requires careful control of valve dwell times to avoid excessive pressure buildup within the cartridge container. This method also increases the velocity of the smoke injected into the flow. The flow of the smoke has been blocked for periods ranging from 30 to 80 milliseconds, depending on the system volume and the length of time the valve is allowed to remain open between valve closings.

Ward, Donald T.↗
NASA NTRS2022

In Situ Measurements of Surface Texture with Virtual Environments Support Science-Driven Human Surface Operations on the Moon and Beyond

Visualization tools enabling real-time scientific analysis are important for supporting future astronaut operations on the lunar surface. Such tools can be built into virtual environments to support scientific investigations, as well as situational awareness, real-time decision making, and efficient communication between astronauts and ground and support systems. Understanding how these tools can be optimized for science is essential for upcoming Artemis missions. In this contribution, we discuss how measurements of surface texture at multiple length scales can greatly enhance in situ science on/of the Moon, and eventually Mars, asteroids, and beyond. Roughness measurements at various wavelengths directly support objectives defined in the Artemis Science Plan, including (O1) “understanding planetary processes,” (O2) “understanding volatile cycles,” and (O3) “interpreting the impact history of the Earth-Moon system” . Key scientific analyses enabled by texture measurements at different length scales include: ● Sub-centimeter scales: Texture measurements can help constrain lava flow crystallinity, lava rheology, emplacement flow dynamics, and cooling histories (O1). Measurements of lacunarity (voids in fractal fill space) can shed light on eruptive volatile content, residence time of migrating volatiles, and near-surface volume available for micro-cold trapping of volatiles (O1, O2). ● Centimeter–meter scales: Texture measurements can be used for the differentiation of individual lava flows, the reconstruction of local stratigraphies and emplacement sequences, characterization of post-emplacement surface modification processes (O1, O3). Derived roughness (polarization) metrics can be used in the detection of water ice and characterization of ice properties (e.g., purity, grade, depth, abundance). ● Hectometer–Kilometer scales: Texture measurements can be used to differentiate major geologic surface units and surface structures (O1), constrain the presence of abundant ground ices (O2), and analyze surface modification and estimate surface age (O3). Real-time measurements of surface texture across these multiple length scales will enable efficient sample identification and scientific investigations by future astronauts. To support these investigations and the objective classification of surface texture, virtual environments employed by astronauts should be able to instantaneously convert raw data into processed data (e.g., digital terrain and elevation models) and derived metrics (e.g., RMS, std, Hurst, CPR) and perform statistical analyses (e.g., PCA, outliers, correlation matrices). Such tools are being developed and tested by the Resource Exploration and Science of our Cosmic Environment (RESOURCE) team, a node of NASA’s Solar System Exploration Research Virtual Institute (SSERVI), and are an excellent example of the powerful synergies of human and robotic ground assets critical in the return of humans to the Moon.

Ariel N. Deutsch↗
DOE OSTI2019

Subtask 1.1 – Advanced Characterization of Unconventional Oil and Gas Reservoirs to Enhance CO2 Storage Resource Estimates

The Energy & Environmental Research Center (EERC) developed methods and approaches for advanced characterization of unconventional oil and gas reservoirs to enhance CO2 storage resource estimates using samples collected from organic-rich shales of the Bakken Formation. The project was funded under U.S. Department of Energy (DOE) Cooperative Agreement No. DE-FE0024233, with nonfederal cost share provided by Hitachi High Technologies America. The project participants included the EERC, Hitachi High Technologies America, the team at the National Energy Technology Laboratory’s (NETL’s) computed tomography (CT) scanning laboratory in Morgantown, West Virginia, and several NETL staff who have worked on CO2 storage resource estimation. Investigations on the use of subsurface geologic formations for CO2 storage have been ongoing for many years. In the last decade, improvements in directional drilling technology and hydraulic fracturing of organic-rich shale and other tight rock formations have opened new opportunities for CO2 storage, with the potential cobenefit of enhancing gas and oil recovery. Gas-bearing, organic-rich shales have been of particular interest as a carbon storage target because their ability to store large volumes of CO2 is implied by their ability to store methane; however, estimating the CO2 storage resource of organic-rich shales and other unconventional reservoirs is challenging because of the lack of detailed geologic and petrophysical data that are needed to develop improved volumetric equations. According to NETL’s method for estimating regional-scale CO2 storage resources, the key factors that control CO2 storage in geologic formations include void space within induced fractures, natural fractures and matrix pores, as well as the occurrence and distribution of organic matter (OM) and clay minerals. A key goal of this effort was to use field emission scanning electron microscopy (FESEM) analysis of fine-grain unconventional reservoir rocks coupled with advanced image analysis and machine learning techniques to quantify the occurrence and distribution of porosity, OM, and clay minerals. By better understanding the potential migration pathways of CO2 within tight rocks and quantifying the occurrence of sorptive components (OM and clays) along potential flow pathways, methods to estimate the CO2 storage resource of such rocks can be refined. The Bakken petroleum system (Bakken) was used as the source of rock samples for this effort because it includes organic-rich shales with contrasting degrees of thermal maturity as well as tight, nonshale clastic rocks with varying clay mineral content. To better understand the key components of interest in tight rock samples, the EERC worked closely with Hitachi to optimize the FESEM settings necessary to resolve the sample mineralogy, porosity and OM content at micro- to nanoscale resolution. The EERC also worked closely with Hitachi to apply and improve its advanced mineral identification and characterization system (AMICS) software so that it is better suited for mineral identification at the high resolution necessary to resolve features of interest in fine-grained rocks such as organic-rich shales. The EERC also applied advanced image analysis software packages, including Ilastik, to classify and quantify the key components of interest within the samples using FESEM imagery. These advanced techniques allowed for quantification and differentiation of the void space within fractures, pores and OM, as well as a technique to determine and quantify the key mineral and OM components surrounding pore spaces. Being able to better quantify the void space and the potential accessibility of CO2 to sorptive minerals or OM within tight rocks allows for improved estimation of the CO2 storage resource potential of organic-rich shales and other unconventional formations. The information derived from the image analysis techniques was used to refine NETL’s method for estimating regional-scale CO2 storage resources in collaboration with staff from NETL. An additional goal of this work was to better understand rates of CO2 permeation, flow pathways, and sorption in organic-rich shales and nonshale reservoir rocks and to evaluate potential changes in CO2 flow over time. The EERC conducted long-duration (several weeks) CO2 flow-through testing at reservoir conditions using Bakken shale and nonshale samples to assess CO2 permeation rates and retention/sorption within core plugs to represent bulk matrix behavior. Conventional isothermal adsorption testing over a range of pressure conditions was also performed on pulverized Bakken samples to better assess potential maximum sorptive capacity of the samples. The results of this testing illustrated the higher sorptive capacity of the Bakken shales when compared to the sorptive capacity of the nonshale reservoir rock, which contains very little OM when compared to the Bakken shales. The results also illustrate the differences in permeation behavior as a function of fracture flow versus matrix flow. In an attempt to visualize CO2 permeation and flow behavior in Bakken shale samples over time, the EERC collaborated with NETL’s CT Lab in Morgantown, West Virginia. A variety of tests were performed at the NETL lab on unfractured and fractured samples of the Bakken shale. The unfractured Bakken shale samples did not result in CO2 migration within the samples over the testing duration; thus imagery of CO2 migration pathways within the matrix of Bakken samples was not obtained. The results of the flow-through testing within fractured samples illustrate changes in CO2 permeation rates over time, with a trend of decreasing CO2 flow over time. The results of this effort suggest that FESEM analysis coupled with advanced image analysis and machine learning techniques are effective in characterizing and quantifying the mineralogy, porosity, fractures, and OM content of 2-D imagery at the high resolution needed to identify features of interest in tight rocks. In addition, the image analysis techniques allow for identification and quantification of different porosity types as well as the key mineral components and/or OM that CO2 would contact in the reservoir as it migrates through pores and fracture networks. This is key information needed to better estimate the CO2 storage resource potential of organic-rich shales and other tight reservoirs. Additional work is needed to determine the sorptive capacity of the individual sample components (OM and clay minerals) at reservoir conditions so that this approach and the refined CO2 storage resource equation can be applied to other shale and unconventional reservoirs. In addition, subsequent work should be performed to better understand the number of images needed to use stereological approaches to estimate the volume of sample components using 2-D imagery, to explore machine learning techniques to differentiate between induced vs. natural fractures, and to estimate the connectivity of pores and fracture networks using 2-D imagery. Being able to derive this information at the microscale would facilitate scaling up of the sample characterization data used for reservoir modeling and simulation to estimate the CO2 storage resource potential of tight formations.

03 NATURAL GAS↗