Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir potential”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Big Data Meets Geothermal Exploration (CRADA Final Report)

As part of the Cyclotron Road program, Zanskar Geothermal & Minerals, Inc. investigated the application of micro-earthquake and ambient noise seismology methods to imaging and characterizing the structural characteristics and hydrothermal flux of subsurface faults. Significant advances in what could be resolved were enabled by two major developments in seismology: 1) the availability of large-n arrays of low-cost seismometers, and 2) the availability of increased computational power and semi-automated data reduction algorithms. In tandem, these advances may improve the signal-to-noise ratio and spatial precision of the data collected and enable higher-resolution characterization of subsurface fracture systems and their spatio-temporal evolution. These tools supported efforts to reduce dry-hole risk and to improve wellfield productivity for geothermal resource development. In particular, two applications of these advances were evaluated: 1) fracture-seismic imaging, which was used to detect ambient emissions from fluid-filled fractures, and 2) reservoir tomography, which used information about travel paths, source locations, and source parameters of micro-earthquakes to identify areas of enhanced permeability. Integration of these methods provided guidance for siting wells and served as prior constraints for reservoir models, informing forecasts of power potential and production and injection strategies aimed at minimizing temperature decline and improving overall resource productivity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Structure and Stability of Iron Fluoride at High Pressure–Temperature and Implication for a New Reservoir of Fluorine in the Deep Earth

Fluorine (F) is the most abundant halogen in the bulk silicate Earth. F plays an important role in geochemical and biological systems, but its abundance and distribution in the terrestrial mantle are still unclear. Recent studies suggested that F reservoirs in the deep mantle are potentially hosted in terrestrial oxide minerals, especially in aluminous bridgmanite. However, the knowledge about the formation and stability field of fluoride in the Earth’s interior is rare. In this study, we combine in situ laser-heated diamond anvil cell, synchrotron X-ray diffraction, and first-principles structure search to show that a new tetragonal structure of FeF 3 is stable at pressures of 78–130 GPa and temperatures up to ~1900 K. Simulation predicted the tetragonal phase takes a much denser structure due to the rotation of FeF 6 octahedral units. The equations of states of tetragonal FeF3 are determined by experiment and verified by simulation. Our results indicate that FeF 3 can be a potential key phase for storing F in the Earth’s lower mantle and may explain some mantle-derived magma with high F concentration.

58 GEOSCIENCES↗
NASA NTRS1977

Space industrialization

Space industrialization is defined as the use of space flight and the space environment for commercial or utilitarian purposes in contrast to other uses such as gains in basic scientific knowledge, national defense, or exploration. Some unique attributes of space that make it amenable to industrial use include overview of the earth, the 'zero gravity' effect, potential for near perfect vacuum, unlimited reservoir for disposal of waste products, availability of essentially uninterrupted flow of solar energy, and the 'perpetual motion' characteristic of orbital mechanics. The role of human participation in assembling and maintaining the large sophisticated systems that will be required for future space industrialization needs is considered.

Disher, J. H.↗
DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗
DOE OSTI2024

Iron-sandstone synergy: Advancing in-situ hydrogen production from natural gas via electromagnetic heating

Here, to address the escalating demands for decarbonization in the petroleum industry, a carbon-zero technology, known as in-situ hydrogen (H 2 ) production via electromagnetic (EM)-assisted catalytic heating, has recently been proposed for generating and extracting clean H 2 directly from petroleum reservoirs. Although preliminary techno-economic analyses show significant potential of this emerging technology for clean and affordable hydrogen, the fundamentals of natural gas conversion to H 2 in the presence of reservoir rocks are poorly understood. In this study, we explore the synergy between sandstone and artificial iron-based catalysts in enhancing in-situ H 2 production from methane (CH 4 ) cracking under EM irradiation. The dynamic behaviors of sandstone under EM heating are comprehensively investigated, including its thermal behaviors, thermal runaway (TR) phenomenon, gas generation during TR, and energy consumption. We found that sandstone demonstrates an evident natural catalytic effect for promoting CH 4 conversion to H 2 , enabling H 2 production starting at about 394 °C. The natural catalytic role of iron minerals in sandstone is elucidated using various advanced characterization techniques. Remarkably, when adding iron catalysts into the sandstone, the highest H 2 concentration and CH 4 conversion reaches 91 mol.% and 80%, respectively, at a temperature of 666 °C, while they are 50 mol.% and 35%, respectively, for the sample consisting of iron catalysts and quartz at the same level of temperature. This result indicates a strong iron-sandstone synergy and a potential to stimulate H 2 production by leveraging this synergy. Throughout the experimental process, the generation of carbon oxides (CO and CO 2 ) is negligible. These findings pave a pathway towards future pilot for carbon-zero in-situ H 2 production from sandstone gas reservoirs.

03 NATURAL GAS↗
DOE OSTI2025

Emerging applications: Neuromorphic computing and reservoir computing

The emergence of doped hafnium oxide (HfO 2 )-based ferroelectric films has enabled highly scalable and silicon-compatible ferroelectric devices, opening new frontiers in neuromorphic and reservoir computing. Among these, ferroelectric field-effect transistors (FeFETs) are particularly promising due to their analog memory characteristics and unique polarization dynamics. These properties make FeFETs ideal candidates for artificial synapses in neuromorphic architectures, supporting deep neural networks and spiking neural networks based on leaky-integrate-and-fire (LIF) mechanisms. Beyond neuromorphic computing, FeFETs also play a crucial role in physical reservoir computing, leveraging their intrinsic nonlinear and history-dependent behavior for efficient real-time learning. This approach offers significant advantages for time-series processing and edge artificial intelligence (AI) applications, addressing the growing need for energy-efficient computing. As a result, this article explores the principles, key demonstrations, and future potential of FeFET-based neuromorphic and reservoir computing, highlighting their impact on next-generation AI hardware.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Investigating CO2 Reservoirs at Gale Crater and Evidence for a Dense Early Atmosphere

One of the most compelling features of the Gale landing site is its age. Based on crater counts, the formation of Gale crater is dated to be near the beginning of the Hesperian near the pivotal Hesperian/Noachian transition. This is a time period on Mars that is linked to increased fluvial activity through valley network formation and also marks a transition from higher erosion rates/clay mineral formation to lower erosion rates with mineralogies dominated by sulfate minerals. Results from the Curiosity mission have shown extensive evidence for fluvial activity within the crater suggesting that sediments on the floor of the crater and even sediments making up Mt. Sharp itself were the result of longstanding activity of liquid water. Warm/wet conditions on early Mars are likely due to a thicker atmosphere and increased abundance of greenhouse gases including the main component of the atmosphere, CO2. Carbon dioxide is minor component of the Earth's atmosphere yet plays a major role in surface water chemistry, weathering, and formation of secondary minerals. An ancient martian atmosphere was likely dominated by CO2 and any waters in equilibrium with this atmosphere would have different chemical characteristics. Studies have noted that high partial pressures of CO2 would result in increased carbonic acid formation and lowering of the pH so that carbonate minerals are not stable. However, if there were a dense CO2 atmosphere present at the Hesperian/Noachian transition, it would have to be stored in a carbon reservoir on the surface or lost to space. The Mt. Sharp sediments are potentially one of the best places on Mars to investigate these CO2 reservoirs as they are proposed to have formed in the early Hesperian, from an alkaline lake, and record the transition to an aeolian dominated regime near the top of the sequence. The total amount of CO2 in the Gale crater soils and sediments is significant but lower than expected if a thick atmosphere was present at the Hesperian/Noachian boundary. Likewise, the absence of carbonates suggests that CO2- weathering processes similar to those present on Earth were not dominant. Instead it is possible that more exotic CO2 deposition has occurred driven by atmospheric photochemistry and/or degradation of organic carbon.

Niles, P. B.↗
DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

The large thermal gradients and high subsurface temperatures of the western region of the U.S. hold great potential for the implementation of enhanced geothermal systems (EGS). The development of these potential EGS resources requires stimulation of the reservoir to enhance permeability and it has been widely reported that a substantial amount of water will be required should conventional hydraulic stimulation be used. This presents a huge challenge and a high risk to the geothermal development because the water stress1 in these areas is already high or extremely high. The use of foam, a gas/liquid mixture predominantly composed of gas, in fracturing is considered and explored in this project as a potential approach to address water concerns with hydraulic stimulation in the development of EGS. This project, led by Oak Ridge National Laboratory (ORNL) in collaboration with Temple University, was awarded in an open lab call in 2018, and was part of the DOE GTO waterless stimulation initiative. The goal of the project was to demonstrate the feasibility of foam fracturing for EGS development through two primary tasks: Task 1: Laboratory study of the effectiveness of foam fracturing for representative geological materials, including cyclic pressurization using foam (led by ORNL) and Task 2: High temperature foam material selection and characterization (led by Temple University). In FY19, ORNL finished the critical review on serval issues associated with foam fracturing and the implementation of the proposed tasks in a lab study (Wang, et al., 2019), and completed the foam fracturing testing using cement as a model material (Wang, et al, 2020a). The work at ORNL was geared up to develop a brand-new foam testing system in FY20. The purchase of main components for the new system was finished in the first half of the FY20. The assembly of the foam testing system and foam fracturing testing were completed in the second half of the FY20 (Wang, et al., 2021a). Task 1 required the development of a test system which can be used to perform hydraulic fracturing of geological specimens with both water and foamed liquids at pressure up to 6,000 psi (41.4 MPa). The system possesses several capabilities that conventional injection systems lack for hydraulic fracturing. In addition to its ability to generate foam with controlled quality, it is capable of cycling pressure levels between specified values up to frequencies of 50 Hz. The latter capability was developed to evaluate the hypothesis that cyclic loading of samples would produce enhanced fracturing. The system consists of two sections: one for foam generation and another for foam injection. The foam is generated through separate control and pressurization of liquid and gas phases with controlled flow rates. The injection section is equipped with a low-flow Coriolis flow that monitors the density of foam to ensure the injection is in the range of target foam quality2. Experimental results of foam fracturing are reported for cylindrical granite specimens using water and aqueous N2 foam as the fracturing fluids. All experiments were performed for unconfined conditions. The effects of injection mode (i.e., monotonic vs cyclic pressurization) on breakdown pressure and failure response sample were investigated using water alone as a fracturing fluid and foams with a range of compositions. It was found that in the case of monotonic injection, the breakdown pressure of granite specimens tended to be slightly higher when fracturing with foam. Additionally, with a foam quality of 90%, the water use can be reduced by 50 to 84%, depending on hole size. On the other hand, it was observed that the breakdown pressure can be brought down to 70% of the monotonic breakdown pressure by using low cycle fatigue. Finally, discussions are presented regarding injectivity and water use reduction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Using Machine Learning to Identify Optimum Prospects for Offshore Geologic CO 2 Storage and Enhanced Oil Recovery in the Gulf of Mexico

The SECARB Offshore Partnership is a government-industry partnership focused on assembling the knowledge base required for secure, long-term, large-scale CO 2 subsea storage in saline formations and mature oil reservoirs. The SECARB Offshore Partnership is evaluating potential storage opportunities in the Central and Eastern planning areas of the outer continental shelf (OCS) of the Gulf of Mexico. The evaluation focuses on active and depleted oil and gas fields and potentially associated CO 2 -enhanced oil recovery (CO 2 -EOR), as well as, deep saline storage resources. As part of Subtask 3.0 (Offshore Storage Resources Characterization) and Subtask 4.0 (Risk Assessment, Simulation, and Modeling), Oklahoma State University (OSU) is developing a machine learning system that employs the SAS Institute’s Viya platform with visual data mining and machine learning to identify and characterize optimum prospects for enhanced oil recovery and stacked storage in saline reservoirs. This report summarizes this effort to date and discusses the many variables that can be used to characterize these storage objectives, as well as, the design vision of the Viya machine learning system.

02 PETROLEUM↗
DOE OSTI2024

Geothermal Reservoir Characterization at Wattenberg Area in Denver Basin

The geothermal potential of the Denver-Julesburg (D-J) Basin is largely overlooked. This study focuses on characterizing the geothermal potential in the greater Wattenberg area within the Denver Basin, using available well data and subsurface information and gas wells drilled to above 2500 m deep. The review reveals that the main structures comprise ENE-trending wrench faults, along with NNE, N, and NNW trending normal faults. A thermal anomaly is present bonded by Lafayette and Longmont wrench faults and overlying vitrinite reflectance anomaly. The geometry of thermal anomaly suggests that it is influenced by Longmont wrench fault and Lafayette fault, as well as the basin margin fault. Bottom hole temperatures are of 160 Degrees Celsius, with thermal gradients of and 62 Degrees Celsius/km. Reservoir temperature estimates from thermal gradients suggest a range of >200 Degrees Celsius at depths of 3 km.

Denver basin↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
DOE OSTI2022

Modeling Subsurface Performance of a Geothermal Reservoir Using Machine Learning

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells—increasing or decreasing the fluid flow rates across the wells—and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. In this study, we describe a new approach combining reservoir modeling and machine learning to produce models that enable such a strategy. Our computational approach allows us, first, to translate sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy, and second, to find optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an “open-source” reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 h, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 s. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs.

15 GEOTHERMAL ENERGY↗
NASA NTRS2013

Simplified Mantle Architecture and Distribution of Radiogenic Power

The mantle components that represent the source region of ocean island basalts (OIB) and feed hotspot volcanism are predicted to contain 160 +/- 20 (2σ(sub m)) ng/g Th, a heat‐producing element. This critical model composition indicates that the OIB source region (OSR) comprises a significant amount of recycled oceanic crust and constitutes 19(sup +3)(sub −2) (2σ(sub m))% of the mantle by mass. The mass fraction of this reservoir supports a mantle architecture with a basal thermochemical layering at an average depth of 2000 +/- 100 (2σ(sub m)) km or two thermochemical piles that extend up to midmantle levels. The hotspot source described here generates 10 pW/kg of radiogenic heat and supplies 7.3 TW to the planet's total surface heat flux. Given that the silicate portion of the Earth produces some 20.4 TW of radiogenic power, with 7.2 TW derived from the continental crust, the mantle source responsible for mid‐ocean ridge volcanism provides only 5.9 TW of radiogenic power (or <2 pW/kg). As a result, the source of hotspots generates >5x more radiogenic heat than the source of mid‐ocean ridges, thus contributing to the energetics that drive mantle convection and potentially the formation of long‐lived plumes via bottom heating of the modern mantle. The potential for a sequestered or unsampled mantle reservoir would impact the relative mass fractions of the source regions of OIB and mid‐ocean ridge volcanism but not the compositional model of the OSR presented here.

Arevalo, Ricardo, Jr.↗
DOE OSTI2021

SECARB Regional Project Assessment

The SECARB Cranfield Early Test Project was designed to safely demonstrate large-scale, long-term CO 2 injection and storage in a CO 2 -enhanced oil recovery (EOR) operation and associated saline reservoir. Storage during EOR operations has the potential to offset commercial costs, holding a significant promise for future large-scale deployment in the southeastern United States. Denbury Onshore LLC began CO 2 -EOR operations in July 2008 at the Cranfield oilfield located east of Natchez, Mississippi (Figure 1). The location offered an ideal setting to investigate storage monitoring at a commercial CO 2 -EOR operation. The field team, led by the Gulf Coast Carbon Center of the Bureau of Economic Geology at the University of Texas at Austin, deployed the largest suite of monitoring, verification, accounting (MVA) technologies to collect data to assess long-term CO 2 utilization and storage. The CO 2 injection interval is within the 3,000-meter-deep fluvial lower Tuscaloosa Formation, a large and regionally extensive saline reservoir, at a broad four-way structural closure that defines the Cranfield oilfield. Tests were conducted in the oil-producing area as well as in the non-oil-producing brine aquifer. The CO 2 used during EOR operation was transported via pipeline from the Jackson Dome, a natural source of CO 2 . The study focused on four integrated research areas: (1) the High-Volume Injection Test area (HiVIT); (2) the Detailed Area of Study (DAS); (3) the Geomechanical area; and (4) the near-surface observatory, also called the “P-site.” CO 2 injection occurred at the HiVIT and the DAS (Figure 2 and Figure 3). The DAS was located down-dip of the oil-water interface of the Cranfield oilfield. The well layout included one injection well and two down-dip observation wells. The team was the first of the regional partnerships to monitor more than 1 million metric tons of CO 2 . At the end of the project, the team monitored more than 11 million metric tons of CO 2 injected and 5 million metric tons of CO 2 stored. The Cranfield Project was the fifth project worldwide, and the first such project in the United States, to reach this CO 2 injection volume while being monitored to demonstrate storage effectiveness.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Genomic and Phenotypic Characterization of Shiga Toxin-Producing Escherichia albertii Strains Isolated from Wild Birds in a Major Agricultural Region in California

Escherichia albertii is an emerging foodborne pathogen. To better understand the pathogenesis and health risk of this pathogen, comparative genomics and phenotypic characterization were applied to assess the pathogenicity potential of E. albertii strains isolated from wild birds in a major agricultural region in California. Shiga toxin genes stx2f were present in all avian strains. Pangenome analyses of 20 complete genomes revealed a total of 11,249 genes, of which nearly 80% were accessory genes. Both core gene-based phylogenetic and accessory gene-based relatedness analyses consistently grouped the three stx2f-positive clinical strains with the five avian strains carrying ST7971. Among the three Stx2f-converting prophage integration sites identified, ssrA was the most common one. Besides the locus of enterocyte effacement and type three secretion system, the high pathogenicity island, OI-122, and type six secretion systems were identified. Substantial strain variation in virulence gene repertoire, Shiga toxin production, and cytotoxicity were revealed. Six avian strains exhibited significantly higher cytotoxicity than that of stx2f-positive E. coli, and three of them exhibited a comparable level of cytotoxicity with that of enterohemorrhagic E. coli outbreak strains, suggesting that wild birds could serve as a reservoir of E. albertii strains with great potential to cause severe diseases in humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗
DOE OSTI2021

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗