Numerical simulation of migration-accumulation of oil resources in porous media is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis puts forward a mathematical model, a careful parallel operator splitting-up implicit iterative scheme, parallel arithmetic program, parallel arithmetic information and alternating-direction mesh subdivision. For the actual situation of Tanhai region of Shengli Petroleum Field, our numerical simulation test results and the actual conditions are coincident. For the model problem (nonlinear coupled system) optimal order estimates in $l^2$ norm are derived to determine the errors. We have successfully solved the difficult problem in the fields of permeation fluid mechanics and petroleum geology.