SEARCH · Engineering Papers

Results for “sequestration”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Machine Learning-Driven Quantification of CO2 Plume Dynamics at Illinois Basin Decatur Project Sites Using Microseismic Data

This study utilizes machine learning to quantify CO2 plume extents by analyzing microseismic data from the Illinois Basin Decatur Project (IBDP). Leveraging a unique dataset of well logs, microseismic records, and CO2 injection metrics, this work aims to predict the temporal evolution of subsurface CO2 saturation plumes. The findings illustrate that machine learning can predict plume dynamics, revealing vertical clustering of microseismic events over distinct time periods within certain proximities to the injection well, consistent with an invasion percolation model. The buoyant CO2 plume partially trapped within sandstone intervals periodically breaches localized barriers or baffles, which act as leaky seals and impede vertical migration until buoyancy overcomes gravity and capillary forces, leading to breakthroughs along vertical zones of weakness. Between different unsupervised clustering techniques, K-Means and DBSCAN were applied and analyzed in detail, where K-means outperformed DBSCAN in this specific study by indicating the combination of the highest Silhouette Score and the lowest Davies–Bouldin Index. The predictive capability of machine learning models in quantifying CO2 saturation plume extension is significant for real-time monitoring and management of CO2 sequestration sites. The models exhibit high accuracy, validated against physical models and injection data from the IBDP, reinforcing the viability of CO2 geological sequestration as a climate change mitigation strategy and enhancing advanced tools for safe management of these operations.

Iyegbekedo, Ikponmwosa↗

DOE OSTI2025

Effect of CO 2 Concentration on the Microbial Activity of Orenia metallireducens (Strain Z6) in Surface Inert Materials

Carbon dioxide (CO 2 ) sequestration has garnered widespread attention as a key strategy for mitigating CO 2 emissions and combating the greenhouse effect. However, the mechanisms underlying the interactions between CO 2 , widespread siliceous minerals and biological processes remain unclear. The present study explored the potential impacts of different CO 2 concentrations on microbial activity, environmental conditions and their feedback on the fate of CO 2 . A total of 20 experimental conditions was created, with the variables including different natural and synthetic siliceous minerals (e.g., quartz sand and a type of commercial glass beads), the presence or absence of the iron-reducing microorganism Orenia metallireducens (strain Z6) and varying CO 2 concentrations (0%, 20%, 50%, 100%) in the presence of ferrihydrite and pyruvate. Geochemical, microbial and mineralogical analyses revealed that elevated CO 2 concentrations significantly inhibited microbial Fe(III) reduction and pyruvate metabolism. Interestingly, compared to cultures without mineral amendments or those with glass beads alone, the addition of quartz sand enabled strain Z6 to better withstand the environmental stress caused by elevated CO 2 , promoting pyruvate fermentation and iron reduction. In addition to an increased pH, the formation of siderite, hematite and vivianite was also observed in the bioactive systems. Although both glass beads and quartz sand were primarily composed of silica, differences in the mineral structure, elemental composition and acid neutralization capacity rendered quartz sand more chemically active and unexpectedly led to greater CO 2 sequestration.

CO2 stress↗

DOE OSTI2026

Marginal Soils Index Analysis & Geospatial Data

This data package contains output files associated with Mongird et al. (in prep) organized into four dataset directories. Each dataset is described in more detail below. 1. Marginal Soils Index Analysis Description: This folder contains a csv file with land needs and availability by state, power generating technology type, and scenario in 2050 when suitable siting areas are additionally constrained to areas with increasing levels of soil marginality. Files: msi_constrained_siting_availability_2050.csv Variables: Scenario - Projected 2050 scenario name State - US state abbreviation Technology - Generating technology type solar = solar photovoltaic gas_cc_re = natural gas combined cycle (recirculating cooling) wind = onshore wind gas_cc_ccs_re = natural gas combined cycle with carbon capture sequestration (recirculating cooling) gas_cc_dry = natural gas combined cycle with (dry cooling) gas_cc_pond = natural gas combined cycle with (pond cooling) coal_conv_ccs_re = conventional coal with carbon capture sequestration (recirculating cooling) Req_Capacity_MW - The amount of rated capacity required in 2050 of the given technology type in the given state and under the given scenario from the capacity expansion plan Req_Capacity_Factor - The assumed capacity factor (fraction between 0 and 1) for the given technology type in the given state and under the given scenario by the capacity expansion plan Req_Land_km2 - The amount of land required (in km-squared) to host the required generating capacity that is capable of meeting the specified capacity factor for the given technology type in the given state and under the given scenario Req_Energy_TWh - Product of Req_Capacity_MW, Req_Capacity_Factor, and 8760/1e6 for the given technology type in the given state and under the given scenario MSI_Case - The level of MSI that siting the given technology is additionally constrained to, where >0 means siting is additionally constrained to suitable land areas that have an MSI value greater than 0 >=1 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 1 >=2 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 2 >=3 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 3 Soil Attribute Rasters Description: This folder contains geospatial raster files for individual soil parameters upscaled to the listed grid resolution (30m or 1 km). 1 km resolution files are a spatial average of non-missing 30m resolution values. All raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. Files: avg_cond_raster_ .tif - Average conductivity of the saturation extract across all soil horizons within a depth of 40 inches, measured in mmhos/cm max_cond_raster_ .tif - Maximum conductivity of the saturation extract across all soil horizons within a depth of 40 inches, measured in mmhos/cm min_ph_raster_ .tif- Min pH values across all soil horizons within a depth of 40 inches. avg_ph_raster_ .tif - Average pH value across all soil horizons within a depth of 40 inches. max_ph_raster_ .tif- Max pH value across all soil horizons within a depth of 40 inches. erosion_factor_raster_ .tif - Product of k-factor and percent slope flood_freq_raster_ .tif - Number of months of the year during which the area is commonly, frequently, or very frequently flooded. max_sar_raster_ .tif - Maximum sodium adsorption ratio across all horizons within a depth of 40 inches rock_frac_raster_ .tif - Fraction of the upper 6 inches of soil composed of rock fragments larger than 3 inches. temp_regime_raster_ .tif - Soil temperature regime with the following key: 0 = pergelic 1 = gelic 2 = cryic 3 = frigid 4 = isofrigid 5 = mesic 6 = isomesic 7 = thermic 8 = isothermic 9 = hyperthermic 10 =isohyperthermic Marginal Soils Index Rasters Description: This folder contains geospatial raster files of the Marginal Soils Index at the listed grid resolution (30m or 1 km). 1 km resolution files are a spatial average of 30m resolution. Both raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. A value of 0 indicates that there were no soil attributes present that indicate marginal soil. NA values indicate that data was unavailable or bodies of water. Files: marginal_soils_index_30m_raster.tif marginal_soils_index_1km_raster.tif Marginal Soils Index Resource Potential Rasters Description: This folder contains geospatial raster files of the Marginal Soils Index + Resource Potential (MSI+RP) score at 1km resolution for geothermal, solar, and wind technologies. Raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. NA values indicate that the location is not suitable for siting the given technology due to policy, environmental, socioeconomic, topological, and other constraints regardless of soil marginality level. Areas with values greater than or equal to zero represent the product of the normalized MSI value and the normalized resource potential value. Files: geothermal_msi_ep_score_raster.tif solar_msi_ep_score_raster.tif wind_msi_ep_score_raster.tif Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Agriculture↗

DOE OSTI2026

Marginal Soils Index Analysis & Geospatial Data

This data package contains output files associated with the Mongird et al. paper entitled "Can US power grid expansion avoid prime agricultural lands?" and is organized into four dataset directories. Each dataset is described in more detail below. 1. Marginal Soils Index Analysis Description: This folder contains a csv file with land needs and availability by state, power generating technology type, and scenario in 2050 when suitable siting areas are additionally constrained to areas with increasing levels of soil marginality. Files: msi_constrained_siting_availability_2050.csv Variables: Scenario - Projected 2050 scenario name State - US state abbreviation Technology - Generating technology type solar = solar photovoltaic gas_cc_re = natural gas combined cycle (recirculating cooling) wind = onshore wind gas_cc_ccs_re = natural gas combined cycle with carbon capture sequestration (recirculating cooling) gas_cc_dry = natural gas combined cycle with (dry cooling) gas_cc_pond = natural gas combined cycle with (pond cooling) coal_conv_ccs_re = conventional coal with carbon capture sequestration (recirculating cooling) Req_Capacity_MW - The amount of rated capacity required in 2050 of the given technology type in the given state and under the given scenario from the capacity expansion plan Req_Capacity_Factor - The assumed capacity factor (fraction between 0 and 1) for the given technology type in the given state and under the given scenario by the capacity expansion plan Req_Land_km2 - The amount of land required (in km-squared) to host the required generating capacity that is capable of meeting the specified capacity factor for the given technology type in the given state and under the given scenario Req_Energy_TWh - Product of Req_Capacity_MW, Req_Capacity_Factor, and 8760/1e6 for the given technology type in the given state and under the given scenario MSI_Case - The level of MSI that siting the given technology is additionally constrained to, where >0 means siting is additionally constrained to suitable land areas that have an MSI value greater than 0 >=1 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 1 >=2 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 2 >=3 means siting is additionally constrained to suitable land areas that have an MSI value greater than or equal to 3 Soil Attribute Rasters Description: This folder contains geospatial raster files for individual soil parameters upscaled to the listed grid resolution (30m or 1 km). 1 km resolution files are a spatial average of non-missing 30m resolution values. All raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. Files: avg_cond_raster_ .tif - Average conductivity of the saturation extract across all soil horizons within a depth of 40 inches, measured in mmhos/cm max_cond_raster_ .tif - Maximum conductivity of the saturation extract across all soil horizons within a depth of 40 inches, measured in mmhos/cm min_ph_raster_ .tif- Min pH values across all soil horizons within a depth of 40 inches. avg_ph_raster_ .tif - Average pH value across all soil horizons within a depth of 40 inches. max_ph_raster_ .tif- Max pH value across all soil horizons within a depth of 40 inches. erosion_factor_raster_ .tif - Product of k-factor and percent slope flood_freq_raster_ .tif - Number of months of the year during which the area is commonly, frequently, or very frequently flooded. max_sar_raster_ .tif - Maximum sodium adsorption ratio across all horizons within a depth of 40 inches rock_frac_raster_ .tif - Fraction of the upper 6 inches of soil composed of rock fragments larger than 3 inches. temp_regime_raster_ .tif - Soil temperature regime with the following key: 0 = pergelic 1 = gelic 2 = cryic 3 = frigid 4 = isofrigid 5 = mesic 6 = isomesic 7 = thermic 8 = isothermic 9 = hyperthermic 10 =isohyperthermic Marginal Soils Index Rasters Description: This folder contains geospatial raster files of the Marginal Soils Index at the listed grid resolution (30m or 1 km). 1 km resolution files are a spatial average of 30m resolution. Both raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. A value of 0 indicates that there were no soil attributes present that indicate marginal soil. NA values indicate that data was unavailable or bodies of water. Files: marginal_soils_index_30m_raster.tif marginal_soils_index_1km_raster.tif Marginal Soils Index Resource Potential Rasters Description: This folder contains geospatial raster files of the Marginal Soils Index + Resource Potential (MSIxRP) score at 1km resolution for geothermal, solar, and wind technologies. Raster files use the USA Contiguous Albers Equal Area Conic (ESRI:102003) projection. NA values indicate that the location is not suitable for siting the given technology due to policy, environmental, socioeconomic, topological, and other constraints regardless of soil marginality level. Areas with values greater than or equal to zero represent the product of the normalized MSI value and the normalized resource potential value. Files: geothermal_msi_rp_score_raster.tif solar_msi_rp_score_raster.tif wind_msi_rp_score_raster.tif Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Agriculture↗

DOE OSTI2021

Combined Techno-Economic Analysis and Life Cycle Assessment of an Integrated Direct Air Capture System with Advanced Algal Biofuel Production

The continuous increase in carbon dioxide (CO2) concentration in the atmosphere since the First Industrial Revolution correlates convincingly with the ongoing rise of the Earth's global average temperature contributing to climate change. As a response, carbon capture and sequestration (CCS) technologies are being implemented to mitigate anthropogenic CO2 emissions by capturing CO2 from high emitting point sources, such as power plants, refineries, and cement factories, and subsequently buried in geological formations underground for long term storage. Recently direct air capture (DAC) has emerged as a promising alternative technology that captures CO2 directly from the atmosphere for use or sequestration. This study investigates the potential synergistic benefits of integrating a solid amine-based DAC system with advanced algal biofuel production in photobioreactors (PBRs). DAC utilization allows the removal of atmospheric CO2 while also decoupling algae production facilities from anthropogenic point CO2 sources and avoiding the cost and logistics challenges of transporting CO2 long distances to remote facilities. Techno-economic analysis and life cycle assessment are performed to assess the economic and environmental benefits of heat and mass integration between the DAC and the PBR for biofuel production. The DAC-PBR system integration also considers on-site flue gas handling options, DAC capital utilization, the tradeoff between centralization or decentralization of key unit operations, and the PBR array. This presentation will discuss how optimizing DAC-PBR process integration can enhance the algal biofuel's economic and environmental sustainability and the prospect of DAC enabling the circular carbon economy.

algal biofuel production↗

DOE OSTI2022

The Power of Being Negative: Producing H2 and Sequestered Carbon from Biomass and Waste Resources

Hydrogen can serve as an enabler for decarbonizing our economy, especially for the transportation, manufacturing, and chemicals sectors. Low-carbon H2 can be generated through a number of pathways, including (1) steam methane reforming coupled with carbon capture and sequestration, (2) water electrolysis utilizing renewable electricity, (3) methane pyrolysis, and (4) biomass and waste conversion. However, only the biomass and waste conversion pathway offers a route to carbon-negative H2 when coupled with carbon capture and sequestration (CCS), thus providing both a service to our environment through carbon dioxide removal from the atmosphere and a valuable downstream product in H2.

biomass↗

DOE OSTI2022

Monitoring Carbon Storage Sites with Time-Lapse Gravity Surveys

Monitoring technologies are deployed on carbon storage sites to demonstrate that the injected CO2 is contained within the storage complex during and long after the active life of the project and that the CO2 behavior is conforming to expectations. Detection of gravity variations accompanying mass redistributions in the subsurface caused by fluid movements in reservoirs provides a unique means for monitoring the dynamics of a carbon sequestration site. Time-lapse gravity surveys have a long history of effectively monitoring temporal density changes in the subsurface. However, monitoring of carbon storage sites with time-lapse gravity surveys using surface or borehole measurements remains limited, although it gained more attention during the last decade. With the recent progresses in instrument capabilities, and the increased number of fluid storage sites (e.g., wastewater, CO2, natural gas, etc.), the gravity method (with its unique ability to remotely estimate both the mass of an injected fluid and its location) has a promising future in monitoring. This chapter provides background information regarding the gravity method and discusses its modeling and application to various carbon storage sites.

Carbon Sequestration, time-lapse gravity, monitori↗

DOE OSTI2025

Upconversion of non-recycled MSW paper fractions into biochar via slow pyrolysis and life cycle analysis: Pathways to net negative GHG emission

This study presents an integrated and sustainable approach to valorizing non-recycled municipal solid waste (MSW), a heterogeneous and underutilized waste stream destined for landfilling, by converting it into valuable biochar resources. Specifically, we investigated the upcycling of nonrecycled paper waste based on compositional analysis into four major fractions: high cellulose, high lignin, high contamination, and high ash content papers. These fractions were then homogenized and subjected to slow pyrolysis. The high cellulose fraction (36.1 %) was the most abundant, and contained 66.7 % cellulose, while the high lignin fraction showed the highest lignin (12.1 %) and carbon content (44 %), resulting in highest energy value of 17.4 MJ kg −1 . Biochar yields ranged from 25.6 % to 35.6 %, with the high ash fraction producing the highest yield and alkalinity (pH ≈ 11.2) due to its higher mineral content. Elemental analysis revealed enhanced carbon content up to 76.9 % and reduced oxygen and hydrogen, confirming effective carbonization. The high lignin-derived biochar showed the highest aromatic carbon content (82.8 %) and greater structural stability, while contaminated and ash-rich fractions exhibited dense, low-porosity surfaces due to the presence of contaminants and minerals. Spectroscopic analysis revealed degradation of carbohydrates, disappearance of cellulose peaks and formation of aromatic and mineral derived phases. The scaled life cycle process yielded a global warming potential (GWP) of 119.3 kg CO 2 -eq per ton of dry paper waste, offset by soil carbon sequestration of − 556.41 kg CO 2 -eq, resulting in a net impact of − 427.36 kg CO 2 -eq. This represents a net carbon removal exceeding by ~186 % the emissions associated with landfilling paper waste with electricity generation.

09 BIOMASS FUELS↗

DOE OSTI2022

Pore-Scale Microenvironments Control Anthropogenic Carbon Mineralization Outcomes in Basalt

Thin sections and hand samples from 50 sidewall cores from the Wallula Basalt Pilot Demonstration, a basaltic carbon sequestration demonstration, provided the opportunity for the in-depth analysis of carbon mineralization induced by the injection of supercritical CO 2 . In this study, we used optical petrography and scanning electron microscopy to characterize the physical and chemical characteristics of the basalt components influenced by carbon mineralization reactions from all available hand samples and thin sections within the three CO 2 injection zones and caprock flow interiors. We found extensive carbonate mineralization, mostly in the form of nodules that were shown to be chemically zoned: Ca-dominant in the core regions and Ca-bearing Fe-dominant in the outer regions. Carbonate mineralization also took the form of fracture-filling carbonate cement, and acicular aragonite was also observed. Overall, we clarified the structural and paragenetic relationships between newly formed minerals, identifying a new fibro-palagonite-like, poorly crystalline silicate phase that grew on the carbonate nodules and pore-lining zeolites. Here we observed Fe-dominant carbonate precipitates surrounding acicular aragonite and rhombohedral Ca-carbonate cores, whereas previous studies of these zoned nodules did not observe these structures. A comprehensive accounting of the carbon mineralization products is vital to understand and predict the behavior of supercritical CO 2 in the subsurface given both the diversity of the host rock between and within injection zones, especially considering that the morphology and chemistry of the diverse precipitates are influenced by the pore-scale microenvironments of the basalt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Water availability modulates maximum canopy heights of low-elevation Amazonian second-growth forests

Tropical second-growth forests of the Amazon sequester large amounts of carbon and are important carbon sinks, contributing substantially to climate change mitigation, biodiversity conservation, and providing crucial ecosystem services. Deforestation due to selective logging and shifting cultivation is expanding second-growth forest areas in tropical forest regions, which if well managed, regenerate rapidly over time. Maximum forest canopy height is an important metric of biomass and carbon accumulation in second-growth forests and is strongly influenced by water availability. The water limitation hypothesis explains the positive influence of water availability on maximum tree heights and has been examined and demonstrated at a small-scale using field data, and at a global scale, with limited accuracy, using remote sensing data in tropical ecosystems. However, this hypothesis concerning maximum canopy height has not been much studied at regional and national scales for tropical second-growth forests. In this study, we leveraged NASA GEDI spaceborne lidar data across the Brazilian Amazon and derived second-growth forest relative height metrics for delineating the influence of water availability, second-growth forest age, and topographic elevation on maximum canopy height. Water availability was found to significantly influence the maximum canopy height of second-growth forest trees, of age range from 30 to 35 years, at elevations less than 500 m and maximum precipitation thresholds of 1500 mm. Our results indicate that changing precipitation patterns or increased drought conditions under different climate change regimes could impact forest structure, plant communities, ecosystem functioning, and carbon sequestration capabilities of tropical second-growth forests in the Amazon.

Carbon sequestration↗

DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗

DOE OSTI2020

Optimal carbon storage reservoir management through deep reinforcement learning

Model-based optimization plays a central role in energy system design and management. The complexity and high-dimensionality of many process-level models, especially those used for geosystem energy exploration and utilization, often lead to formidable computational costs when the dimension of decision space is also large. This work adopts elements of recently advanced deep learning techniques to solve a sequential decision-making problem in applied geosystem management. Specifically, a deep reinforcement learning framework was formed for optimal multiperiod planning, in which a deep Q-learning network (DQN) agent was trained to maximize rewards by learning from high-dimensional inputs and from exploitation of its past experiences. To expedite computation, deep multitask learning was used to approximate high-dimensional, multistate transition functions. Both DQN and deep multitask learning are pattern based. As a demonstration, the framework was applied to optimal carbon sequestration reservoir planning using two different types of management strategies: monitoring only and brine extraction. Both strategies are designed to mitigate potential risks due to pressure buildup. Results show that the DQN agent can identify the optimal policies to maximize the reward for given risk and cost constraints. Finally, experiments also show that knowledge the agent gained from interacting with one environment is largely preserved when deploying the same agent in other similar environments.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

A parametric study of accelerated carbonation in alkali-activated slag

Resistance to carbonation is one important attribute that low-CO{sub 2} cement alternatives must possess, and is particularly crucial for cement alternatives subjected to aggressive CO{sub 2} concentrations such as those used in construction of oil wells and wells for below ground carbon sequestration. Here, a parametric study of alkali-activated slag (AAS) carbonation in aggressive environments has been conducted to examine (i) calcium carbonate polymorphism using X-ray diffraction (XRD) and Fourier-transform infrared (FTIR) spectroscopy, and (ii) the extent of calcium carbonate formation and CO{sub 2} adsorption using thermogravimetric analysis (TGA). A range of AASs have been studied by varying the magnesium content of the slag, the activator type (sodium hydroxide and sodium silicate), the activator concentration, and the curing time prior to carbonation. It was uncovered that both (i) magnesium from the slag and (ii) silica from the activating solution are needed to reduce the propensity for the sodium-containing calcium-alumino-silicate-hydrate gel to undergo decalcification.

36 MATERIALS SCIENCE↗

DOE OSTI2024

The intrinsic mechanical properties of hydromagnesite, Mg 5 (CO 3 ) 4 (OH) 2 ·4H 2 O, a key phase of reactive MgO carbonate cement

To potentially enable CO 2 sequestration, reactive MgO carbonate cement is emerging as an alternative binder to Portland cement. Understanding the mechanical properties of its binding phase is critical for understanding the strength development and performing materials design for reactive MgO cement systems; however, the intrinsic mechanical properties of hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ·4H 2 O), a key binding phase, remain unexplored. Here the present study utilized synchrotron-based high-pressure X-ray diffraction to determine the unit cell-scale, intrinsic mechanical properties of hydromagnesite for the first time. Up to hydrostatic loading of 7.7 GPa, the bulk modulus of hydromagnesite was determined as 59 GPa or 71 GPa fitted using the second-order or third-order Birch-Murnaghan equation of state, which we contextualize with binding phases in various cement systems. The experiment results are applicable in materials design of low-carbon concrete and valuable for the validation and calibration of atomistic models.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

DOE OSTI2022

In situ inelastic neutron scattering of mixed CH 4 –CO 2 hydrates

An abundant source of CH 4 can be found in natural hydrate deposits. Recent demonstration of CH 4 recovery from hydrates via CO 2 exchange has revealed the potential as a fuel source that also provides a medium for carbon sequestration. It is vital to understand the structural and dynamic impacts of guest variation in CH 4 , CO 2 , and mixed hydrates and link the results to the stability of various deposits in nature, harvesting methane, and sequestering CO 2 . Molecular vibrations are examined in CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates at 5 and 190 K and Xe hydrates for comparison. Inelastic neutron scattering (INS) is an ideal spectroscopy technique to observe the dynamic modes in the hydrate structure and enclathrated CH 4 , as it is extremely sensitive to 1 H. The presence of CO 2 in hydrates tightens the lattice. It introduces more active librational modes to the host lattice, while hindering the motion of CH 4 in mixed CH 4 -CO 2 hydrate at 5 K. At 190 K, a large broadening of the CH 4 librational modes indicates disorder in the structure leading to dissociation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Techno-economic life cycle assessment of CO 2 -EOR operations towards net negative emissions at farnsworth field unit

Optimizations of CO 2 Water Alternating Gas(WAG)- systems with multi-objectives of incremental recovery and maximization of CO 2 storage are challenging. Here, the incorporation of a total Greenhouse gas (GHG) life cycle assessment is mostly ignored leading to inaccurate estimation of overall net carbon emissions of their operations. In this study, the effect of a total GHG life cycle assessment on a multi-objective CO 2 -WAG optimization with integrated techno-economic assessment (TEA) which factors carbon tax credit is conducted. A life cycle assessment (LCA) was conducted utilizing a 20 -year optimized post history matched data from a high fidelity reservoir simulation model. Using data generated from the optimum result, a techno-economic life cycle analysis was further conducted. The first scenario classified as the base model had an estimated 81% of purchased CO 2 sequestered. The results through a comprehensive techno-economic LCA model yielded a net estimate of 73% of purchased CO 2 . The optimized forecasted model which considered key operational and reservoir factors such as WAG ratio, injection rates and periods, and well specification resulted in an improved sequestration of 92% of purchased CO 2 . However, this also dropped to 84% after taking it through LCA. These results clearly indicate a significant amount of net CO 2 is not accounted for when operations are not analyzed through LCA. From the LCA, direct flaring volumes of CO 2 , energy consumption and efficiency of unit equipment were noticed to be the major causes of these reductions. Considering ten main sources of energy as source of energy generation, a comparative techno-eco LCA was conducted. The results confirmed a lower net volume and NPV for energy sources with higher carbon footprints and vice versa. Thus, a total LCA of CO 2 -WAG greatly influences net storage factor of purchased CO 2 and hence project NPV where tax credit/incentives per ton of CO 2 sequestered is considered. Although operational conditions are optimized for best results, there are significant factors that leads to minimization of net storage factor. This study therefore provides an insightful information for optimizing CO 2 -WAG multi-objectives to achieve minimum GHG emission.

02 PETROLEUM↗