Modelling, vol.6, no.3, 2025 (Scopus)
Underground production and injection operations result in mechanical compaction and mineral chemical reactions that alter porosity and permeability. These changes impact the flow and, eventually, the long-term sustainability of reservoirs utilized for CO2 sequestration and geothermal energy. Even though mechanical and chemical deformations in rocks take place at the pore scale, it is important to investigate their impact at the continuum scale. Rock deformation can be examined using intergranular pressure solution (IPS) models, primarily for uniaxial compaction. Because the reaction rate parameters are estimated using empirical methods and the assumption of constant mineral saturation indices, these models frequently overestimate the rates of compaction and strain by several orders of magnitude. This study presents a new THMC algorithm by combining thermo-mechanical computation with a fractal approach and hydrochemical computations using PHREEQC to evaluate the pressure solution. Thermal stress and strain under axisymmetric conditions are calculated analytically by combining a derived hollow circle mechanical structure with a thermal resistance model. Based on the pore scale, porosity and its impact on the overall excessive stress and strain rate in a domain are estimated by applying the fractal scaling law. Relevant datasets from CO2 core flooding experiments are used to validate the proposed approach. The comparison is consistent with experimental findings, and the novel analytical method allows for faster inspection compared to numerical simulations.