Multiphysics phase-field modeling for thermal cracking and permeability evolution in oil shale matrix during in-situ conversion process

In-situ conversion process (ICP) is a promising recovery technology for economically effective exploitation of unmatured oil shale. High temperature treatment during ICP facilitates the chemical conversion of kerogen content in oil shale, while markedly improving the accessibility of pyrolyzed fluid through thermally-induced microcracks. Although extensively studied experimentally, the microstructure development in shale is not fully understood due to the lack of consideration for multiphysics effects during ICP at the mineral scale. This paper presents a numerical study on thermal cracking and permeability evolution in oil shale matrix during in-situ conversion process. To this end, a multiphysics phase-field model for investigating thermo-mechanical response, chemical reaction, and seepage behavior is developed, while arbitrary crack growth in heterogeneous shale matrix is naturally predicted using the phase-field method for fracture. Implementing the proposed numerical model, thermal cracking of heterogeneous granodiorite is first simulated, from which the numerical outcome regarding crack morphology is reasonably consistent with experimental data. Permeability evolution of oil shale matrix is found to be attributed to early-stage thermal crack propagation and late-stage kerogen decomposition by correlating multiple variables. Anisotropic permeability is also observed and investigated by examining crack morphology, pore-space interconnectivity and fluid permeation. Further analysis reveals that the shale matrix microstructure, in-situ stress and heating temperatures play vital roles in influencing thermal cracking and permeation behaviors of the shale matrix. The results provide a unique perspective to understanding thermal cracking and permeability augmentation of heterogeneous rocks.