Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

РАСЧЕТ ТЕМПЕРАТУРНОГО РЕЖИМА МНОГОЛЕТНЕМЕРЗЛЫХ ГРУНТОВ С УЧЕТОМ РАДИАЦИОННОГО ИЗЛУЧЕНИЯ В ИНФРАКРАСНОЙ ОБЛАСТИ СПЕКТРА И ТЕРМИЧЕСКОГО СОПРОТИВЛЕНИЯ СНЕЖНОГО ПОКРОВА

Марков Е.В. 1 Пульников С.А. 1 Гербер А.Д. 2
1 ФГБОУ ВО «Тюменский государственный нефтегазовый университет»
2 ООО НИЦ «Геотехэкспертиза»
Выявлено, что методика расчета температурного режима многолетнемерзлых грунтов методом численного решения дифференциального уравнения теплопроводности, закрепленная в нормативно-технической документации РФ, не учитывает радиационный теплообмен в инфракрасной области спектра. Анализ научно-технической литературы показал, что неучёт данного вида теплообмена приводит к существенному искажению результатов моделирования теплового взаимодействия подземных сооружений с многолетнемерзлыми грунтами. Выведено уравнение для задания граничных условий третьего рода на поверхности грунта, учитывающее теплообмен в инфракрасной области спектра и задающее снежный покров методом эквивалентного термического сопротивления, доказана единственность его решения. Показано, что теплообмен в инфракрасной области спектра вносит существенный вклад в температурный режим многолетнемерзлых грунтов и его учет необходим для обеспечения высокой степени достоверности результатов моделирования.
многолетнемерзлый грунт
граничные условия третьего рода
инженерные изыскания
подстилающая поверхность грунта
теплообмен в инфракрасной области спектра
снежный покров
эквивалентное термическое сопротивление
теорема Штурма
1. Алескерова З.Ш., Пульников С.А., Сысоев Ю.С., Казакова Н.В. Оценка эффективности геотехнического мониторинга магистральных газопроводов по качеству продуцируемой информации // Известия высших учебных заведений. Нефть и газ. – 2015. – № 3. – С. 81–86.
2. Горковенко А.И., Вагнер В.В., Пульников С.А. Результаты исследований взаимодействий вибронагруженных магистральных подземных газопроводов с окружающими грунтами // Известия высших учебных заведений. Нефть и газ. – 2007. – № 4. – С. 73–77.
3. Карнаухов М.Ю., Лазарев С.А., Пульников С.А., Сысоев Ю.С. Исследование процесса формирования эксплуатационного положения протяженных участков МГ Уренгой – Челябинск в сложных гидрогеологических условиях // Газовая промышленность. – 2015. – № S724 (724). – С. 53–58.
4. Кушнир С.Я., Карнаухов М.Ю., Пульников С.А., Сысоев Ю.С. Анализ пространственных перемещений магистральных газопроводов с определением граничных зон // Известия высших учебных заведений. Нефть и газ. – 2011. – № 5. – С. 72–75.
5. Кушнир С.Я., Пульников С.А., Малюшин Н.А., Сенив Д.М. Пространственная устойчивость подземных магистральных трубопроводов // Проблемы сбора, подготовки и транспорта нефти и нефтепродуктов. – 2012. – № 1. – С. 51–56.
6. Кушнир С.Я., Пульников С.А., Серебренников А.А., Сенив Д.М. Устойчивость призм обвалования подземных магистральных нефтегазопроводов // Известия высших учебных заведений. Нефть и газ. – 2010. – № 4. – С. 65–70.
7. Кушнир С.Я., Пульников С.А., Сысоев Ю.С., Карнаухов М.Ю. Аналитическая задача определения удлинения газопровода в области аркообразования // Проблемы сбора, подготовки и транспорта нефти и нефтепродуктов. – 2011. – № 4. – С. 74–80.
8. Кушнир С.Я., Пульников С.А., Сысоев Ю.С. Пространственная устойчивость подземного магистрального газопровода на обводненных участках трассы // Известия высших учебных заведений. Нефть и газ. – 2012. – № 1. – С. 72–76.
9. Марков Е.В., Пульников С.А., Гербер А.Д. Проблемы задания граничных условий при моделировании теплового взаимодействия «горячих» трубопроводов с многолетнемерзлыми грунтами // Фундаментальные исследования. – 2015. – № 2–10. – С. 2106–2110.
10. Марков Е.В., Пульников С.А., Сысоев Ю.С. Анализ типовых методик расчета глубины протаивания вечномерзлых грунтов под трубопроводами // Современные проблемы науки и образования. – 2015. – № 1. – С. 142.
11. Мерданов Ш.М., Егорова Т.Ю., Киселев П.В. Создание реологического уравнения состояния снега // Проблемы функционирования систем транспорта. – 2010. – С. 229.
12. Пульников С.А. Взаимодействие вибронагруженных магистральных газопроводов с окружающими грунтами: дис. ... канд. техн. наук. – Тюмень, 2007. – 173 с.
13. Пульников С.А., Сысоев Ю.С., Гербер А.Д., Карнаухов М.Ю. Методика обработки данных мониторинга пространственного положения эксплуатируемого участка подземного магистрального газопровода // Проблемы сбора, подготовки и транспорта нефти и нефтепродуктов. – Уфа: Изд-во ИПТЭР, 2013. – Вып. 3 (93). – С. 58–66.
14. РСН 67-87. Инженерные изыскания для строительства. Составление прогноза измерений температурного режима вечномерзлых грунтов численными методами. – 1987. – 40 с.
15. Тарасенко А.А., Сильницкий П.Ф., Тарасенко Д.А. Проблемы проектирования теплоизоляции оснований вертикальных стальных цилиндрических резервуаров, возводимых в криолитозоне // Труды десятой международной конференции по мерзлотоведению TICOP. – 2012. – С. 300–301.

Заполярные территории относятся к районам со сложными природно-климатическими и инженерно-геологическими условиями. Корректный расчет температурного режима грунтов криолитозоны необходим для обеспечения требуемого уровня надежности сооружений. Наибольшую опасность представляют термически неустойчивые просадочные многолетнемерзлые грунты (ММГ) с высокой влажностью [1–8; 11–13; 15].

В нормативно-технической документации РФ утверждена методика расчета температурного режима ММГ, которая учитывает влияние температурного режима воздуха, солнечной радиации, скорости ветра и снежного покрова, но не учитывает теплообмен подстилающей поверхности грунта с атмосферой в инфракрасной области спектра [14]. Проведенные ранее исследования показывают, что инфракрасное излучение оказывает существенное влияние на результаты прогнозных теплотехнических расчетов [9, 10].

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

Задачи:

– проанализировать исходные данные, необходимые для выполнения расчета температурного режима ММГ;

– вывести формулу для задания граничных условий третьего рода на поверхности грунта;

– обосновать необходимость введения поправки на теплообмен посредством инфракрасного излучения.

Моделирование температурного режима ММГ (рисунок) выполняется путем численного решения дифференциального уравнения теплопроводности методом конечных элементов или конечных разностей:

markov01.wmf (1)

где ρ(T) – плотность грунта; T – температура; cр(T) – теплоемкость грунта; t – время; K(T) – теплопроводность грунта.

На нижней и боковых поверхностях моделируемого массива грунта выбирается граничное условие типа «нулевой тепловой поток». На поверхности грунта задается граничное условие:

markov02.wmf (2)

где Q – тепловой поток от грунта в атмосферу; Qк – конвективный тепловой поток; Qс – поток тепла от солнечной радиации; Qр – радиационный поток тепла в ИК области спектра; a – коэффициент конвективной теплоотдачи подстилающей поверхности грунта; qi – поток тепла от солнечной радиации; εi – степень черноты подстилающей поверхности грунта; σ0 – постоянная Стефана – Больцмана; bi – коэффициент излучения атмосферы; индекс i относится к i-му месяцу в году.

Выражение (2) полностью описывает тепловой поток в периоды времени, когда снежный покров отсутствует. Однако учет снежного покрова методом эквивалентного термического сопротивления усложняет задачу. Появляется дополнительная неизвестная величина – температура поверхности снега. Далее мы займемся выводом формулы, задающей тепловой поток от поверхности грунта в атмосферу через снежный покров.

pic_64.tif

Схема для моделирования температурного режима многолетнемерзлых грунтов и основные граничные условия: Hs – высота снежного покрова; Hгр – глубина моделируемого грунтового основания, соответствующая границе нулевой годовой амплитуды температуры; Lгр – ширина моделируемого грунтового основания; Tпов – температура подстилающей поверхности грунта; Tв – температура атмосферного воздуха; Tгр – температура грунта; ИГЭ-1, 2, 3 – инженерно-геологические элементы

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

Величина теплового потока через снежный покров в случае стационарного теплового потока имеет следующий вид (3):

markov03.wmf (3)

где λi – коэффициент теплопроводности снега; δi – толщина снежного покрова.

Теперь мы имеем систему из двух уравнений (2) и (3) и две неизвестные величины: Q и Tпов. Приравняем уравнения (2) и (3) друг к другу, исключив из рассмотрения величину теплового потока Q:

markov04.wmf (4)

Как видно, уравнение (4) является полиномиальным, четвертой степени относительно неизвестной переменной Tпов. Сгруппируем коэффициенты при одинаковых степенях в формуле (4):

markov05.wmf (5)

Заметим, что все числовые множители мономов в уравнении (5) являются положительными. В соответствии с физическим смыслом задачи нам необходимо найти температуру подстилающей поверхности грунта (температуру поверхности снега в зимний период времени), которая может принимать значения в диапазоне от 0 до ∞.

В связи с громоздкостью окончательного выражения для Tпов мы не будем получать решение в явном виде. Мы лишь докажем, что решение уравнения (5) существует и единственно для данной задачи. Для доказательства применим теорему Штурма. Найдем все полиномы штурма и запишем число смены знаков в полиномах штурма для исследуемых интервалов (табл. 1).

Зная число смены знаков в полиномах Штурма, определим количество действительных положительных корней уравнения N:

N = w(inf) – w(0) = 2 – 1 = 1. (6)

Таким образом, доказано существование единственного положительного корня уравнения (5), а значит и однозначная разрешимость системы уравнений (2)–(3). Благодаря тому, что степень уравнения (5) не превышает четырех, становится возможным найти аналитическое решение, например, по формулам Декарта – Эйлера. Решение для Tпов подставляется в уравнение (3) или (4) и находится величина теплового потока Q.

Теперь рассмотрим влияние снежного покрова на соотношение тепловых потоков различного рода. Исходные данные и результаты расчетов приведены в табл. 2.

Таблица 1

Определение числа смены знаков в полиномах Штурма

Граница диапазона

Полиномы штурма

0

inf

markov06.wmf

+

markov07.wmf

+

+

markov08.wmf

+

markov09.wmf

Число смены знаков, w

2

1

Таблица 2

Оценочный расчет слагаемых конвективного и радиационного теплообмена

Основные виды теплообмена

Последовательность расчета

Конвективный тепловой поток

ai (Tпов – Tв, i)

Поток тепла от солнечной радиации qi

Тепловой поток ИК излучения

markov10.wmf

Исходные данные

ai = 30 Вт/(К∙м2)

(Tгр – Tв, i) = 5 К

qi = 0 Вт/м2 (полярная ночь)

εi = 0,9 (для снега)

σ0 = 5,67∙10–8 Вт/(м2∙К4)

Tв, i = 248 К

Tгр = 253 К

bi = 0,89

Термическое сопротивление снежного покрова, δi/λi, м2/Вт

0,1…5

Температура поверхности снега, К

248,66…247,39

Величина теплового потока от грунта в атмосферу, Вт/м2

20,01…–18,23

0

23,32…19,34

Вклад в суммарный тепловой поток, %

46,2…48,51)

0

53,8...51,51)

Примечание. 1 – вклад в суммарный тепловой поток вычислен от суммы модулей величин конвективного и радиационного теплового потока.

Из табл. 2 видно, что абсолютная величина теплового потока от инфракрасного излучения почти всегда преобладает над конвективным тепловым потоком. Более того, при увеличении толщины снежного покрова конвективный тепловой поток меняет свой знак на противоположный. Это обусловлено тем, что радиационная температура атмосферы ниже температуры воздуха. Как следствие, наблюдается переохлаждение приповерхностных слоев снега до температур, более низких, чем температура атмосферного воздуха. Такой эффект может быть обнаружен только в том случае, если в расчете учитывается радиационный тепловой поток в инфракрасной области спектра по закону Стефана – Больцмана. Его учет является обязательным условием достоверного описания процесса теплопередачи от грунта в атмосферу. Это является очередным подтверждением выводов работ [9, 10].

Выводы

Выявлено, что методика расчета температурного режима ММГ, закрепленная в нормативно-технической документации РФ, не учитывает радиационный теплообмен в инфракрасной области спектра. Анализ научно-технической литературы показал, что неучёт данного вида теплообмена приводит к существенному искажению результатов моделирования теплового взаимодействия подземных сооружений с ММГ [9, 10]. Выведена система уравнений для задания граничных условий третьего рода на поверхности грунта, учитывающая теплообмен в инфракрасной области спектра и задающая снежный покров методом эквивалентного термического сопротивления. Доказана единственность решения системы уравнений для данной задачи. Показано, что теплообмен в инфракрасной области спектра вносит существенный вклад в температурный режим ММГ и его учет необходим для обеспечения высокой степени достоверности результатов моделирования.

Рецензенты:

Якубовский Ю.Е., д.т.н., заведующий кафедрой «Прикладная механика», ФГБОУ ВО «Тюменский государственный нефтегазовый университет», г. Тюмень;

Тарасенко А.А., д.т.н., профессор кафедры ТУР, ФГБОУ ВО «Тюменский государственный нефтегазовый университет», г. Тюмень.


Библиографическая ссылка

Марков Е.В., Пульников С.А., Гербер А.Д. РАСЧЕТ ТЕМПЕРАТУРНОГО РЕЖИМА МНОГОЛЕТНЕМЕРЗЛЫХ ГРУНТОВ С УЧЕТОМ РАДИАЦИОННОГО ИЗЛУЧЕНИЯ В ИНФРАКРАСНОЙ ОБЛАСТИ СПЕКТРА И ТЕРМИЧЕСКОГО СОПРОТИВЛЕНИЯ СНЕЖНОГО ПОКРОВА // Фундаментальные исследования. – 2015. – № 11-1. – С. 100-104;
URL: http://www.fundamental-research.ru/ru/article/view?id=39291 (дата обращения: 20.04.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074