Рассмотрены постановка и численное моделирование совместной задачи, включающей динамику образования слоя вечной мерзлоты под воздействием сезонных изменений отрицательной среднегодовой температуры. Отмечено, что предложенная постановка задачи и методика расчета позволяют одновременно определить динамику зоны оттаивания многолетнемерзлых пород и температурный режим добываемой жидкости в любой момент времени.