For the section coupled system of multilayer dynamics of fluids in porous media, a parallel scheme modified by the characteristic finite difference fractional steps is proposed for a complete point set consisting of coarse and fine partitions. Some tech- niques, such as calculus of variations, energy method, twofold-quadratic interpolation of product type, multiplicative commutation law of difference operators, decomposition of high order difference operators, and prior estimates, are used in theoretical analysis. Optimal order estimates in 12 norm are derived to show accuracy of the second order approximation solutions. These methods have been used to simulate the problems of migration-accumulation of oil resources.
