Important geoscience and environmental applications such as geologic carbon storage, environmental remediation, and unconventional oil and gas recovery are best understood in the context of reactive flow and multicomponent transport in the subsurface environment. The coupling of chemical and microbiological reactions with hydrological and mechanical processes can lead to complex behaviors across an enormous range of spatial and temporal scales. These coupled responses are also strongly influenced by the heterogeneity and anisotropy of the geologic formations. Reactive transport processes can change the pore morphology at the pore scale, thereby leading to nonlinear interactions with advective and diffusive transport, which can strongly influence larger-scale properties such as permeability and dispersion. Therefore, one of the greatest research challenges is to improve our ability to predict these processes across scales (DOE 2007). The development of pore-scale experimental and modeling methods to study reactive processes involving mineral precipitation and dissolution, and biofilm dynamics allows more fundamental investigation of physical behavior so that more accurate and robust upscaled constitutive models can be developed for the continuum scale.

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.