This paper presents a new numerical model of reactive polymer flow in heterogeneous porous media. A moment representation of the log-normal polymer molecular weight distribution is used to model polymer as a multi-component species. Three leading moments are used to simulate the polymer transport and reaction processes in a two-dimensional porous medium. The 2D, multiphase polymer flow model is based on a mass-transport equation for multi-component species and is coupled with kinetic models of the gelation process using an operator splitting scheme. The sensitivity of various parameters and constitutive equations is presented.