Статья посвящена методам решения одномерной радиальной задачи теплопередачи в окружающие скважину многолетнемерзлые породы (ММП). Приведена математическая модель указанной задачи, для ее решения рассмотрены четыре численных метода: метод энтальпий с использованием явной схемы, метод ловли фронта в узел сетки с использованием неявной схемы, метод ловли фронта в узел сетки с использованием шеститочечной симметричной схемы, и квазистационарный подход. Получено автомодельное решение, позволяющее рассчитывать радиус протаивания в ММП и используемое для оценки вышеприведенных численных методов (в статье приведено сравнение автомодельного решения и результатов расчетов радиуса протаивания по методу энтальпий). Также представлен сравнительный анализ численных решений рассматриваемой задачи теплопередачи в ММП с использованием различных численных методов. Итогами работы являются построенное автомодельное решение для радиальной задачи Стефана и выбор наиболее предпочтительного численного метода, обеспечивающего высокую точность и наибольшую среди рассмотренных методов скорость расчетов радиуса протаивания и распределения температуры в окружающих скважину мерзлых породах.