Рассматриваются вопросы математического моделирования трехфазной фильтрации несмешивающихся флюидов - нефти, газа и воды. Системе дифференциальных уравнений сохранения ставится в соответствие полностью явная численная схема на суперэлементной расчетной сетке. Предлагается методика учета работы добывающих скважин при забойном (и пластовом) давлении ниже давления насыщения в рамках суперэлементной модели (при больших размерах ячеек и существенном изменении давления внутри ячейки). Получены удовлетворительные результаты при сопоставительном тестировании с коммерческим программным продуктом.