The combined flooding in oil recovery is a multiphase chromatographic process.Modeling for multiphase

multicomponent nonlinear chromatography was applied to simulation of the surfactant/polymer flooding process.By means of the method of characteristics

the characteristic velocities of saturation and concentration waves were solved.Then

integrating the eigenvectors led to wave profiles of saturation and concentration in physical spatial domain.The simulation results indicated that the concentration shock of polymer appeared earlier than that of surfactant; moreover due to solutes’movement

saturation wave could develop into multi-shock

which enabled the oleic flow to form oil bank rapidly.In addition

the influencing factors

such as concentration

flow velocity

adsorption constant and partitioning coefficient were also discussed.