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