We study a two-scale coupled system consisting of a macroscopic elliptic equation and a microscopic parabolic equation. This system models the interplay between a gas and liquid close to equilibrium within a porous medium with distributed microstructures. We use formal homogenization arguments to derive the target system. We start by proving well-posedness and inverse estimates for the two-scale system. We follow up by proposing a Galerkin scheme which is continuous in time and discrete in space, for which we obtain well-posedness, a priori error estimates and convergence rates. Finally, we propose a numerical error reduction strategy by refining the grid based on residual error estimators.