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 CO2 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. 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 CO2 injection were used to constrain the inverse model. Forward simulation results with representative parameter values from the literature match both low and high surface uplifts reasonably well and predicted the maximum change in the pore pressure. In particular, forward modeling results with estimated Biot's coefficients for reservoir and caprock layers, match the observed uplift well, highlighting the significance of Biot's coefficient in coupled reservoir and geomechanical models. Parameter estimation with 12 parameter sets for both low and high uplift data demonstrates that multiple sets of parameters can match the observed data equally well and the inclusion of the pore pressure data is critically important to constrain the parameter solution during inverse modeling. For a majority of cases, estimation results for both low and high uplift data show the vertical intrinsic permeability and Young's modulus of the reservoir remained close to 13 mD (1.3×10−14 m2) and 10 GPa, respectively, suggesting that these parameters may represent the actual effective properties. Additionally, higher correlations between reservoir permeability and caprock's Biot's coefficient with high surface uplift data were observed consistently under the pore pressure constraint, suggesting the inclusion of the pore pressure constraint is required to estimate the proper values of coupled flow and geomechanical properties associated with different surface uplift data. Overall, this study suggests that given limited data, including Biot's coefficient, in addition to permeability and Young's modulus can enhance parameter estimation of the geomechanical response during GCS.
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.
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. As a result, 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.