Математическая модель распространения в атмосфере парниковых газов, порожденных стационарным точечным источником Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

УДК 502.36

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

ИСТОЧНИКОМ

С.А. ЧЕРНЯВСКИЙ

Статья представлена доктором технических наук, профессором Кузнецовым В.Л.

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

Ключевые слова: турбулентная диффузия, стационарный источник, парниковые газы, абсорбция.

Введение

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

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

1. Постановка задачи

Рассматривается задача о распространении в приземном слое атмосферы парниковых газов, выбрасываемых промышленным предприятием. Относительно источника выбросов экологического загрязнения, моделируемого точечным источником, будем считать известными его «производительность» - I и его пространственное расположение. Ветровое поле V = V (г) будем считать плоскослоистым, постоянным на каждой высоте и не меняющимся во времени.

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

2. Моделирование процесса распространения парниковых газов в атмосфере

Основным уравнением, характеризующим распространение парниковых газов в атмосфере, является локальный закон сохранения вещества, который представлен в виде

— + )+1п = 0. (1)

дt Х 1

Здесь п - концентрация парникового газа; 1 - коэффициент, описывающий уменьшение концентрации

парникового газа за счёт абсорбции; J - плотность потока диффундирующего вещества, представляющего собой сумму потоков диффузии

J = J1 + J 2 (2)

^ = п(х, у, z, t)• V (х, у, z) (3)

J2 = ~(к 1 + К2 )§а(п(х, у, z, t)) (4)

где J1 - поток молекулярной диффузии; J 2 - поток турбулентной диффузии; К1 - коэффициент молеку-

лярной диффузии, который отражает перенос примеси в нижний слой атмосферы; К2 - коэффициент

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

Перенос парникового газа в нижний слой незначителен в отличие от переноса в верхние слои атмосферы, поэтому К1 » 0 .

С учётом подстановок в уравнение (1) соотношений (2) - (4) и сделанных математических преобразований имеем следующее дифференциальное уравнение

дп дп дп дп д 2п д 2п д 2п

----+ V-----+ V-----+ V — = к —- + к —- + к —- -1п ,,ч

дt х дх У ду 2 д2 х дх У ду 2 д2 . (5)

Уравнение (5) представляет собой уравнение турбулентной диффузии. Его решение получим при следующих приближениях.

Ветровое поле стационарно, вектор V = V (г) параллелен поверхности земли и направлен вдоль оси х в выбранной нами системе координат. На первом этапе будем полагать, что осадки отсутствуют (1 = 0). Относительно турбулентности атмосферы предположим её изотропность и однородность, т.е.

к » к » к = к .

х у г

Смоделируем источник экологического загрязнения точечным источником мощности М , расположенного в точке с координатами (0,0, Н) , через которую в атмосферу выбрасывается парниковый газ.

С учётом сделанных допущений уравнение (5), описывающее появление и распространение парниковых газов, имеет вид

— + Ух —---------к Ап = Мб(у)б(2 - Н )ії(ґ)

ді х дх •

(6)

Естественно предположить, что на бесконечном удалении от источника концентрация выброшенных парниковых газов стремится к нулю, т.е. п ® 0 при |г | ® ¥ .

Решения уравнения (6), описывающего эволюцию распространения парникового газа при разовом выбросе имеют вид

^ (г - Н )2 (г+Н )2 ^

п1 (х, у, 2,і ) =

М

п2\ х

(х, у, г, і ) =

8(крі)

М

-в

(х-Ухі)2 -у_ 4 кі 4кі

в

4кі

+ в

4кі

( х-уі )2 -4 кі 4 кі

V

( - )2 2 -

8(кжі)

_( 2 -Н )2 4 кі

- в

_( 2 + Н )2 4 кі

V

(7)

(8)

/

где і - время, прошедшее после выброса.

Решение (7) имеет место при распространении парникового газа над земной, а (8) - над водной поверхностями.

Анализ решения (8) показывает, что водная поверхность хорошо абсорбирует парниковый газ, т.е.

ііш Ш п 2 (х, у, г, і }<іхйуйі =

М

V'

1

1іш { V

і 8(клі)

-(х-У )2

4кі 4кі

Ґ -Ьн£ в 4к + в

(.2 + Н )2

4 кі

ёхёуё2 =0.

В случае непрерывного действия источников, функцию стационарного распределения концентрации парниковых газов в атмосфере можно получить, используя соотношения (7) или (8).

В качестве примера рассмотрим вывод выражения для стационарной концентрации парниковых газов над земной поверхностью.

Нетрудно видеть, что

п

(х, у, 2) = ІI(і,)п(х, у, 2, і - і'}йі' = І п1 (х, у, 2, т)ёт.

в

в

в

Значение интеграла, дающего функцию распределения концентрации парникового газа для случая стационарных выбросов при отсутствии атмосферных осадков, имеет вид

п*(х У, 2) =

М

4лк

(

1

-у/х2 + у2 +(2-Н )2

1

х2 + у2 + (2 - Н)2

-^х2 + у2 +(2+Н )2 2к

і

х2 + у2 + (2 + Н )2

(9)

При наличии атмосферных осадков в нестационарном случае при разовом выбросе функция распространения концентрации парникового газа может быть представлена в виде

п(х, у, 2, t) = п1 (х, у, 2, t)е

п(х, у, 2, і):

М

(х -Ухі )2 + у 2 4кі

С (я -Н )2 (г+Н )2

4кі + в 4кі

-1

8(кжі)

V

Соответствующее стационарное решение пх в случае атмосферных осадков находится по формуле, аналогичной (9)

( , и........»---V . , П------—Г7--, >

1(х y, 2) =

М

4рк

1

-Ч (х2 +у2 +(2-Н)2 )(у

-^х2 + у2 +(2 - Н)2

1

- „АІ*2 +у2+(2+н)2 Xу

■^х2 + у2 +(2+Н)2

(10)

Коэффициент абсорбции X зависит от интенсивности осадков и влияет на уменьшение парниковых газов из атмосферы. В работе в качестве примера осадков взята морось, воздействующая на парниковый газ, интенсивность которой 0,01 с-1.

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

т(ух ) = -

3. Результаты численного исследования

^у^ + 4к1

2к

(11)

В качестве парникового газа выбран С02, поскольку это вещество выбрасывают предприятия всех отраслей промышленности, и оно наиболее сильно влияет на изменение среднегодовой температуры. Для

С02 коэффициент диффузии к = 1,4• 10-5 м2/с при температуре воздуха Т = 0 °С и нормальном давлении.

У 4

—- 2- .....

б 0- <щ||5

1 1 -2- -4-

Рис. 1. Распространение С02 при разных значениях скорости ветра: а - распространение С02 при скорости ветра V,. = 4 м/с; б - распространение С02 при скорости ветра V,. = 17 м/с

На рис. 1 представлены контурные графики распространения С02 в атмосфере при разных значениях скоростей ветра на высоте Н = 2 = 15 м (высота жизненно-важного слоя атмосферы) с заданной производительностью источника I = 10 кг/с. Данные графики при наличии и отсутствии осадков подобны. Это связано с тем, что морось, которая воздействует на С02, не вносит существенного вклада в процесс абсорбции С02.

У

+

в

в

в

в

в

+

в

в

У

а

X

X

На рис. 2 представлена графическая зависимость изменения эффективного коэффициента диффузии (11) с увеличением скорости ветра. При скорости ветра ух > 20 м/с влияние С02 на атмосферу практически не ощутимо.

Заключение

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

диффузии от скорости ветра

ЛИТЕРАТУРА

1. Берлянд М.Е. Современные проблемы атмосферной диффузии и загрязнения атмосферы. - Л.: Гидрометеоиздат, 1975.

2. Монин А. С. Атмосферная диффузия. - М.: Институт океанологии им. П.П. Ширшова РАН, 1959.

3. Степаненко С.Н., Волошин В.Г., Типцов С.В.. Решение уравнения турбулентной диффузии для стационарного и точечного источника. - Одесса: Одесский государственный экологический университет, 2007.

THE MATHEMATICAL MODEL OF THE DISTRIBUTION OF GREENHOUSE GASES IN THE ATMOSPHERE, GENERATED BY THE STATIONARY POINT SOURCE

Chernyavskiy S.A.

Examine a stationary solution of the equation of turbulent diffusion of greenhouse gases in the isothermal atmosphere with a flat-flaky wind field, under the assumption small hilly landscape areas in the presence of precipitation. It is shown that the presence of atmospheric precipitation can be taken into account with the help of introduction of effective diffusion coefficient, which in the general case is a function of wind speed.

Key words: eddy diffusion, stationary source, greenhouse gases, absorption.

Сведения об авторе

Чернявский Сергей Анатольевич, 1989 г.р., окончил МГТУ ГА (2011), соискатель кафедры безопасности полётов и жизнедеятельности МГТУ ГА, ведущий технолог центра информационных технологий и средств обучения, область научных интересов - математические аспекты инженерной экологии и безопасности полётов.

Рис. 2. Зависимость эффективного коэффициента