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