Аннотация:This work is devoted to developing the self-consistent thermo-hydro-chemo-mechanical reactive transport model to predict and describe natural and industrial petroleum processes at different scales.We develop a version of the front tracking approach for multicomponent multiphase flow in order to treat spontaneous splitting of discontinuities. We revisit the solution for the Riemann problem and systematically classify all possible configurations as functions of initial concentrations on both sides of the discontinuity. We validate the algorithm against finite volume high-resolution technics and high-order spectral finite elements.To calculate the parameters of phase equilibria, we utilize an approach based on the direct minimization of the Gibbs energy of a multicomponent mixture. This method ensures the consistency of the thermodynamic lookup tables. The core of the algorithm is the non-linear free-energy constrained minimization problem, formulated in the form of a linear programming problem by discretization in compositional space.The impact of the complex rheological response of porous matrix on the morphology of fluid flow and shear deformation localization is considered. Channeling of porosity waves and shear bands morphology and their orientation is investigated for viscoelastoplastic both shear and bulk rheologies.