Engineering Papers⌕ Search

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.

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

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessment of CO 2 trapping mechanisms in partially depleted oil-bearing sands

The goal of this work is to evaluate various CO 2 sequestration mechanisms occurring in the Morrow B Sandstone in the Farnsworth Unit. A history-matched numerical simulation model was created using extensive geological, petrophysical, and operational data collected from the field. The numerical model is competent to investigate the impact of residual, structural-stratigraphic, solubility, and mineral trapping mechanisms on the fluid transportation dynamics and petrophysical property variations. The model forecasts the field response of 20 years of WAG injections. Afterward, all wells were shut-in, and the reservoir was allowed to evolve for 1000 years to investigate the fate of injected CO 2 . In this paper, we assess the impacts of various trapping mechanisms on oil recovery and CO 2 storage efficacy. By analyzing the results reported from the numerical simulation model, the in-situ fluid composition and mineralogy changes are also investigated. More importantly, we seek to confirm the petrophysical property variations due to the CO 2 injection with observations from laboratory measurements. The experiences gained from this study provide valuable insights regarding physiochemical storage induced by the CO 2 injection activities and serve as a benchmark case for future CO 2 enhanced oil recovery (EOR) projects involving reactive solute transport.

42 ENGINEERING↗
DOE OSTI2020

Water content of carbon dioxide at hydrate forming conditions

There is an interest to ensure sub-saturated water content in lines containing carbon dioxide in applications such as enhanced oil recovery and carbon sequestration, to reduce risks of hydrate blockage and corrosion. The water content of carbon dioxide at various temperatures and pressures has been measured in the past, but there is no consistent set of measurements that could be used for carbon dioxide storage and transportation design work. The solubility of water in a carbon dioxide rich gas phase at hydrate forming conditions was measured in this work. Pressures ranged from 12.06 to 29.30 bar along two isotherms, 1 °C and -7 °C, all within the gaseous carbon dioxide and hydrate stability zone. For the first time in these types of measurements, the solid phase was also characterized and confirmed to be carbon dioxide hydrate via X-ray computed tomography, simultaneous with water content measurements of the gas phase. Once carbon dioxide hydrate conversion had reached a maximum value (65% estimated by X-ray computed tomography), the equilibrium water content was measured. Prior to reaching this maximum carbon dioxide hydrate conversion, the water content in carbon dioxide was observed to decrease as liquid water converted to carbon dioxide hydrate. Furthermore, this slow conversion to hydrate, metastability of the hydrate phase, or unexpected phases may be responsible for the large discrepancy between prior data sets for similar carbon dioxide water content measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

DICE-gas turbine compound reheat combined cycle

Coal-fired Direct Injection Carbon Engine – Gas Turbine (DICE-GT) Compound-Reheat Combined Cycle (CRCC), is a combined cycle power plant comprising a multiplicity of coal-fired reciprocating internal combustion engines (RICE), which is commonly referred to by the acronym DICE (Direct Injection Carbon Engine); a natural gas-fired gas turbine (turbine and combustor); a heat recovery steam generator (HRSG); a steam turbine generator; and an integrally geared and intercooled centrifugal air compressor. Finally, the concept uses coal-water slurry/fluid in the DICE and natural gas (or hydrogen) in the gas turbine for unmatched efficiency (well above 50% net LHV) and modularity (block sizes of 120, 240 and 360 MW) with easy adaptability to carbon capture and sequestration with minimal additional cost, complexity and performance hit.

42 ENGINEERING↗
DOE OSTI2020

Optimization of energy requirements for CO 2 post-combustion capture process through advanced thermal integration

The energy optimization modeling work described here was performed to determine efficiency improvements that could be achieved for existing coal-fired power plants to retrofit a partial CO 2 capture from the post-combustion flue gas for carbon sequestration through thermal integration. The work presented includes optimization of the mono-ethanol amine (MEA)-based post-combustion CO 2 capture to reduce energy requirements that could be achieved at existing power plants by thermal integration of the steam turbine cycle, boiler, CO 2 compression train and post-combustion CO 2 capture process to offset efficiency and capacity losses that would be incurred by retrofit or implementation of post-combustion CO 2 capture. Furthermore, partial CO 2 capture, involving treatment of less than 100% of the flue gas leaving the plant and modular design of the CO 2 scrubbing system, was also investigated. Thermal integration of the steam turbine cycle with boiler and CO 2 compression train improved cycle and plant performance and offset, in part, the negative effects of post-combustion CO 2 capture. The best-analyzed integration options improved gross power output by 5% and net unit efficiency by 1.57%, relative to the conventional MEA process. Operating with 40% CO 2 capture increased gross power output by 11.6–14% (depending on the MEA thermal integration option), relative to the conventional MEA integration and 90% CO 2 capture. The improvement in net unit performance is larger compared to the improvement in turbine cycle performance because of the CO 2 compression work, which is also reduced by partial CO 2 capture.

42 ENGINEERING↗
DOE OSTI2021

Mechanisms and kinetics of citrate-promoted dissolution of a uranyl phosphate mineral

The formation of uranyl phosphate precipitates is a remediation strategy to limit the mobility of uranium in contaminated soils. However, exposure to organic ligands, like the plant exudate citrate, can remobilize the uranium. The purpose of this study is to provide a more thorough comprehension of citrate-promoted dissolution of a uranyl phosphate mineral, chernikovite [(H 3 O)(UO 2 )(PO 4 )·3H 2 O], by determining the extent of uranium release from chernikovite at a wide range of citrate concentrations. Here we have quantified the kinetics of dissolution and proposed potential mechanisms of chernikovite dissolution to gain a better understanding of the fate of uranyl phosphate precipitates in the environment. Batch dissolution and continuously stirred tank reactor (CSTR) experiments indicate that increasing citrate concentrations from 0.1 mM to 50 mM increases the concentration of dissolved uranium. However, at citrate concentrations of 10 mM and greater, the effectiveness of the ligand to enhance the dissolution decreases. Density functional theory modeling and Raman spectroscopy indicate a strong interaction between citrate and the uranyl ions at the surface of the uranyl phosphate. The interaction between citrate and uranyl suggests the formation of an alteration layer of uranyl-citrate surface complexes, which could impede uranium release at citrate concentrations of 10 mM or higher. Together, these results indicate that the citrate-promoted dissolution of chernikovite is a fast reaction that is hindered by a combination of surface saturation and a secondary-phase precipitation reaction at higher concentrations of citrate and reaction times greater than or equal to 72 min based on the hydraulic residence time of CSTR experiments. These findings improve our understanding of the stability of uranyl phosphate for use in sequestration of uranium groundwater plumes at contaminated sites.

58 GEOSCIENCES↗
DOE OSTI2022

Cadmium isotope systematics in sedimentary carbonate: Extending the utility of the cadmium isotope palaeo-productivity proxy

The cadmium (Cd) isotope compositions of seven Holocene (<11.7 ka) carbonate sediments, including the bulk carbonate, as well as the > 140 μm (zooplankton), <63 μm and < 20 μm (phytoplankton) fractions, and individual species of foraminifera isolated from the bulk carbonate were obtained for sub-samples processed using the Cd-cleaning and Mg/Ca cleaning methods. The respective Cd isotope systematics of each carbonate archive were then compared with the Cd isotope signature within the ocean using modern seawater data. We show that a constant Cd isotope fractionation factor of -0.25 ± 0.03 ‰ (2SE, n = 4) describes the partitioning of Cd from seawater into the carbonate phase ($Δ^{114}Cd_{CdC O_3 - seawater}$), allowing for accurate reconstructions of the past seawater nutrient regime. The Cd isotope signature of the < 63 µm and < 20 µm (phytoplankton) fractions can display comparable values to those recorded in the modern surface ocean, and can be used to model the magnitude of Cd uptake and utilisation in the past ocean. Individual species of planktic and benthic foraminifera, which inhabit different depths of the water column, can provide a full surface to deep water column reconstruction of the Cd nutrient regime, while infaunal benthic foraminifera can be used to interrogate the Cd isotope systematics of sediment porewaters. Bulk sediment and the > 140 µm (zooplankton) fraction also retain important information relating to the overlying water column Cd isotope systematics and can serve as essential archives for palaeo-productivity reconstructions where the separation of individual components from the bulk sediment is not possible. This new Cd isotope palaeo-productivity toolkit can considerably aid reconstruction efforts and vastly improves the ability to quantify the role of phytoplankton in modulating carbon sequestration through transitions in Earth’s climate evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Morphological and crystallographic controls in the replacement of calcite and aragonite by cerussite and otavite

Mineral replacement reactions are essential solid-fluid interactions in natural and industrial processes including metasomatism, diagenesis, and metamorphism, and sequestration of metal toxins from polluted water. In this paper, we explore the morphological evolution of aragonite and calcite (polymorphs of CaCO 3 ) during their replacement by cerussite (PbCO 3 , an isomorph of aragonite) and otavite (CdCO 3 , an isomorph of calcite) in acidic aqueous solutions (with initial pH 3 to 4). Observations using scanning electron and synchrotron-based transmission X-ray microscopies reveal the formation of pseudomorphic shells (~5 to 10 µm thick) of cerussite and otavite when Pb 2+ is reacted with either calcite or aragonite and when Cd 2+ is reacted with aragonite. The formation of pore space at the substrate-precipitate interface and within the precipitate is found to be key in promoting chemical exchange between the dissolving phase and the solution. In contrast, minimal reactivity is observed when Cd 2+ is reacted with calcite. These results demonstrate that the lack of substrate-precipitate epitaxy (either because of inconsistent symmetries or due to large lattice constant mismatch) is a critical factor that enables replacement of primary CaCO 3 minerals by Pb- or Cd-bearing secondary minerals while the presence of epitaxial relationships inhibit those reactions. When calcite and aragonite coexist, the replacement of aragonite by cerussite and otavite occurs preferentially over calcite. We postulate that the relative reactivity of calcite and aragonite to Pb 2+ and Cd 2+ can be determined by interplay between morphological and crystallographic properties of CaCO 3 materials as well as solution chemistry. Our results demonstrate that coupled dissolution-precipitation processes of carbonate minerals play an important role in controlling mobility and bioavailability of lead and cadmium in environments where acidic metal-contaminated water interacts with carbonate-rich rocks and soils.

58 GEOSCIENCES↗
DOE OSTI2025

The role of Ca-bridged organic matter in an alkaline soil, as revealed by multimodal chemical imaging

Mineral–organic matter (OM) studies have predominantly focused on acidic soils that are abundant in iron (Fe) oxides and aluminum (Al) oxides. We have probed mineral–OM interactions in an alkaline or calcareous soil of the Aridisols class. Unlike the role of Fe and Al, the role of Ca-minerals (particularly calcite), which are ubiquitous in alkaline soils, in OM sequestration is not well understood. Multiple recent model studies with aqueous Ca2+ or synthetic calcite and a suite of OM compounds have shown Ca-OM assemblages to be spatially correlated with calcite at the microscale. To study the chemical state of both Ca and Fe and their competing role in soil organic matter (SOM) stabilization, we performed laboratory characterization using x-ray diffraction, Mössbauer spectroscopy, x-ray photoelectron spectroscopy, scanning electron microscopy, and scanning transmission electron microscopy, alongside synchrotron-based microscale chemical imaging using scanning transmission x-ray microscopy combined with near-edge x-ray absorption fine structure. Ca mineral–organic associations were found to be ubiquitous in this system and are likely critical for understanding SOM stabilization/degradation in alkaline soils. From our findings on mineralogy, speciation, and the nature of Ca-OM bridging, we identified differences in C and Ca chemistry based on the relative location of OM to Ca minerals. The OM near the calcite crystal was enriched in lipid and protein moieties, Ca-OM next to Fe minerals displayed a strong contribution from aromatic compounds, while on the surface of microbes, the carbonate was believed to be of microbial in origin, as also suggested by preliminary works reporting on the formation of amorphous calcite or nano-calcite. In Ca-OM admixed with carbonate, it was difficult to distinguish Ca-associated OM from amorphous calcite or nano-calcite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗
DOE OSTI2022

Evaluation of accuracy and convergence of numerical coupling approaches for poroelasticity benchmark problems

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. Furthermore, this work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

58 GEOSCIENCES↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2020

Mechanistic modeling of CO 2 well leakage in a generic abandoned well through a bridge plug cement-casing gap

Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.

58 GEOSCIENCES↗
DOE OSTI2020

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗
DOE OSTI2021

A screening model for predicting injection well pressure buildup and plume extent in CO 2 geologic storage projects

We present a new screening model for predicting injection well pressure buildup and CO 2 plume migration for CO 2 geologic sequestration projects. The model requires only limited information and is quite accurate when compared to detailed simulation results. Such models can assist project developers during the early days of project planning (e.g., for 45Q related projects), and also help regulators perform simple checks against detailed numerical models. The screening model consists of two correlations: one for injectivity index in terms of the slope of the CO 2 fractional flow curve, and the second for total storage efficiency within the footprint of plume as a function of gravity number and slope of the CO 2 fractional flow curve. Using these two correlations, as well as a knowledge of some basic reservoir characteristics and estimates of fluid properties from standard correlations, the injection-well pressure buildup and CO 2 plume extent in the formation can be readily estimated. Finally, the new correlations show a good match with the results of the underlying simulations, and also provide good agreement with independent calculations for two example problems.

42 ENGINEERING↗
DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗
DOE OSTI2024

An experimental investigation on the CO 2 storage capacity of the composite confining system

Assuring secure containment of stored CO 2 is of paramount importance—for climate change mitigation, for permitting, and for reassuring the public. Regional seals, such as those sealing petroleum accumulations, have proven to be effective for securing CO 2 . However, the goal of permanent sequestration can also be satisfied with “composite confining systems” consisting of multiple, possibly discontinuous flow barriers that, in aggregate, create a system with very high permeability anisotropy and effectively retard the vertical migration of CO 2 . This study focuses on investigating the barrier characteristics necessary for effective containment of CO 2 , using both physical flow experiments and modified invasion percolation simulations. The simulations are calibrated to the physical experiments and are used to further extend the analysis. Results show that for a composite confining system, a) even barriers with low capillary entry pressure contrast to the underlying flow unit can divert rising CO 2 ; b) curved or anticlinal barrier topography can enhance CO 2 trapping; c) fining-upward gradations make little difference to barrier effectiveness or CO 2 retention; d) longer barriers retain more CO 2 regardless of barrier topography. Finally, field-scale simulations have demonstrated the importance of barrier length for increasing CO 2 storage capacity. Furthermore, the results presented can be used to inform the development of new screening criteria for characterization and effectiveness of composite confining systems.

58 GEOSCIENCES↗
DOE OSTI2025

Biogenic carbon capture at pulp mills via sodium spiking and oxy-fuel calcination

Over 13 million metric tons of biogenic carbon dioxide (CO 2 ) are mineralized yearly in United States (US) pulp mill recovery boilers as molten sodium carbonate. The mineralized CO 2 is released downstream in the rotary lime kiln where it can be captured at a relatively low cost for permanent sequestration. The use of biocarbon and bioenergy at pulp mills enables the possibility for atmospheric carbon removal when biogenic CO2 is captured and sequestered. We demonstrate the feasibility of capturing biogenic CO 2 via sodium spiking coupled with oxy-fuel calcination in the rotary lime kiln at existing kraft pulp mills. Sodium (Na) spiking elevates the Na ion loading in the kraft chemical looping process by replacing chlorine-based bleaching with a highly alkaline bleaching sequence. Chemical pulping processes are modeled and simulated to understand the technical limits of implementing sodium spiking in existing pulp mills. Each 1 % of sodium added to the bleaching operations increases the rate of CO 2 mineralization by 4 %, and the maximum increase in sodium content in the kraft process is 9 %. Without sodium spiking, estimated CO 2 capture costs are $131 and $107/mt CO 2 for air and oxy-fuel combustion, respectively. With the implementation of sodium spiking, the cost of CO 2 captured decreases by 27 % and 31 %, showing costs of $96 and $74/mt CO 2 for air and oxy-fuel combustion, respectively. Oxy-fuel combustion reduces the costs of CO 2 capture compared to air combustion, but merged deployment with sodium spiking seems to be the most cost-effective pathway to integrate carbon capture.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗