В данной работе построен приближенный метод решения уравнения для насыщенности в задаче двухфазной неравновесной фильтрации. Это уравнение относится к уравнению типа конвекции-диффузии с преобладанием конвекции и с дополнительным членом, содержащим производную решения третьего порядка. Из-за гиперболического характера уравнения его решение сопровождается рядом трудностей, которые приводят к необходимости тщательного выбора метода решения. На основе вычислительных экспериментов проведено сравнение трех классических стабилизированных методов конечных элементов (SUPG, GLS и USFEM).
Динамика протекания фильтрационных течений многофазной жидкости нелинейным образом зависит как от структурно-механических свойств жидкости, так и свойств окружающего скелета. Исследование процесса течения многофазной жидкости в пористой среде наиболее полно проведено в предположении о локальном фазовом равновесии. Однако в реальных пластовых условиях существенное влияние на процесс фильтрации имеет свойство запаздывания насыщенности фазы, изучение которого привело к возникновению теории неравновесной фильтрации. Необходимость учета данного явления при разработке нефтяных месторождений обсуждается во многих работах [1, 2]. В настоящей работе рассматривается модель двухфазной неравновесной фильтрации с обобщенным законом неравновесности