Imaging techniques for the analysis of porous structures have revolutionized our ability to quantitatively characterize geomaterials. Digital representations of rock from CT images and physics modeling based on these pore structures provide the opportunity to further advance our quantitative understanding of fluid flow, geomechanics, and geochemistry, and the emergence of coupled behaviors. Additive manufacturing, commonly known as 3D printing, has revolutionized production of custom parts with complex internal geometries. For the geosciences, recent advances in 3D printing technology may be co-opted to print reproducible porous structures derived from CT-imaging of actual rocks for experimental testing. The use of 3D printed microstructure allows us to surmount typical problems associated with sample-to-sample heterogeneity that plague rock physics testing and to test material response independent from pore-structure variability. Together, imaging, digital rocks and 3D printing potentially enables a new workflow for understanding coupled geophysical processes in a real, but well-defined setting circumventing typical issues associated with reproducibility, enabling full characterization and thus connection of physical phenomena to structure. Here we report on our research exploring the possibilities that these technologies can bring to geosciences for coupled multiscale experimental and numerical analysis using 3D printed fractured rock specimens.
For assessing energy-related activities in the subsurface, it is important to investigate the impact of the spatial variability and anisotropy on the geomechanical behavior of shale. The Brazilian test, an indirect tensile-splitting method, is performed in this work, and the evolution of strain field is obtained using digital image correlation. Experimental results show the significant impact of local heterogeneity and lamination on the crack pattern characteristics. For numerical simulations, a phase field method is used to simulate the brittle fracture behavior under various Brazilian test conditions. In this study, shale is assumed to consist of two constituents including the stiff and soft layers to which the same toughness but different elastic moduli are assigned. Microstructural heterogeneity is simplified to represent mesoscale (e.g., millimeter scale) features such as layer orientation, thickness, volume fraction, and defects. The effect of these structural attributes on the onset, propagation, and coalescence of cracks is explored. The simulation results show that spatial heterogeneity and material anisotropy highly affect crack patterns and effective fracture toughness, and the elastic contrast of two constituents significantly alters the effective toughness. However, the complex crack patterns observed in the experiments cannot completely be accounted for by either an isotropic or transversely isotropic effective medium approach. This implies that cracks developed in the layered system may coalesce in complicated ways depending on the local heterogeneity, and the interaction mechanisms between the cracks using two-constituent systems may explain the wide range of effective toughness of shale reported in the literature.
Dissolved CO2 in the subsurface resulting from geological CO2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO2 seeps from a natural CO2 sequestration analog, Crystal Geyser, Utah. Observations along the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO2 seeps.
Shales have small pore and throat sizes ranging from nano to micron scales, low porosity and limited permeability. The poor permeability and complex pore connectivity of shales pose technical challenges to (a) understanding flow and transport mechanisms in such systems and, (b) in predicting permeability changes under dynamic saturation conditions. This study presents quantitative experimental evidence of the migration of water through a generic shale core plug using micro CT imaging. In addition, in-situ measurements of gas permeability were performed during counter-current spontaneous imbibition of water in nano-darcy permeability Marcellus and Haynesville core plugs. It was seen that water blocks severely reduced the effective permeability of the core plugs, leading to losses of up to 99.5% of the initial permeability in experiments lasting 30 days. There was also evidence of clay swelling which further hindered gas flow. When results from this study were compared with similar counter-current gas permeability experiments reported in the literature, the initial (base) permeability of the rock was found to be a key factor in determining the time evolution of effective gas permeability during spontaneous imbibition. With time, a recovery of effective permeability was seen in the higher permeability rocks, while becoming progressively detrimental and irreversible in tighter rocks. These results suggest that matrix permeability of ultra-tight rocks is susceptible to water damage following hydraulic fracturing stimulation and, while shut-in/soaking time helps clearing-up fractures from resident fluid, its effect on the adjacent matrix permeability could be detrimental.
Dissolved CO2 in the subsurface resulting from geological CO2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This work is motivated by observations of CO2 seeps from a natural CO2 sequestration analog, Crystal Geyser, Utah. Observations along the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. The functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO2 seeps.
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.
Multiphase flow in capillary regimes is a fundamental process in a number of geoscience applications. The ability to accurately define wetting characteristics of porous media can have a large impact on numerical models. In this paper, a newly developed automated three-dimensional contact angle algorithm is described and applied to high-resolution X-ray microtomography data from multiphase bead pack experiments with varying wettability characteristics. The algorithm calculates the contact angle by finding the angle between planes fit to each solid/fluid and fluid/fluid interface in the region surrounding each solid/fluid/fluid contact point. Results show that the algorithm is able to reliably compute contact angles using the experimental data. The in situ contact angles are typically larger than flat surface laboratory measurements using the same material. Wetting characteristics in mixed-wet systems also change significantly after displacement cycles.
Four sets of nonreactive solute transport experiments were conducted with micromodels. Each set consisted of three experiments with one variable, i.e., flow velocity, grain diameter, pore-aspect ratio, and flow-focusing heterogeneity. The data sets were offered to pore-scale modeling groups to test their numerical simulators. Each set consisted of two learning experiments, for which all results were made available, and one challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the transverse dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice Boltzmann (LB) approach, and one used a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used the learning experiments to appropriately discretize the spatial grid representations. For the continuum modeling, the required dispersivity input values were estimated based on published nonlinear relations between transverse dispersion coefficients and Peclet number. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values, resulting in reduced dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models, which account for the micromodel geometry and underlying flow and transport physics, needed up to several days on supercomputers to resolve the more complex problems.
Yoon, Hongkyu; Lee, Jonghyun; Kitanidis, Peter K.; Werth, Charles J.; Valocchi, Albert J.
Characterizing subsurface properties is crucial for reliable and cost-effective groundwater supply management and contaminant remediation. With recent advances in sensor technology, large volumes of hydrogeophysical and geochemical data can be obtained to achieve high-resolution images of subsurface properties. However, characterization with such a large amount of information requires prohibitive computational costs associated with “big data” processing and numerous large-scale numerical simulations. To tackle such difficulties, the principal component geostatistical approach (PCGA) has been proposed as a “Jacobian-free” inversion method that requires much smaller forward simulation runs for each iteration than the number of unknown parameters and measurements needed in the traditional inversion methods. PCGA can be conveniently linked to any multiphysics simulation software with independent parallel executions. In this paper, we extend PCGA to handle a large number of measurements (e.g., 106 or more) by constructing a fast preconditioner whose computational cost scales linearly with the data size. For illustration, we characterize the heterogeneous hydraulic conductivity (K) distribution in a laboratory-scale 3-D sand box using about 6 million transient tracer concentration measurements obtained using magnetic resonance imaging. Since each individual observation has little information on the K distribution, the data were compressed by the zeroth temporal moment of breakthrough curves, which is equivalent to the mean travel time under the experimental setting. Only about 2000 forward simulations in total were required to obtain the best estimate with corresponding estimation uncertainty, and the estimated K field captured key patterns of the original packing design, showing the efficiency and effectiveness of the proposed method.
Singh, Rajveer; Yoon, Hongkyu; Sanford, Robert A.; Katz, Lynn; Fouke, Bruce W.; Werth, Charles J.
We investigated the ability of Pseudomonas stutzeri strain DCP-Ps1 to drive CaCO3 biomineralization in a microfluidic flowcell (i.e., micromodel) that simulates subsurface porous media. Results indicate that CaCO3 precipitation occurs during NO3– reduction with a maximum saturation index (SIcalcite) of ~1.56, but not when NO3– was removed, inactive biomass remained, and pH and alkalinity were adjusted to SIcalcite ~ 1.56. CaCO3 precipitation was promoted by metabolically active cultures of strain DCP-Ps1, which at similar values of SIcalcite, have a more negative surface charge than inactive strain DCP-Ps1. A two-stage NO3– reduction (NO3– → NO2– → N2) pore-scale reactive transport model was used to evaluate denitrification kinetics, which was observed in the micromodel as upper (NO3– reduction) and lower (NO2– reduction) horizontal zones of biomass growth with CaCO3 precipitation exclusively in the lower zone. Our model results are consistent with two biomass growth regions and indicate that precipitation occurred in the lower zone because the largest increase in pH and alkalinity is associated with NO2– reduction. CaCO3 precipitates typically occupied the entire vertical depth of pores and impacted porosity, permeability, and flow. This study provides a framework for incorporating microbial activity in biogeochemistry models, which often base biomineralization only on SI (caused by biotic or abiotic reactions) and, thereby, underpredict the extent of this complex process. Furthermore, these results have wide-ranging implications for understanding reactive transport in relevance to groundwater remediation, CO2 sequestration, and enhanced oil recovery.
Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore-scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore-scale experiments of single- and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction-dependent characteristics of granular packs and their impact on fluid flow. A series of drainage and imbibition cycles were conducted on a water-wet, soda-lime glass bead pack under varying confining stress conditions. Simultaneously, X-ray micro-CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro-CT images were segmented using a gradient-based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. Furthermore, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.
Multiscale characteristics of anisotropic, heterogeneous pore structure and compositional (e.g., kerogen, clay, cement, etc) distribution profoundly influence the hydro, mechanical, and chemical response of shale materials during stimulation and production. In this work the impact of these lithologic heterogeneities on physical, chemical, and mechanical properties is investigated over a micron to core scale of shale samples for Cretaceous Mancos Shale. Principal macroscopic lithofacies at a decimeter scale are petrographically examined. Thin sections (∼2-3cm) impregnated with fluorochromes are examined using laser scanning confocal microscopy and optical microscopy with different filters to characterize micro-facies (i.e., texture patterns) and using electron microprobe to identify the mineralogical distribution. Advanced multiscale image analysis for texture classification will be used to identify key features of samples which will be further analyzed using dual focused ion beam-scanning electron microscopy, aberration corrected-scanning TEM and energy dispersive X-ray spectrometry for nano-pore and organic-pore structures and mineralogies at nano scale. This characterization will be examined against experimental data including acoustic emission and nano-indentation measurements of elastic properties using focused ion-beam milled pillars. Finally, multiscale 3-D image stacks will be segmented to rigorously test the scale of a representative elementary volume based on multiple measures from image analysis and pore-scale simulations.
We develop a 2-D pore scale model of coupled fluid flow, reactive transport, and calcium carbonate (CaCO 3) precipitation and dissolution. The model is used to simulate transient experimental results of CaCO 3 precipitation and dissolution under supersaturated conditions in a microfluidic pore network (i.e., micromodel) in order to improve understanding of coupled reactive transport systems perturbed by geological CO 2 injection. In the micromodel, precipitation is induced by transverse mixing along the centerline in pore bodies. The reactive transport model includes the impact of pH upon carbonate speciation and a CaCO 3 reaction rate constant, the effect of changing reactive surface area upon the reaction, and the impact of pore blockage from CaCO 3 precipitation on diffusion and flow. Overall, the pore scale model qualitatively captured the precipitate morphology, precipitation rate, and maximum precipitation area using parameter values from the literature. In particular, we found that proper estimation of the effective diffusion coefficient (D eff) and the reactive surface area is necessary to adequately simulate precipitation and dissolution rates. In order to match the initial phase of fast precipitation, it was necessary to consider the top and bottom of the micromodel as additional reactive surfaces. In order to match a later phase when dissolution occurred, it was necessary to increase the dissolution rate compared to the precipitation rate, but the simulated precipitate area was still higher than the experimental results after ∼30 min, highlighting the need for future study. The model presented here allows us to simulate and mechanistically evaluate precipitation and dissolution of CaCO 3 observed in a micromodel pore network. This study leads to improved understanding of the fundamental physicochemical processes of CaCO 3 precipitation and dissolution under far-from-equilibrium conditions. Copyright 2012 by the American Geophysical Union.
This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers.
Mercury intrusion porosimetry (MIP) is an often-applied technique for determining pore throat distributions and seal analysis of fine-grained rocks. Due to closure effects, potential pore collapse, and complex pore network topologies, MIP data interpretation can be ambiguous, and often biased toward smaller pores in the distribution. We apply 3D imaging techniques and lattice-Boltzmann modeling in interpreting MIP data for samples of the Cretaceous Selma Group Chalk. In the Mississippi Interior Salt Basin, the Selma Chalk is the apparent seal for oil and gas fields in the underlying Eutaw Fm., and, where unfractured, the Selma Chalk is one of the regional-scale seals identified by the Southeast Regional Carbon Sequestration Partnership for CO2 injection sites. Dual focused ion - scanning electron beam and laser scanning confocal microscopy methods are used for 3D imaging of nanometer-to-micron scale microcrack and pore distributions in the Selma Chalk. A combination of image analysis software is used to obtain geometric pore body and throat distributions and other topological properties, which are compared to MIP results. 3D data sets of pore-microfracture networks are used in Lattice Boltzmann simulations of drainage (wetting fluid displaced by non-wetting fluid via the Shan-Chen algorithm), which in turn are used to model MIP procedures. Results are used in interpreting MIP results, understanding microfracture-matrix interaction during multiphase flow, and seal analysis for underground CO2 storage.
Dissolved CO2 from geological CO2 sequestration may react with dissolved minerals in fractured rocks or confined aquifers and cause mineral precipitation. The overall rate of reaction can be limited by diffusive or dispersive mixing, and mineral precipitation can block pores and further hinder these processes. Mixing-induced calcite precipitation experiments were performed by injecting solutions containing CaCl2 and Na2CO3 through two separate inlets of a micromodel (1-cm x 2-cm x 40-microns); transverse dispersion caused the two solutions to mix along the center of the micromodel, resulting in calcite precipitation. The amount of calcite precipitation initially increased to a maximum and then decreased to a steady state value. Fluorescent microscopy and imaging techniques were used to visualize calcite precipitation, and the corresponding effects on the flow field. Experimental micromodel results were evaluated with pore-scale simulations using a 2-D Lattice-Boltzmann code for water flow and a finite volume code for reactive transport. The reactive transport model included the impact of pH upon carbonate speciation and calcite dissolution. We found that proper estimation of the effective diffusion coefficient and the reaction surface area is necessary to adequately simulate precipitation and dissolution rates. The effective diffusion coefficient was decreased in grid cells where calcite precipitated, and keeping track of reactive surface over time played a significant role in predicting reaction patterns. Our results may improve understanding of the fundamental physicochemical processes during CO2 sequestration in geologic formations.