In this paper, taking account of the major physical and chemical mechanisms, such as: for polymer, shearing propertypermeability reduction, adsorption, inaccessible porous volumes, for gel, gelation speed, water viscosity changing with gel, permeability reduction, adsorption and retention in reservoir rocks, a three-dimensional, three-phase (oleic, vapor, aqueous) and six-component mathematical model has been established for polymer flooding and profile control. By use of full implicit finite difference method and calling PETSc linear solving system, the full implicit polymer flooding and profile control simulation software has been developed on PC-Linux environment based on black oil simulator, water flooding, polymer flooding and profile control simulation methods are integrated and applied into practice.