Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
SEARCH · Engineering Papers
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.
Multidimensional, experimental and modeling evaluation of permeability evolution, the Caney Shale Field Lab, OK, USA
Carbon capture and storage (CCS) is a technological strategy to reduce CO 2 emissions from hydrocarbon-based energy production and non-energy industrial emitters, like fertilizer and cement production. To expand the geographic opportunities for CCS, offshore geologic settings are increasingly considered for saline CO 2 storage, particularly where legacy oil and gas infrastructure may be repurposed for CO 2 transport and injection. In this context, the Gulf of Mexico (GoM), United States (U.S.), may offer tremendous opportunity for CCS because the surrounding states are among the largest CO 2 producers in the U.S. and there has been extensive oil and gas development over the last half century. Nevertheless, there remains significant uncertainty in how GoM reservoir/seal systems will respond to industrial-scale CO 2 injections, particularly in the context of meso-scale (1–10s of meters) permeability variations that are difficult to identify prior to injection. Such permeability variations can have significant impacts on CO 2 migration and fluid pressure propagation, both of which govern the overall effectiveness of the CCS project. This research uses ensemble simulation methods to quantify the influence of spatially variable and a priori unknown permeability fields on CO 2 plume development and fluid pressure propagation when CO 2 is injected at 1 MMt per year for ten years. Results show a number of characteristic patterns that reveal the substantial influence of both near- and far-field permeability on reservoir performance. By focusing on basin characteristics that are common among offshore basins worldwide, results from this study suggest that offshore CCS may be a feasible carbon management strategy in regions where offshore basins are proximal to industrial CO 2 sources.
The Petra Nova Project (Project) is a commercial scale post-combustion carbon capture project developed by a joint venture between NRG Energy, Inc. (NRG) and JX Nippon Oil Exploration (EOR) Limited (JX). The Project is designed to separate and capture carbon dioxide (CO2) from an existing coal-fired unit’s flue gas slipstream at NRG’s W.A. Parish Electric Generating Station (WAP) located southwest of Houston, Texas. The captured CO2 is dried, compressed, and transported via an 81-mile pipeline to the West Ranch oilfield (West Ranch) in Jackson County, Texas, where it is injected to boost oil production. The Project, which is partially funded by a grant (Grant) from the United States (U.S.) Department of Energy (DOE) under the Clean Coal Power Initiative (CCPI) Round 3, uses the Kansai Mitsubishi Carbon Dioxide Recovery advanced amine-based CO 2 absorption technology (KM-CDR Process®), which was jointly developed by Mitsubishi Heavy Industries, Ltd. (MHI) and the Kansai Electric Power Co. Inc., to treat and capture at least ninety percent (90%) of the CO 2 from a 240-megawatt equivalent (MWe) flue gas slipstream off of Unit 8 at WAP. When operating at full capacity, the Project captures approximately 5,200 short tons of CO 2 per day, which would otherwise be emitted into the atmosphere, representing the largest commercial scale deployment of post-combustion CO2 capture technology at a coal power plant to date. Under the Grant, the Project was managed in 3 phases: (1) Phase 1: Project Definition / Front End Engineering Design (FEED) (2) Phase 2: Detailed Engineering, Procurement & Construction (3) Phase 3: Demonstration and Monitoring. On December 29, 2016, commercial operation of the Project was achieved, ending Phase 2 and starting Phase 3, a 3-year demonstration period running from January 1, 2017 through December 31, 2019. The key objectives of Phase 3 were to (a) demonstrate the specific advanced technologies constructed during Phase 2 and (b) monitor the injected CO 2 at West Ranch to demonstrate technologies and protocols for monitoring, verification, and accounting (MVA). As of the end of Phase 3, Petra Nova captured 3,904,978 short tons of CO 2 (3,542,537 metric tons) that was transported to West Ranch. To support the DOE obligation to monitor, verify, and account for the sequestered CO 2 at West Ranch, Petra Nova contracted with the Bureau of Economic Geology (in the Jackson School of Geosciences at The University of Texas at Austin) to (a) design a monitoring program, (b) draft an MVA Plan for DOE review and approval, and (c) working with Petra Nova and the operator of West Ranch to manage and report on the MVA activity. This report discussed the technical aspects of the project during each of the 3 phases of the project as identified above.
Growing interest in offshore geologic carbon sequestration (GCS) motivates risk assessment of large-scale subsea CO 2 well blowouts or pipeline ruptures. For major leaks of CO 2 from wells or pipelines, significant fluxes of CO 2 may occur from the sea surface depending on water depth. In the context of risk assessment of human health and safety, we have used previously simulated coupled well-reservoir and water column model results as a source term for dense gas dispersion of CO 2 above the sea surface. The models are linked together by one-way coupling, that is, output of one model is used as input to the next model. These first-of-their-kind coupled flow results are applicable to assessing the hazard of CO 2 to people at and downwind of the sea surface location of emission. Hazard is quantified by plotting the downwind dispersion length (DDL), which we define in the study as the distances from the emission source to the point at which the emitted CO 2 has been diluted to 5% and 1.5% in air by volume. Here results suggest that large-scale blowouts in shallow water (10 m) may cause hazardous CO 2 plumes extending on the order of several hundred meters downwind. Details of the modeling show DDL has a maximum for windspeed (at an elevation of 10 m) of approximately 5 m/s, with smaller DDL for both weaker and stronger winds. This is explained by the fact that wind favors transport but also causes dispersion; therefore there is a certain wind speed that maximizes DDL.
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.
Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.
Geological CO 2 sequestration (GCS) can help mitigate global warming and enhance methane recovery from coal beds. However, few studies have linked the effects of CO 2 to surface chemistry changes controlling wetting behavior in deep coal beds. Contact angles (CAs) of CO 2 /N 2 -high volatile bituminous coal-water systems were measured under different temperatures and pressures. The surface chemistry and physical structure of coals were characterized to investigate changes in physicochemical properties and their relations with wettability after reactions. For N 2 treatment, the time-dependence of static and dynamic CAs were insignificant, ranging within 4°. For gaseous CO 2 treatment, the static CAs and the average advancing angles increased slightly. With supercritical (sc) CO 2 , both the static and dynamic CAs increased significantly, and θ adv changed to intermediate-wet (92°). Reactions with minerals exposed to scCO 2 resulted in greater surface roughness and heterogeneity, greater contact angle hysteresis and more surface sites occupied by scCO 2 rather than H 2 O. Increases in hydrophobic functional groups and decreases in hydrophilicity were shown by FTIR spectra, reflecting the shedding of polar oxygen-containing functional groups, reduction of hydrogen bonds, and increasing percentage of hydrocarbons. XRD patterns obtained following scCO 2 -treatment showed that crystallite growth and molecular polymerization were higher toward graphite-like. The calculated structural parameters of functional groups and crystallites both showed elevated coal rank. Changes in crystallite structure, notably higher carbon content and decreased negative surface charge, are unfavorable for water-wetting. Finally, this study contributes to understanding surface chemistry changes responsible for decreased wettability during CO 2 -enhanced coal bed methane recovery and GCS in coal reservoirs.
The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.
Geological storage of carbon dioxide (CO 2 ) in depleted gas reservoirs and deep saline aquifers is a key part of global decarbonization efforts. As carbon capture and storage advances toward commercial-scale deployment, the credibility and scalability of laboratory experiments are increasingly vital for guiding safe and effective field implementation. This review offers a comprehensive, cross-scale evaluation of experimental methodologies, including core flooding, high-pressure, high-temperature systems, microfluidic visualization, and emerging systems such as multilayer commingled/compartmentalized core flooding, 3D-printed micromodels, and AI-powered digital twins. These innovations are demonstrated to enhance representativeness, reproducibility, and real-time insight, thereby addressing the limitations of conventional workflows. A critical analysis of methodological gaps, such as inconsistent pressure–temperature conditions, oversimplified brine chemistry, and a lack of standardization, reveals experimental sources of scale translation errors and performance uncertainty. By comparing the unique challenges of depleted gas reservoirs (such as low water saturation and legacy well leakage) to those of saline aquifers (including pressure buildup and caprock integrity), this review identifies formation-specific priorities for experimental design. Novel contributions include a synthesis of best practices, integration strategies for model calibration, and recommendations for standardizing core handling, saturation procedures, and reporting protocols. Furthermore, this work serves as a guide for developing robust, field-relevant experimental strategies that can increase the deployment and regulatory acceptance of CO 2 storage technologies at scale.
Geological CO 2 sequestration in deep saline aquifers has been extendedly investigated to reach the goal of carbon neutral as large amount of CO 2 can be reduced in a short time. Here, this study investigates the feasibility of water-alternating-gas (WAG) injection and brine extraction on enhancing CO 2 storage in a target deep saline aquifer, Minnelusa Sandstone in the Powder River Basin (PRB) of Wyoming. An integrated numerical model is developed to account for four scenarios: (1) CO 2 continuous injection (CI) through one injection well, (2) WAG injection through one injection well, (3) CO 2 CI through one injection well and brine extraction through one producer, (4) WAG injection through one injection well and brine extraction through one producer. Simulation results suggest that WAG injection and brine extraction corporately make it possible to enhance CO 2 injectivity while securing CO 2 storage safety. For instance, WAG injection considerably reduces structural trapping contribution while enhances dissolved and residual trapping contributions. As for brine extraction, it can decrease the maximum averaged reservoir pressure by 37% and 22% for CI and WAG injection, respectively. Besides, sensitive analyses of the operational parameters for the fourth scenarios are performed. Results reveal that, when the amount of total CO 2 injection has been predetermined, CO 2 injection time and rate within each period of WAG should be small whereas it is preferable to have large water injection time and rate. This is to secure desirable dissolved and residual trapping contributions (store CO 2 with safety). Sensitive degree results confirm that CO 2 injection rate, water injection time, and CO 2 injection time of WAG injection exhibit significant effects on CO 2 storage whereas the impacts of producer bottom-hole pressure (BHP), water injection rate, and CO 2 -water injection time are relatively weaker. This study not only sheds light on achieving double-win goal: enhance CO 2 injectivity and store CO 2 with safety, but also can be a critical reference for other CO 2 storage in deep saline aquifers.
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.
Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.
Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.
This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.
Geological features play a pivotal role in determining the feasibility of deploying CO₂ direct air capture (DAC) technologies, primarily because they influence the availability of cost-effective energy sources, such as natural gas and geothermal energy, and also due to the potential for CO₂ sequestration. Many regions face challenges due to variable weather conditions including seasonal temperature fluctuations, high or low humidity, and sub-ambient temperatures. These extremes can reduce DAC performance or even lead to catastrophic events. Aqueous solvents considered for DAC systems are particularly vulnerable to seasonal variations in colder climates, where the solvent may underperform or freeze. It is therefore essential to investigate the CO₂ capture efficiency of aqueous solvents across a broad range of environmental temperatures, spanning sub-zero to hot conditions (>30 °C). In this study, DAC operation is examined using a high-flux solvent–air crossflow contactor under two major weather scenarios: (i) cold conditions below 0 °C and (ii) hot conditions above 30 °C. A parametric study is conducted to investigate the contactor performance regarding CO₂ removal efficiency, uptake capacity, and reaction kinetics versus temperature when the air velocity through the contactor exceeds 1 m/s. The efficacy of the contactor is systematically investigated using various anti-freeze amino-acid solvent formulations. A mass-transfer mechanistic model is developed to assess the process performance over a wide temperature range and propose scalable design guidelines. Machine learning is also employed to identify key parameters affecting the CO₂ capture efficiency. It is shown that air velocity and temperature are the primary factors influencing CO₂ uptake. Based on performance data obtained under subfreezing temperatures, a technoeconomic analysis is conducted to evaluate the feasibility of using aqueous solvents in seasonal cold regions. In conclusion, the findings of this study provide valuable insights into siting considerations for deploying solvent-based DAC, thereby contributing to the advancement of sustainable carbon removal solutions.
Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.