We propose a new iterative numerical scheme designed for massively parallel processing for an immiscible displacement in a naturally fractured reservoir. The procedure is based on a domain decomposition technique applied to a mixed finite element approximation of the problem; the domain is decomposed into individual elements. Numerical experiments are presented to illustrate its performance on a CM-5 system.