В работе представлен численный алгоритм для моделирования двухфазного течения в образцах горных пород, содержащих как открытые поры, так и поры, заполненные микропористым материалом. Для моделирования потоков используется уравнение Навье-Стокса-Бринкмана. Перенос фаз определяется уравнением Кана-Хиллиарда в открытых порах и уравнением Баклея-Леверетта в пористом материале. В работе предложена единая конечно-разностная аппроксимация двух уравнений переноса, удовлетворяющая естественным условиям сопряжения. Однако уравнение Кана-Хиллиарда требует дополнительного граничного условия, которое должно быть выполнено на границе раздела, чтобы обеспечить угол смачивания.