En este trabajo se estudió numéricamente el flujo agua-petróleo en una medio poroso heterogéneo, con énfasis en la interface entre dos capas homogéneas de la matriz porosa. La heterogeneidad considerada consistió en una función discontinua de capilaridad sobre la interface. La ecuación diferencial fue resuelta utilizando un esquema tipo upwind completamente implícito basado en el método de volúmenes finitos. La incógnita principal fue la saturación del agua. Se evaluó el impacto sobre la saturación de la interface de diversos parámetros tales como: la porosidad, la saturación inicial de agua, la razón de inyección de agua, la fuerza gravitacional, y el tamaño del grano del material. Se concluye que una completa caracterización de la matriz porosa, junto con una adecuada combinación de parámetros de control y modelado, son condiciones fundamentales para definir condiciones de extracción adecuadas. El algoritmo estudiado posee un gran potencial para ser utilizado en etapas tempranas de diseño y planeación para la extracción del petróleo.