Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Economic feasibility of geologic carbon storage demands sustaining large storage rates without damaging caprock seals. Reactivation of pre-existing or newly formed fractures may provide a leakage pathway across caprock layers. In this study, we apply an equivalent continuum approach within a finite element framework to model the fluid-pressure-induced reactivation of pre-existing fractures within the caprock, during high-rate injection of super-critical CO2 into a brine-saturated reservoir in a hypothetical system, using realistic geomechanical and fluid properties. We investigate the impact of reservoir to caprock layer thickness, wellbore orientation, and injection rate on overall performance of the system with respect to caprock failure and leakage. We find that vertical wells result in locally higher reservoir pressures relative to horizontal injection wells for the same injection rate, with high pressure inducing caprock leakage along reactivated opening-mode fractures in the caprock. After prolonged injection, leakage along reactivated fractures in the caprock is always higher for vertical than horizontal injection wells. Furthermore, we find that low ratios of reservoir to caprock thickness favor high excess pressure and thus fracture reactivation in the caprock. Injection into thick reservoir units thus lowers the risk associated with CO2 leakage.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

An enriched finite element method is described for capillary hydrodynamics including dynamic wetting. The method is enriched via the Conformal Decomposition Finite Element Method (CDFEM). Two formulations are described, one with first-order accuracy and one with second-order accuracy in time. Both formulations utilize a semi-implicit form for the surface tension that is shown to effectively circumvent the explicit capillary time step limit. Sharp interface boundary conditions are developed for capturing the dynamic contact angle as the fluid interface moves along the wall. By virtue of the CDFEM, the contact line is free to move without risk of mesh tangling, but is sharply captured. Multiple problems are used to demonstrate the effectiveness of the methods.

Abstract not provided.

Abstract not provided.

Coupled reservoir and geomechanical simulations are significantly important to understand the long-term behavior of geologic carbon storage (GCS) systems. In this study, we performed coupled fluid flow and geomechanical modeling of CO₂ storage using available field data to (1) validate our existing numerical model and (2) perform parameter estimation via inverse modeling to identify the impact of key geomechanical (Young's modulus and Biot's coefficient) and hydrogeological (permeability and anisotropy ratio) properties on surface uplift and the pore pressure buildup at In Salah in Algeria. Furthermore, two sets of surface uplift data featuring low and high uplifts above two injection wells and the maximum change in the pore pressure due to CO₂ injection were used to constrain the inverse model.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Geologic carbon storage in deep saline aquifers is a promising technology for reducing anthropogenic emissions into the atmosphere. Dissolution of injected CO2 into resident brines is one of the primary trapping mechanisms generally considered necessary to provide long-term storage security. Given that diffusion of CO2 in brine is woefully slow, convective dissolution, driven by a small increase in brine density with CO2 saturation, is considered to be the primary mechanism of dissolution trapping. Previous studies of convective dissolution have typically only considered the convective process in the single-phase region below the capillary transition zone and have either ignored the overlying two-phase region where dissolution actually takes place or replaced it with a virtual region with reduced or enhanced constant permeability. Our objective is to improve estimates of the long-term dissolution flux of CO2 into brine by including the capillary transition zone in two-phase model simulations. In the fully two-phase model, there is a capillary transition zone above the brine-saturated region over which the brine saturation decreases with increasing elevation. Our two-phase simulations show that the dissolution flux obtained by assuming a brine-saturated, single-phase porous region with a closed upper boundary is recovered in the limit of vanishing entry pressure and capillary transition zone. For typical finite entry pressures and capillary transition zone, however, convection currents penetrate into the two-phase region. This removes the mass transfer limitation of the diffusive boundary layer and enhances the convective dissolution flux of CO2 more than 3 times above the rate assuming single-phase conditions.

Fracture and fragmentation are extremely nonlinear multiscale processes in which microscale damage mechanisms emerge at the macroscale as new fracture surfaces. Numerous numerical methods have been developed for simulating fracture initiation, propagation, and coalescence. Here, we present a computational approach for modeling pervasive fracture in quasi-brittle materials based on random close-packed Voronoi tessellations. Each Voronoi cell is formulated as a polyhedral finite element containing an arbitrary number of vertices and faces. Fracture surfaces are allowed to nucleate only at the intercell faces. Cohesive softening tractions are applied to new fracture surfaces in order to model the energy dissipated during fracture growth. The randomly seeded Voronoi cells provide a regularized discrete random network for representing fracture surfaces. The potential crack paths within the random network are viewed as instances of realizable crack paths within the continuum material. Mesh convergence of fracture simulations is viewed in a weak, or distributional, sense. The explicit facet representation of fractures within this approach is advantageous for modeling contact on new fracture surfaces and fluid flow within the evolving fracture network. Applications of interest include fracture and fragmentation in quasi-brittle materials and geomechanical applications such as hydraulic fracturing, engineered geothermal systems, compressed-air energy storage, and carbon sequestration.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.