A new constitutive model for underground chemical explosions in partially saturated porous rock
Explore the source record for details and available documents.
Vorobiev, Oleg Y. [External Affiliation]↗
SEARCH · Search NASA
Results for “SATURATION”
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.
Performance of Convection-Permitting and Convection-Parameterized Models in Reproducing the Extreme Precipitation Intensity Relationship with Surface Conditions
Here, this study investigates the warm-season extreme precipitation–temperature scaling relationship in CONUS404, a convection-permitting (4 km) Weather Research and Forecasting (WRF) Model simulation over the conterminous United States for the past four decades, and compares it with the WRF-Thermodynamic Global Warming (WRF-TGW) historical simulation at a coarser resolution (12 km) using parameterized convection. We also analyze the NCEP stage IV and NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) datasets as observational benchmarks. We examine how extreme precipitation intensity (EPI) varies with temperature and saturation deficit over representative regions based on hourly data. The stage IV and IMERG data show a similar pattern of EPI variation with temperature and saturation deficit, except that the EPI peak is lower in IMERG than in stage IV. Under dry and hot conditions, EPI decreases too rapidly with elevated saturation deficit in both CONUS404 and WRF-TGW compared to observations, but the performance of CONUS404 is superior to WRF-TGW. When the near-surface atmosphere is saturated or close to saturated, both CONUS404 and WRF-TGW produce higher peak values of EPI relative to the observational references; IMERG exhibits scaling rates close to the Clausius–Clapeyron (C–C) relationship, while CONUS404, WRF-TGW, and stage IV all demonstrate super-C–C scaling behaviors. Despite marked warming over the past four decades, in both CONUS404 and WRF-TGW, the scaling relationship between EPI and temperature in a saturated atmosphere remains stable and robust. This indicates a strong potential for the EPI–temperature scaling rate under saturation to be used as an emergent constraint in reducing uncertainties of future extreme precipitation projection.
Atmosphere↗
Multi-stage heat pump with two-phase injection for cold climate applications implementing natural refrigerants
Heat pump adoption in cold climates is limited by various technical challenges, including heat capacity degradation, compressor deterioration, and frost accumulation on the outdoor heat exchanger. Furthermore, wide-spread adoption of heat pumps with high global warming potential presents a possible environmental concern. Thus, we present a heat pump concept termed the “saturation cycle,” which uses numerous compression stages to improve energy efficiency and injects two-phase refrigerant between compression stages to dramatically reduce compressor discharge temperatures. Given the growing interest in heat pumps, there is a need for a straightforward framework for parametric evaluation of natural refrigerants in high-efficiency heat pump cycle architectures, like the saturation cycle, in adverse climates. This work develops a thermodynamic modeling framework and focuses on the application of natural refrigerants to explore the performance trends at extremely cold climates down to −25 °C. This study shows that a two-stage saturation cycle and four-stage saturation cycle can reduce the compressor discharge superheating by up to 51–57 % or 73–80 %, respectively, depending on the refrigerant. Furthermore, at −25 °C, a two-stage R-290 saturation cycle heat pump theoretically can use 25 % less power than a single-stage heat pump or can provide up to 59 % more heating capacity when using the same compressor power as the baseline single-stage heat pump. The exergetic parametric analyses show that two-phase refrigerant injection fundamentally improves the thermodynamic performance of the heat pump, reducing overall exergy destruction by 10 % compared to a vapor-injection heat pump, compared at −25 °C. Thus, the saturation cycle is a promising heat pump cycle architecture for cold climate applications due to its superior energy savings (or capacity improvement) and substantial de-superheating capabilities.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
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↗
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↗
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↗
Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers
Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.
Engineering↗
Informing field-scale CO 2 storage simulations with sandbox experiments: The effect of small-scale heterogeneities
Small-scale heterogeneities can significantly affect the fate of the CO 2 plume and trapping during CO 2 migration. We conducted geologic carbon storage field-scale simulations to investigate the impact of small-scale heterogeneities on plume dynamics and trapping performance. Small-scale heterogeneities have been shown to increase the amount of trapped CO 2 during buoyancy-driven flow. The trapped CO 2 saturation is validated by previous sandbox experimental work during buoyancy-driven flow in realistic heterogeneous domains and is implemented through the critical CO 2 saturation parameter (i.e., the first non-zero value in the drainage CO 2 relative permeability curve). Depending on the type and degree of heterogeneity, various critical CO 2 saturation values are exhibited. Furthermore, we investigated the effect of small-scale heterogeneities when multiple capillary pressure models are employed. This study demonstrates that an increase in critical CO 2 saturation reduces the CO 2 plume size and lateral extent, accompanying an increase in residual trapping and a decrease in solubility trapping. Lastly, we show that independent of the capillary pressure model used, an increase in critical saturation leads to similar CO 2 plume dynamics distribution and trapping performance. These results emphasize the importance of quantifying the effect of small-scale heterogeneity as they affect the large-scale behavior of the CO 2 plume.
58 GEOSCIENCES↗
Mechanisms of electromagnetic field control on mineral scaling in brackish water reverse osmosis: Combined homogenous and heterogeneous nucleation
Electromagnetic field (EMF) treatment has emerged as a promising approach for scaling control due to its cost-effectiveness, simplicity, and low energy consumption. However, there is a limited understanding of the mechanisms by which applied EMF impacts mineral scaling in reverse osmosis (RO) systems. This has led to inconclusive and varied results and uncertainties regarding its effectiveness. This study elucidates the impacts of EMF on homogenous and heterogeneous nucleation and membrane performance during RO desalination of different feedwaters. Our results reveal that EMF exhibits greater efficacy in treating near-saturated water (SI∼0), especially when coupled with extended hydraulic flushing (HF). For saturated brackish water desalination, heterogeneous scaling predominantly occurs on membrane surfaces, with the effectiveness of EMF in inhibiting scaling primarily attributed to the hydration effect. In supersaturated solutions, EMF promotes bulk precipitation due to the magnetohydrodynamic effect, quickly blocking membrane pores. Thus, when the saturation reaches a certain high level during RO desalination, magnetohydrodynamic EMF effects can accelerate flux decline caused by homogeneous scaling. In conclusion, this work provides an efficient method for predicting EMF efficiency, emphasizing the importance of saturation conditions and HF cleaning duration in determining membrane performance, suggesting these show promise for improving undersaturated or near-saturated feedwater desalination via RO.
42 ENGINEERING↗
3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope
Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.
58 GEOSCIENCES↗
Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers
In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.
carbon dioxide (CO2)↗
Greenhouse gas flux response in biochar- and compost-amended urban soils under simulated soil hydrologic dynamics
Understanding greenhouse gas emission dynamics in lawn soils is essential for improving climate change mitigation strategies in urban and suburban environments. This project measured fluxes of carbon dioxide (CO 2 ), methane (CH 4 ), and nitrous oxide (N 2 O) fluxes from turfgrass soil columns amended with biochar, compost, biochar-compost blend, control (no amendment) in a controlled greenhouse mesocosm experiment. This project simulated contrasting water saturation regimes consisting of normal irrigation with sprinkler, transient half or full saturation by water table manipulation, and subsequent drying phase under both sod or seeded grass conditions. Compost-amended columns (compost or biochar-compost blend) exhibited higher average CO 2 fluxes (6 µmol m-2 s-1) compared to biochar or control columns (4.5 – 5.0 µmol m-2 s-1) across all saturation levels under sod conditions. Seeded grass conditions generally resulted in less CO 2 emissions. The CO 2 fluxes were positively correlated with air temperature and negatively correlated with soil moisture. Biochar-amended columns retained high soil moisture (~95 %) throughout experiments, demonstrating superior moisture retention compared to compost. CH 4 and N 2 O fluxes exhibited temporal increases (5-8 nmol m-2 s-1) during saturation and drainage phases, indicating their sensitivity to hydrologic conditions. These findings suggest that temperature and amendment types are primary driver of CO 2 emissions, while CH 4 and N 2 O fluxes are more responsive to water saturation dynamics in lawn soils.
54 ENVIRONMENTAL SCIENCES↗
Laboratory evaluation of cyclic underground hydrogen storage in the Temblor sandstone of the San Joaquin Basin, California
Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.
08 HYDROGEN↗
Time‐And‐Space Averaging Applied to Intermittent Multiphase Flow Experiments
Abstract Various researchers have studied fluctuations in pore‐scale phase occupancy during multiphase flow in porous media using synchrotron‐based X‐ray microcomputed tomography (micro‐CT). However, the impact of these fluctuations on the concept of a representative volume is not yet fully understood. In this study, we performed spatial and temporal averaging of multiphase flow experiments visualized with synchrotron‐based micro‐CT, focusing on oil saturation as the key parameter to determine a representative time‐and‐space average. Our findings revealed that a saturation value representative of both time and space was achieved during fractional flow experiments in drainage mode with fractional flows of 0.8, 0.5, and 0.3. Furthermore, we computed a range of relative permeabilities on the basis of whether momentaneous saturation or time‐and‐space averaged saturation was utilized for direct simulation. Our results highlighted the importance of time‐and‐space averaging in determining a representative relative permeability and indicated that the temporal and spatial scales covered in a typical micro‐CT flow experiment were sufficient to obtain a representative saturation value for sandstone rock under intermittent flow conditions.
Environmental Sciences & Ecology↗
Quantitative analysis of zonal flow influence on turbulent plasmas driven by trapped electrons
The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.
Song, Jiheon [Hanyang University, Seoul (Korea, Re↗
Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration
In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.
Holcomb, Paul↗
Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics
Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.
58 GEOSCIENCES↗
The quantum evolutions of the diffractive transverse-momentum dependent gluon distribution
Using the Colour Glass Condensate description of electron-nucleus collisions at high energy, we study the diffractive production of a pair of jets with transverse momenta much larger than the nuclear saturation momentum Q s . At leading order in the QCD coupling, the di-jet cross-section exhibits transverse-momentum dependent (TMD) factorisation, with a gluon diffractive TMD distribution (DTMD) which is controlled by gluon saturation and describes the transverse-momentum imbalance between the produced jets. The next-to-leading corrections generate the various quantum evolutions of the diffractive gluon distribution. We focus on the Collins-Soper-Sterman (CSS) evolution which describes the change in the gluon DTMD when increasing the “hard scale” (the typical transverse momentum of the di-jets). We consider two different representations for this equation, one in transverse-momentum space, the other one in transverse-coordinate space. They are not fully equivalent with each other (despite being related by a Fourier transform) because of the respective boundary conditions. These conditions encode the essential physics of gluon saturation together with the effects of two other types of quantum evolution: the BK/JIMWLK evolution over the rapidity gap (“inside the Pomeron”) and the DGLAP evolution outside the rapidity gap (“within the diffractive system”). We demonstrate that, due to gluon saturation, one can compute both the boundary conditions and the CSS solutions mostly from first principles, without the need for a non-perturbative Sudakov. We numerically find a good agreement between the CSS solutions in the two aforementioned representations.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗