We propose a coupled hydro-chemo-mechanical model and its 1D numerical implementation. We demonstrate its application to model filtration of a multicomponent fluid in deforming and reacting host rocks, considering changes in the densities, phase proportions and chemical compositions of coexisting phases. We presented 1D numerical implementation on the example of soapstone formation from serpentinite during H₂O-СО2 fluid filtration with low concentration of СО2 coupled with viscous deformation of mineral matrix, considering MgO-SiO2-H₂O-СО2 system. The numerical results show porosity wave propagation by viscous (de)compaction mechanism accompanied with the formation of an elongated zone with higher filtration properties. After the formation of such a channel, the formation and propagation of reaction fronts occurs associated with the transformation of the mineral composition of the original rock. During H2O-CO2 fluid filtration, starting from 1 weight percent of dissolved CO2, carbonization of hydrated serpentinite starts, specifically antigorite transforms to magnesite and talc.