For nonlinear coupled system with moving boundary values of multilayer dynamics of fluids in porous media, two kinds of upwind finite difference schemes applicable to parallel arithmetic are put forward.Some techniques, such as the change of variables, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted.The estimates in l2 norm are derived to determine the error in the approximate solution.This method has already been applied to the numerical simulation of migration-accumulation of oil resources.