In this paper, we establish a homogenization result for a doubly nonlinear parabolic system arising from the hygro-thermo-chemical processes in porous media taking into account memory phenomena. We present a mesoscale model of the composite (heterogeneous) material where each component is considered as a porous system and the voids of the skeleton are partially saturated with liquid water. It is shown that the solution of the mesoscale problem is two-scale convergent to that of the upscaled problem as the spatial parameter goes to zero. KEYWORDS coupled transport processes in porous media, homogenization, nonlinear degenerate parabolic system, Robin boundary conditions, two-scale convergence MSC CLASSIFICATION 35K55; 35K65; 35D30; 35B27; 35B40; 35A01Here, represents the averaged mass density of the -phase (eg, solid, liquid water, oil, and gas), J is the mass flux and s stands for the production term. Further, e is the total internal energy of the -phase, q is the heat flux, stands for the volumetric heat source, represents the term expressing energy exchange with the other phases, and the symbol H stands for the specific enthalpy of the -phase.Because of the nonlinear structure, complexity and multi-scale nature of the problems (1) and (2), this system must be solved numerically using suitable computational methods. Nevertheless, the complexity of the microscopic structure of heterogeneous multiphase materials makes detailed numerical simulations very expensive. Therefore, the effective (homogenized) material coefficients are used instead of taking into account properties of individual phases. The Math Meth Appl Sci. 2019;42:6227-6258.wileyonlinelibrary.com/journal/mma