In this paper we propose and study a subgrid model for linear convection-diffusion-reaction problems with fractal rough coefficients. The subgrid model is based on extrapolation model of a modeling residual from coarser scales using a computed solution without subgrid model on a finest scale as reference. We present a priori and a posteriori error estimates, and we show in experiments that a solution with subgrid model on a scale h corresponds to a solution without subgrid model on a scale less than h/4