Рассматривается численное моделирование температурно-фильтрационного режима в теле и основании каменно-земляной плотины в условиях вечной мерзлоты. В качестве примера моделируется каменно-земляная плотина Вилюйской ГЭС-1,2 в двумерной постановке. Для аппроксимация по пространству используется метод конечных элементов, по времени - схема расщепления по физическим процессам