Analysis and forecasting of changes in the soil temperature distribution in the area of the city of Norilsk

• Main Authors: Fedotov, A. A., Kaniber, V.V., Khrapov, P. V.
• Format: Article in Journal/Newspaper
• Language: Russian
• Published: International Journal of Open Information Technologies 2020
• Subjects: Norilsk; norilsk; permafrost
• Online Access: http://injoit.org/index.php/j1/article/view/996

The article studies the initial boundary value problem for a non-stationary one-dimensional heat equation that simulates the distribution of soil temperature in the area of Norilsk. The mathematical model is based on solid-liquid phase transitions. To determine the model parameters, data from the meteorological station and reports on engineering and geological surveys were used to obtain the necessary physical and thermophysical characteristics of the calculated area. For the numerical solution of the problem, the finite volume method (FVM) was used. The calculation was started on January 1 of the first year of observation of the soil temperature distribution. Simulation of the soil temperature distribution was carried out until the moment of reaching the non-stationary periodic mode. In order to analyze the temperature field, graphs of the temperature dependence on the depth for January and July of the 7 selected years of observation were constructed. The study of the results showed that it takes about 50 years for the soil temperature to reach a non-stationary periodic mode at a depth of 20 m. For a steady-state periodic mode, temperature dependences on the depth for each month are constructed and the depth of the active layer (0.5 m) and the depth of zero amplitudes (12 m) are found. The forecast of the ground temperature distribution for 2080 was modeled for two scenarios of Representative Concentration Pathway (RCP) of global warming: moderate RCP2.6 and negative RCP8.5. The RCP2.6 scenario showed an increase in the depth of thawing (freezing) by 0.1 m and a decrease in the depth of zero amplitudes by 1 m, as well as an increase of soil temperature at all depths by an average of 2°C. The results of calculations for the RCP8.5 scenario showed an increase in the depth of the active layer by 0.5 m (2 times compared to the present time) and a decrease in the depth of zero amplitudes by 4 m, the soil temperature increased by 5.7°C on average. In connection with the fuel spill that occurred in May 2020 from a storage tank in Norilsk due to subsidence of the foundation mainstays, the soil temperature distribution was modeled in 1980. The results showed an increase in the depth of the active layer since 1980 by 0.1 m, which could cause subsidence of the foundation piles of the storage tank. Further, a comparison of the long-term average monthly temperatures by decades was made, in which it was found that 2010-2020 has the largest number of the highest values of average monthly temperatures. This study further confirms the degradation of permafrost В работе проведено исследование начально-краевой задачи для нестационарного одномерного уравнения теплопроводности, моделирующей распределение температуры грунта в районе г. Норильска. Математическая модель построена с учетом фазовых переходов – твердое тело-жидкость. Для определения параметров модели использованы данные метеорологической станции и отчеты об инженерно-геологических изысканиях, с помощью которых получены необходимые физические и теплофизические характеристики расчетной области. Для численного решения задачи использовался метод контрольного объема. За начало расчетов принято 1 января первого года наблюдения за температурным режимом грунта. Моделирование температурного режима грунта проводилось до момента выхода на нестационарный периодический режим. С целью анализа изменения температурного поля грунта были построены графики зависимости температуры от глубины для января и июля 7 выбранных годов наблюдения. Исследование результатов показало, что для выхода температуры грунта на глубине до 20 м на нестационарный периодический режим необходимо около 50 лет. Для установившегося периодического режима построены зависимости температуры от глубины для каждого месяца и найдена глубина активного слоя (0,5 м) и глубина нулевых амплитуд (12 м). Смоделирован прогноз температурного режима грунта на 2080 год по двум сценариям Representative Concentration Pathway (RCP) глобального потепления: умеренному RCP2.6 и негативному RCP8.5. Сценарий RCP2.6 показал увеличение глубины оттаивания (промерзания) на 0,1 м и уменьшение глубины нулевых амплитуд на 1 м, а также увеличение температуры по глубине в среднем на 2°C. Результаты расчетов по сценарию RCP8.5 показали увеличение глубины активного слоя на 0,5 м (в 2 раза по сравнению с настоящим временем) и уменьшение глубины нулевых амплитуд на 4 м, температура грунта в среднем выросла на 5,7°C. В связи с произошедшим в мае 2020 года разливом топлива из резервуара в г. Норильске из-за проседания опор фундамента был смоделирован температурный режим грунта в 1980 году. Результаты показали увеличение глубины активного слоя с 1980 года на 0,1 м, которое и могло вызвать просадку свай фундамента резервуара. Далее было проведено сравнение многолетних среднемесячных температур по десятилетиям, в котором было установлено, что 2010-2020 годы имеют наибольшее количество самых высоких значений среднемесячных температур. Данное исследование дополнительно подтверждает деградацию вечной мерзлоты