Search

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. We use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

Fracture injection/shut-in tests are commonly used to measure the state of stress in the subsurface. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of (Hayashi and Haimson, 1991). In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. The calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

Fundamental mechanical and creep behavior of minerals at elevated temperatures is relevant to many subsurface exploration and resource recovery processes, yet remains largely elusive. This study uses instrumented indentation to measure microscale deformation and nanoscale creep of biotite, albite, and quartz, which are the major mineral phases of granitic rocks, at temperatures up to 400 °C. The results show that the elastic modulus of biotite is not only the lowest at room temperature but also decays the most at elevated temperature (i.e., by 25% reduction at 400 °C) compared to that of albite and quartz. The creep deformation within 20 s doubles when the temperature increases from 20 °C to 400 °C for albite, quartz, and biotite, and biotite showed at least three times the overall creep deformation compared to quartz and albite under the same temperature. The indentation creep deformation can be well characterized using a logarithmic time model, showing both transient and constant-rate creep deformation. The secondary creep rate for biotite is about three times that of albite and quartz at the same temperature, and as the temperature increases from 20 °C to 400 °C, the rate of secondary creep becomes approximately 10 times faster for the three tested minerals. Both transient creep deformation and the rate of secondary creep can be empirically correlated to the elastic moduli, which allows quick estimates of the temperature-dependent creep behavior of granitic minerals using their elastic properties, especially when those creep constants are not readily available at elevated temperatures. These results enhance the understanding of temperature-dependent creep deformation at small scales and provide insight into mineral-level damage in granitic rocks due to temperature and long-term deformation.

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150-250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: This increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: A thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. These experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150-250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: This increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: A thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. These experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, we show the classic fingering and transitions in flow regimes at (Formula presented.) (capillary number) from −7.2 to −3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO2 (scCO2) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

Estimation of CO2 storage capacity in fractured porous reservoirs requires a better understanding of the CO2-water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO2-water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. The core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO2-water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO2 saturations in the matrix increaze with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO2 storage efficiency in fractured porous media.

We investigate the influence of foliation orientation and fine-scale folding on the static and dynamic elastic properties and unconfined strength of the Poorman schist. Measurements from triaxial and uniaxial laboratory experiments reveal a significant amount of variability in the static and dynamic Young’s modulus depending on the sample orientation relative to the foliation plane. Dynamic P-wave modulus and S-wave modulus are stiffer in the direction parallel to the foliation plane as expected for transversely isotropic mediums with average Thomsen parameters values 0.133 and 0.119 for epsilon and gamma, respectively. Static Young’s modulus varies significantly between 21 and 117 GPa, and a peculiar trend is observed where some foliated sample groups show an anomalous decrease in the static Young’s modulus when the symmetry axis (x3-axis) is oriented obliquely to the direction of loading. Utilizing stress and strain relationships for transversely isotropic medium, we derive the analytical expression for Young’s modulus as a function of the elastic moduli E1, E3, ν31, and G13 and sample orientation to fit the static Young’s modulus measurements. Regression of the equation to the Young’s modulus data reveals that the decrease in static Young’s modulus at oblique symmetry axis orientations is directly influenced by a low shear modulus, G13, which we attribute to shear sliding along foliation planes during static deformation that occurs as soon as the foliation is subject to shear stress. We argue that such difference between dynamic and static anisotropy is a characteristic of near-zero porosity anisotropic rocks. The uniaxial compressive strength also shows significant variability ranging from 21.9 to 194.6 MPa across the five sample locations and is the lowest when the symmetry axis is oriented 45° or 60° from the direction of loading, also a result of shear sliding along foliation planes during static deformation.

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO2 (scCO2) is used in enhanced oil recovery. We performed scCO2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from −8.1 to −4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO2 saturation increase at both capillary fingering and crossover flow regimes (−8.1 ≤ logCa ≤ − 6.1), the one heterogeneous wetting to scCO2 results in ~0.09 saturation increase only at the crossover flow regime (−7.1 ≤ logCa ≤ − 6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > − 6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO2 saturation and distribution suggest the importance of wettability on CO2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO2 (scCO2) is used in enhanced oil recovery. We performed scCO2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from −8.1 to −4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO2 saturation increase at both capillary fingering and crossover flow regimes (−8.1 ≤ logCa ≤ − 6.1), the one heterogeneous wetting to scCO2 results in ~0.09 saturation increase only at the crossover flow regime (−7.1 ≤ logCa ≤ − 6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > − 6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO2 saturation and distribution suggest the importance of wettability on CO2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiagvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiagvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

Dissolution trapping is one of the most important mechanisms for geological carbon storage (GCS). Recent laboratory and field experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and coupled scCO2 dissolution and water flow, i.e., scCO2 dissolution at local pores/pore throats creating new water-flow paths, which in turn enhance dissolution by increased advection and interfacial area. However, the impacts of pore-scale characteristics on these coupled processes have not been investigated. In this study, imbibition and dissolution experiments were conducted under 40 °C and 9 MPa using a homogeneous/isotropic hexagonal micromodel, two homogeneous elliptical micromodels with low or high anisotropy, and a heterogeneous sandstone-analog micromodel. The four micromodels, initially saturated with deionized (DI)-water, were drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at different rates with logCa (the capillary number) ranging from −6.56 to −4.34. Results show that bypass of scCO2 by displacing water is the dominant mechanism contributing to the residual CO2 trapping, triggered by heterogeneity in pore characteristics or pore-scale scCO2-water distribution. Bypass can be enhanced by pore heterogeneity or reduced by increasing transverse permeability, resulting in relatively low (

Dissolution trapping is one of the most important mechanisms for geological carbon storage (GCS). Recent laboratory and field experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and coupled scCO2 dissolution and water flow, i.e., scCO2 dissolution at local pores/pore throats creating new water-flow paths, which in turn enhance dissolution by increased advection and interfacial area. However, the impacts of pore-scale characteristics on these coupled processes have not been investigated. In this study, imbibition and dissolution experiments were conducted under 40 °C and 9 MPa using a homogeneous/isotropic hexagonal micromodel, two homogeneous elliptical micromodels with low or high anisotropy, and a heterogeneous sandstone-analog micromodel. The four micromodels, initially saturated with deionized (DI)-water, were drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at different rates with logCa (the capillary number) ranging from −6.56 to −4.34. Results show that bypass of scCO2 by displacing water is the dominant mechanism contributing to the residual CO2 trapping, triggered by heterogeneity in pore characteristics or pore-scale scCO2-water distribution. Bypass can be enhanced by pore heterogeneity or reduced by increasing transverse permeability, resulting in relatively low (

Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and a prolonged depletion of residual scCO2. In this study, pore-scale scCO2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO2 into the sandstone-analogue pore network initially saturated by water without dissolved CO2 (dsCO2). During the experiments, time-lapse images of dye intensity, reflecting water pH, were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO2 dissolution and phase equilibrium occurs when scCO2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO2 dissolution at phase interfaces and diffusion of dsCO2 at the pore scale (10–100 µm) observed after scCO2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. This finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO2 at narrow pore throats.

Recently, both core- and pore-scale imbibition experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and a prolonged depletion of residual scCO2. In this study, pore-scale scCO2 dissolution and mass transfer under drainage conditions were investigated using a two-dimensional heterogeneous micromodel and a novel fluorescent water dye with a sensitive pH range between 3.7 and 6.5. Drainage experiments were conducted at 9 MPa and 40 °C by injecting scCO2 into the sandstone-analogue pore network initially saturated by water without dissolved CO2 (dsCO2). During the experiments, time-lapse images of dye intensity, reflecting water pH, were obtained. These images show non-uniform pH in individual pores and pore clusters, with average pH levels gradually decreasing with time. Further analysis on selected pores and pore clusters shows that (1) rate-limited mass transfer prevails with slowly decreasing pH over time when the scCO2-water interface area is low with respect to the volume of water-filled pores and pore clusters, (2) fast scCO2 dissolution and phase equilibrium occurs when scCO2 bubbles invade into water-filled pores, significantly enhancing the area-to-volume ratio, and (3) a transition from rate-limited to diffusion-limited mass transfer occurs in a single pore when a medium area-to-volume ratio is prevalent. The analysis also shows that two fundamental processes – scCO2 dissolution at phase interfaces and diffusion of dsCO2 at the pore scale (10–100 µm) observed after scCO2 bubble invasion into water-filled pores without pore throat constraints – are relatively fast. The overall slow dissolution of scCO2 in the millimeter-scale micromodel can be attributed to the small area-to-volume ratios that represent pore-throat configurations and characteristics of phase interfaces. This finding is applicable for the behavior of dissolution at pore, core, and field scales when water-filled pores and pore clusters of varying size are surrounded by scCO2 at narrow pore throats.

In modeling of geological carbon storage, dissolution of supercritical CO2 (scCO2) is often assumed to be instantaneous with equilibrium phase partitioning. In contrast, recent core-scale imbibition experiments have shown a prolonged depletion of residual scCO2 by dissolution, implying a non-equilibrium mechanism. In this study, eight pore-scale scCO2 dissolution experiments in a 2D heterogeneous, sandstone-analog micromodel were conducted at supercritical conditions (9MPa and 40°C). The micromodel was first saturated with deionized (DI) water and drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at constant flow rates after scCO2 drainage was completed. High resolution time-lapse images of scCO2 and water distributions were obtained during imbibition and dissolution, aided by a scCO2-soluble fluorescent dye introduced with scCO2 during drainage. These images were used to estimate scCO2 saturations and scCO2 depletion rates. Experimental results show that (1) a time-independent, varying number of water-flow channels are created during imbibition and later dominant dissolution by the random nature of water flow at the micromodel inlet, and (2) a time-dependent number of water-flow channels are created by coupled imbibition and dissolution following completion of dominant imbibition. The number of water-flow paths, constant or transient in nature, greatly affects the overall depletion rate of scCO2 by dissolution. The average mass fraction of dissolved CO2 (dsCO2) in water effluent varies from 0.38% to 2.72% of CO2 solubility, indicating non-equilibrium scCO2 dissolution in the millimeter-scale pore network. In general, the transient depletion rate decreases as trapped, discontinuous scCO2 bubbles and clusters within water-flow paths dissolve, then remains low with dissolution of large bypassed scCO2 clusters at their interfaces with longitudinal water flow, and finally increases with coupled transverse water flow and enhanced dissolution of large scCO2 clusters. The three stages of scCO2 depletion, common to experiments with time-independent water-flow paths, are revealed by zoom-in image analysis of individual scCO2 bubbles and clusters. The measured relative permeability of water, affected by scCO2 dissolution and bi-modal permeability, shows a non-monotonic dependence on saturation. The results for experiments with different injection rates imply that the non-equilibrium nature of scCO2 dissolution becomes less important when water flow is relatively low and the time scale for dissolution is large, and more pronounced when heterogeneity is strong.

In modeling of geological carbon storage, dissolution of supercritical CO2 (scCO2) is often assumed to be instantaneous with equilibrium phase partitioning. In contrast, recent core-scale imbibition experiments have shown a prolonged depletion of residual scCO2 by dissolution, implying a non-equilibrium mechanism. In this study, eight pore-scale scCO2 dissolution experiments in a 2D heterogeneous, sandstone-analog micromodel were conducted at supercritical conditions (9MPa and 40°C). The micromodel was first saturated with deionized (DI) water and drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at constant flow rates after scCO2 drainage was completed. High resolution time-lapse images of scCO2 and water distributions were obtained during imbibition and dissolution, aided by a scCO2-soluble fluorescent dye introduced with scCO2 during drainage. These images were used to estimate scCO2 saturations and scCO2 depletion rates. Experimental results show that (1) a time-independent, varying number of water-flow channels are created during imbibition and later dominant dissolution by the random nature of water flow at the micromodel inlet, and (2) a time-dependent number of water-flow channels are created by coupled imbibition and dissolution following completion of dominant imbibition. The number of water-flow paths, constant or transient in nature, greatly affects the overall depletion rate of scCO2 by dissolution. The average mass fraction of dissolved CO2 (dsCO2) in water effluent varies from 0.38% to 2.72% of CO2 solubility, indicating non-equilibrium scCO2 dissolution in the millimeter-scale pore network. In general, the transient depletion rate decreases as trapped, discontinuous scCO2 bubbles and clusters within water-flow paths dissolve, then remains low with dissolution of large bypassed scCO2 clusters at their interfaces with longitudinal water flow, and finally increases with coupled transverse water flow and enhanced dissolution of large scCO2 clusters. The three stages of scCO2 depletion, common to experiments with time-independent water-flow paths, are revealed by zoom-in image analysis of individual scCO2 bubbles and clusters. The measured relative permeability of water, affected by scCO2 dissolution and bi-modal permeability, shows a non-monotonic dependence on saturation. The results for experiments with different injection rates imply that the non-equilibrium nature of scCO2 dissolution becomes less important when water flow is relatively low and the time scale for dissolution is large, and more pronounced when heterogeneity is strong.

Injection of carbon dioxide (CO2) into saline aquifers confined by low- permeability cap rock will result in a layer of CO2 overlying the brine. Dissolution of CO2 into the brine increases the brine density, resulting in an unstable situation in which more-dense brine overlies less-dense brine. This gravitational instability could give rise to density-driven convection of the fluid, which is a favorable process of practical interest for CO2 storage security because it accelerates the transfer of buoyant CO2 into the aqueous phase, where it is no longer subject to an upward buoyant drive. Laboratory flow visualization tests in transparent Hele-Shaw cells have been performed to elucidate the processes and rates of this CO2 solute-driven convection (CSC). Upon introduction of CO2 into the system, a layer of CO2-laden brine forms at the CO2-water interface. Subsequently, small convective fingers form, which coalesce, broaden, and penetrate into the test cell. Images and time-series data of finger lengths and wavelengths are presented. Observed CO2 uptake of the convection system indicates that the CO2 dissolution rate is approximately constant for each test and is far greater than expected for a diffusion-only scenario. Numerical simulations of our system show good agreement with the experiments for onset time of convection and advancement of convective fingers. There are differences as well, the most prominent being the absence of cell-scale convection in the numerical simulations. This cell-scale convection observed in the experiments may be an artifact of a small temperature gradient induced by the cell illumination.