Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geomechanical Modeling”

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 55 records · Page 3

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗

Unlocking the Tight Oil Reservoirs of the Powder River Basin, Wyoming

The project focused on detailed geologic characterization, geomechanical studies, well completion optimization, stimulation monitoring, and field development strategies. A key aspect involved partnerships with industry and academic collaborators such as Occidental Petroleum, Southern Illinois University, Britt Rock Mechanics and Piri Technologies. Data acquisition included drilling, logging, coring, deployment of fiber optics, and microseismic monitoring. The project emphasized feedback loops for continuous model updating and integration of economic evaluations to guide development strategies.

02 PETROLEUM↗

A decade of progress in understanding and managing legacy well integrity for geologic carbon storage

This study reviews a decade of research progress in legacy well integrity and risk management for geologic carbon storage (GCS) to commemorate the 20 th anniversary of the Intergovernmental Panel on Climate Change’s 2005 Special Report on Carbon Capture and Storage. In the past ten years, legacy well research has benefited from global efforts to constrain emissions from abandoned oil and gas wells, a continued focus on well materials performance in the presence of CO 2 -rich fluids, and practical experience gained through GCS implementation. Field measurements of abandoned well emissions show that leakage is not universal or catastrophic but forms a continuum of low-to-moderate fluxes that depend on isolation integrity and environmental attenuation. Materials research has constrained the conditions under which Portland cements exhibit self-sealing and non-sealing behaviors, and has identified the impact of geomechanical properties, non-uniform pathway apertures, multi-phase flow, and impurities in the CO 2 stream, on leakage pathways as important new areas for investigation. GCS projects at brownfield sites have inspired the creation of new workflows that integrate various tools and technologies to manage legacy well leakage risks. GCS implementation has also motivated a push towards scenario-based well modeling that directly informs permit applications. These advances inspire new research questions for the coming decade, particularly around the level of legacy well leakage risk that is environmentally acceptable and tolerable to stakeholders when sequestering millions of tonnes of CO 2 annually.

Carbon capture and storage↗

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor

The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.

Geomechanics↗

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

Hydromechanical Modeling of Fault Rupture in Geologic CO 2 Sequestration: A Comparison of Two Failure Criteria

Abstract The large‐scale implementation of geological carbon sequestration has raised concerns about potential fault activation and induced seismicity, which could compromise storage integrity and pose seismic risks. We theoretically compared two failure criteria, Mohr‐Coulomb (MC) and Modified Cam‐Clay criteria (MCC), to assess fault rupture during CO 2 storage. Both criteria characterize fault behavior in a specific stress regime but differ in reducing the complexity of fault rupture to a few key mechanisms. Using a coupled hydromechanical model, we demonstrate that the choice of a failure criterion and the physics of fault weakening associated with these criteria strongly condition fault response to a given fluid injection. MC mainly relates rupture to friction, while MCC relates rupture to fault poro‐plasticity. Our findings highlight that the selection of a failure criterion, being inherently subjective, can significantly alter the predicted fault behavior during CO 2 storage, thereby impacting the reliability of geomechanical risk assessments.

Cao, M. [Energy Geosciences Division Lawrence Berk↗

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗

Thermo-hydro-mechanical analysis of subsurface ice-based thermal energy storage

Ice-based thermal energy storage systems are widely utilized for cooling and managing peak electrical demand globally, offering daily or weekly storage capabilities for both individual homes and larger office buildings. However, scaling these systems for district-level cooling or integrating them with renewable energy sources presents challenges, especially in accommodating larger volumes and addressing seasonal storage requirements in densely populated urban areas. This paper proposes a novel solution by evaluating subsurface ice-based thermal energy storage, in which the underground is subjected to seasonal freeze/thaw cycles. However, these cycles may influence ground behavior, affecting pore pressure and inducing ground movement. To systematically investigate these challenges, we enhance the TOUGH-FLAC simulator by integrating water/ice phase change capabilities and updating the effective stress–strain constitutive relation. Both modifications are validated against analytical solutions or experimental data. Through numerical simulations spanning a decade with ten seasonal freeze/thaw cycles, we evaluate the performance and long-term stability of a generic subsurface ice-based thermal energy storage system, considering factors such as ground permeability, freezing pipe spacing, freeze/thaw damage, and glycol solution temperature. The simulations indicate that ice formation induces pore pressure variations that drive seasonal surface heave and settlement, controlled by ground permeability, pipe spacing, and glycol solution temperature, along with tensile and localized shear deformation around freeze pipes. This highlights the need for accurate ground property characterization and geomechanical analysis for subsurface ice-based thermal energy storage.

Ground freezing↗

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗

Jackson, L., Johnson, M.B., Latrach, A., Grimes, D., Martinez, C., and Mclaughlin, J.F., 2024, Multidisciplinary geotechnical data collection, curation, and analysis for conformity with the regulatory framework for geologic carbon storage in Wyoming, USA: Geological Society of America Abstracts with Programs. Vol. 56, No. 5, 2024, doi: 10.1130/abs/2024AM-405024

Title: Multidisciplinary Geotechnical Data Collection, Curation, and Analysis for Conformity with the Regulatory Framework for Geologic Carbon Storage in Wyoming, USA. Text: Construction and operation of wells for geologic sequestration of carbon dioxide necessitate that they are permitted under the Environmental Protection Agency’s Underground Injection Control Class VI requirements. Class VI wells conform to stringent requirements to ensure long-term safety and integrity of the storage site and the protection of Underground Sources of Drinking Water. Entities pursuing Class VI permitting must provide comprehensive geologic site characterization, including regional geologic structure and stratigraphy, aquifer information, reservoir and confining unit geomechanical properties, geochemical analyses, assessment of trapping capacity and mechanisms, and a variety of other of multidisciplinary geotechnical data. The Wyoming Class VI Site Characterization Database Project is focused on developing a geologic site characterization database of geotechnical information, which has been compiled and verified from established, public databases/entities and scientific literature to expedite Class VI permitting in Sweetwater County within the Greater Green River Basin of southern Wyoming. The preliminary suite of compiled data from 14,000 wells includes 8,000 wells with logs and 7,250 wells with formation tops, ~70 wells with core data (e.g., X-Ray diffraction, petrographic, and petrophysical data), ~2,500 water analyses, ~740 seismic events data, and ~520 bottom-hole temperature measurements. Future work on—and stemming from—this project will include new core analyses, calculation and interpolation of subsurface temperature gradients, mechanical earth models, geochemical simulations, storage capacity estimation, stratigraphic column generation and correlation, and construction of subsurface maps. Finally, this work will help to inspire and facilitate subsurface data compilation and curation beyond Sweetwater County, Wyoming.

42 ENGINEERING↗

Introduction to Vector Field Visualization

Vector field visualization techniques are essential to help us understand the complex dynamics of flow fields. These can be found in a wide range of applications such as study of flows around an aircraft, the blood flow in our heart chambers, ocean circulation models, and severe weather predictions. The vector fields from these various applications can be visually depicted using a number of techniques such as particle traces and advecting textures. In this tutorial, we present several fundamental algorithms in flow visualization including particle integration, particle tracking in time-dependent flows, and seeding strategies. For flows near surfaces, a wide variety of synthetic texture-based algorithms have been developed to depict near-body flow features. The most common approach is based on the Line Integral Convolution (LIC) algorithm. There also exist extensions of LIC to support more flexible texture generations for 3D flow data. This tutorial reviews these algorithms. Tensor fields are found in several real-world applications and also require the aid of visualization to help users understand their data sets. Examples where one can find tensor fields include mechanics to see how material respond to external forces, civil engineering and geomechanics of roads and bridges, and the study of neural pathway via diffusion tensor imaging. This tutorial will provide an overview of the different tensor field visualization techniques, discuss basic tensor decompositions, and go into detail on glyph based methods, deformation based methods, and streamline based methods. Practical examples will be used when presenting the methods; and applications from some case studies will be used as part of the motivation.

Kao, David↗

Carbon Storage in Fold‐and‐Thrust Belts: An Overlooked Gigatonne Storage Opportunity

This study presents numerical investigations of the trapping characteristics of fold-and-thrust belt structures, defining three carbon capture and storage (CCS) play types that could be used to store commercial volumes (millions of tonnes) of CO 2 . Specifically, we present simulations of CO 2 storage in three fold-and-thrust belt models comprising a thrust-ramp, duplex, and thrust-fold geometry. To constrain these play types in realistic geology, each model is based on a study site, including a novel investigation of a greenfield saline reservoir in Virginia, USA, being considered for commercial carbon storage and two well-characterized petroleum fields: the Wilburton field in Oklahoma, USA, and the Incahuasi field in Bolivia. Our results provide insight into several key parameters, such as the long-term security of injected CO 2 in these geologies and injection strategies for maximizing storage efficiency while reducing pressure-related risk. These results improve the understanding of CCS in fold-and-thrust belt storage sites globally by describing general storage parameters that may be applied to site-specific projects. We find that thrust-ramp geometries may securely trap CO 2 through solubility and hydrodynamic trapping under suitable reservoir conditions, duplex structures may store some quantities of CO 2 but are pressure-constrained, and that thrust-ramp structures may store large quantities of CO 2 by maximizing fetch volume, which simultaneously lowers geomechanical risk by reducing pressure buildup along zones of weakness.

42 ENGINEERING↗