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