Mohammad Amin Amooie

Critical Dynamics of Gravito-Convective Mixing in Geological Carbon Sequestration

When CO2 is injected in saline aquifers, dissolution causes a local increase in brine density that can cause Rayleigh-Taylor-type gravitational instabilities. Depending on the Rayleigh number, density-driven flow may mix dissolved CO2 throughout the aquifer at fast advective time-scales through convective mixing. Heterogeneity can impact density-driven flow to different degrees. Zones with low effective vertical permeability may suppress fingering and reduce vertical spreading, while potentially increasing transverse mixing. In more complex heterogeneity, arising from the spatial organization of sedimentary facies, finger propagation is reduced in low permeability facies, but may be enhanced through more permeable facies. The connectivity of facies is critical in determining the large-scale transport of CO2-rich brine. We perform high-resolution finite element simulations of advection-diffusion transport of CO2 with a focus on facies-based bimodal heterogeneity. Permeability fields are generated by a Markov Chain approach, which represent facies architecture by commonly observed characteristics such as volume fractions. CO2 dissolution and phase behavior are modeled with the cubic-plus-association equation-of-state. Our results show that the organization of high-permeability facies and their connectivity control the dynamics of gravitationally unstable flow. We discover new flow regimes in both homogeneous and heterogeneous media and present quantitative scaling relations for their temporal evolution.

Implicit finite volume and discontinuous Galerkin methods for multicomponent flow in unstructured 3D fractured porous media

Abstract We present a new implicit higher-order finite element (FE) approach to efficiently model compressible multicomponent fluid flow on unstructured grids and in fractured porous subsurface formations. The scheme is sequential implicit: pressures and fluxes are updated with an implicit Mixed Hybrid Finite Element (MHFE) method, and the transport of each species is approximated with an implicit second-order Discontinuous Galerkin (DG) FE method. Discrete fractures are incorporated with a cross-flow equilibrium approach. This is the first investigation of all-implicit higher-order MHFE-DG for unstructured triangular, quadrilateral (2D), and hexahedral (3D) grids and discrete fractures. A lowest-order implicit finite volume (FV) transport update is also developed for the same grid types. The implicit methods are compared to an Implicit-Pressure-Explicit-Composition (IMPEC) scheme. For fractured domains, the unconditionally stable implicit transport update is shown to increase computational efficiency by orders of magnitude as compared to IMPEC, which has a time-step constraint proportional to the pore volume of discrete fracture grid cells. However, when lowest-order Euler time-discretizations are used, numerical errors increase linearly with the larger implicit time-steps, resulting in high numerical dispersion. Second-order Crank–Nicolson implicit MHFE-DG and MHFE-FV are therefore presented as well. Convergence analyses show twice the convergence rate for the DG methods as compared to FV, resulting in two to three orders of magnitude higher computational efficiency. Numerical experiments demonstrate the efficiency and robustness in modeling compressible multicomponent flow on irregular and fractured 2D and 3D grids, even in the presence of fingering instabilities.

Dissolution Trapping of Carbon Dioxide in Heterogeneous Aquifers

The geologic architecture in sedimentary reservoirs affects the behavior of density-driven flow and the dispersion of CO2-rich brine. The spatial organization and connectivity of facies types play an important role. Low-permeability facies may suppress fingering and reduce vertical spreading, but may also increase transverse mixing. This is more pronounced when geologic structures create preferential flow pathways through connected facies types. We perform high-resolution simulations of three-dimensional (3D) heterogeneous formations whose connectivity cannot be represented in two-dimensional models consistent with percolation theory. This work focuses on the importance of 3D facies-based heterogeneity and connectivity on advection-diffusion transport of dissolved CO2. Because the dissolution of CO2 and the subsequent density increase of brine are the driving force for gravitational instabilities, we model the phase behavior with the accurate cubic-plus-association equation-of-state, which accounts for the self-association of polar water molecules and the cross-association between CO2 and water. Our results elucidate how the spatial organization of facies affects the dynamics of CO2 convective mixing. Scaling relations for the evolution of a global dispersion-width provide insights that can be universally applied. The results suggest that the long-term evolution and scaling of dispersion are surprisingly similar for homogeneous and (binary and multiscale) heterogeneous porous media.

Hydrothermodynamic mixing of fluids across phases in porous media

We investigate the coupled dynamics of fluid mixing and viscously unstable flow under both miscible (single-phase) and partially miscible (two-phase) conditions, and in both homogeneous and heterogeneous porous media. Higher-order finite element methods and fine grids are used to resolve the small-scale onset of fingering and tip splitting. An equation of state determines the thermodynamic phase behavior and Fickian diffusion. We compute global quantitative measures of the spreading and mixing of a diluting slug to elucidate key differences between miscible and partially miscible systems. Hydrodynamic instabilities are the main driver for mixing in miscible flow. In partially miscible flow, however, we find that relative permeabilities spread the two-phase zone. Within this mixing zone dissolution and evaporation drive mixing thermodynamically while reducing mobility contrasts and thus fingering instabilities. The different mixing dynamics in systems involving multiple phases with mutual solubilities have important implications in hydrogeology and energy applications, such as geological carbon sequestration and gas transport in hydrocarbon reservoirs.

Multicomponent competitive monovalent cation exchange in hierarchical porous media with multimodal reactive mineral facies

This paper presents a stochastic model for multicomponent competitive monovalent cation exchange in hierarchical porous media. Reactive transport in porous media is highly sensitive to heterogeneities in physical and chemical properties, such as hydraulic conductivity (K), and cation exchange capacity (CEC). We use a conceptual model for multimodal reactive mineral facies and develop a Eulerian-based stochastic theory to analyze the transport of multiple cations in heterogeneous media with a hierarchical organization of reactive minerals. Numerical examples investigate the retardation factors and dispersivities in a chemical system made of three monovalent cations (Na+, K+, and Cs+). The results demonstrate how heterogeneity influences the transport of competitive monovalent cations, and highlight the importance of correlations between K and CEC. Further sensitivity analyses are presented investigating how the dispersion and retardation of each cation are affected by the means, variances, and integral scales of K and CEC. The volume fraction of organic matter is shown to be another important parameter. The Eulerian stochastic framework presented in this work clarifies the importance of each system parameters on the migration of cation plumes in formations with hierarchical organization of facies types. Our stochastic approach could be used as an alternative to numerical simulations for 3D reactive transport in hierarchical porous media, which become prohibitively expensive for the multicomponent applications considered in this work.