Abstract: By using explicit saturation and implicit concentration methods,numerical solution is presented to the polymer drive model characterizing dispersion and adsorption of polymer solution flowing through a porous medium.Saturation equation is solved by using an explicit total variation diminishing(TVD)method.To ensure the stability of concentration equation calculation,Crank-Nicolson difference format is used in spatial discretization,and variable is quasi-linear disposed in temporal discretization.The validity of the method presented is verified by comparison with analytic solutions.The calculated example indicates that dispersion causes dilution and dissipation of polymer solutions,and adsorption results in a loss,thus leading to lagging of concentration propagation.Also,the important polymer flooding mechanism- "oil block" is illustrated in the results of calculation.Under the condition of slug injection,the breakthrough time of oil enrichment zone lies between the breakthrough time of polymer concentration front and polymer concentration peak.

Key words: polymer flooding, numerical simulation, enhanced oil recovery