2
Моделирование процесса переноса и диффузии мелкодисперсных частиц в атмосфере с учетом эрозии почвы
Н. Равшанов1, д-р техн. наук, профессор, заведующий лабораторией
Д.А. Каршиев2, доцент, заведующий кафедрой
Б.Э. Юлдашев1,2, доцент
1Научно-инновационный центр информационно-коммуникационных технологий при Ташкентском университете информационных технологий
2Ташкентский педиатрический медицинский институт
(Узбекистан, г. Ташкен)
Аннотация
В работе рассматривается решение актуальной задачи по прогнозированию экологического состояния окружающей среды, на примере Приаральского региона, где за счет эрозии почвы осушенной части Аральского моря в атмосферу выбрасывается большое количество соли и мелкодисперсных вредных частиц. Приведены математическая модель, ее информационное обеспечение и результаты проведенных вычислительных экспериментов на ЭВМ. Разработанная модель, предназначена для исследования, мониторинга и прогнозирования процесса распространения загрязняющих веществ в атмосфере, и описывается системой дифференциальных уравнений в частных производных с соответствующими начальными и краевыми условиями. Указаны основные параметры, которые играют существенную роль в процессе переноса и диффузии вредных веществ в атмосфере: скорость ветра и его направления; рельеф местности; объем выброса солевых частиц в атмосферу, который зависит от скорости ветра и влажности почвы и т.д. На основе анализа проведенных численных экспериментов сформулированы выводы, служащие основанием для принятия управленческих решения по данной проблеме.
Ключевые слова: математическая модель, численный эксперимент, экология, эрозия почвы, вредные вещества, скорость ветра, мониторинг, прогнозирование.
осушенное море эрозия почва атмосфера
Введение
Растительный покров на поверхности земли, объекты производства, гражданская застройка, рельеф местности и прочие орографические элементы существенным образом влияют на изменение направления и скорости ветра, приводят к турбулентному перемещению ветра в пограничном слое атмосферы. В зависимости от указанных объектов, толщина действия турбулентного перемещения может достигать сотен метров. Зачастую турбулентное перемещение воздушной массы влечет эрозию почвы с выносом мелкодисперсных частиц в атмосферу. Из данного обстоятельства следует, что помимо стационарных и нестационарных источников вредных веществ в атмосфере, поверхность земли, при определенных условиях, также является генератором выброса вредных частиц.
Исследование рассеивания примесей в атмосфере стало одной из актуальных задач в проблеме охраны окружающей природной среды. Решение этой задачи связано с учетом многих факторов, влияющих на рассеивание примесей в атмосфере. К ним относятся метеорологические условия, тип источника, свойства примесей и т.д.
Значительный интерес представляет также изучение выноса в атмосферу и рассеивание легких и тяжёлых (имеющих собственную скорость осаждения) частиц естественного и промышленного происхождения. Существует множество источников загрязнения такого рода, в том числе, породные отвалы, промплощадки, осушенные части Аральского моря, солончаки Барса-кельмес, которые являются основным источником рассеивания примесей от наземных источников.
Проведенные исследования процесса миграции мелкодисперсных частиц от поверхности земли показали, что их количество, в основном, зависит от горизонтальной составляющей скорости ветра над поверхностью земли, состава и физико-механических свойств почвы, а также ее влажности [1]. Эта все дает дополнительную нагрузку на санитарное состояние атмосферы рассматриваемого региона. Решение данной задачи в контексте мониторинга и прогнозирования экологического состояния промышленных регионов является актуальной проблемой.
Нужно подчеркнуть, что развитие информационных технологий дало импульс для проведения комплексных исследований экологически значимых задач с привлечением конструктивных методов математического моделирования и вычислительного эксперимента на ЭВМ с целью принятия управленческих решений.
В настоящее время имеется огромное количество работ по моделированию процессов переноса и диффузии мелкодисперсных частиц.
Подробный анализ этих работ показал, что основными источниками выброса вредных веществ в окружающую среду авторами рассматриваются стационарные и нестационарные объекты: заводы, фабрики по производства сырья и продуктов питания, горноперерабатывающие предприятия, нефте- и газоперерабатывающие объекты, транспортные системы и т.д. Между тем, анализ состояния атмосферы мегаполисов и промышленных регионов показывает, что дополнительным источником загрязнения регионов является подстилающая поверхность земли, где со временем скапливается большое количество выпавших вредных частиц. Поверхность земли, покрытая осевшими вредными частицами, при определенных изменениях погодно-климатических факторов становится источником вредных веществ в атмосферу.
Примером подобного источника выброса вредных веществ в атмосферу для Средней Азии служит осушенная часть Аральского моря, где за счет эрозии почвы в атмосферу выбрасывается большое количество соли и других аэрозольные частиц.
Для решения задач мониторинга и прогнозирования экологического состояния регионов и принятия управленческих решений разработано математическое и программное обеспечение, в котором при моделировании объекта исследования учитываются погодно-климатические условия, эрозия почвы, физико-механические свойства аэрозольные частицы и другие факторы.
Постановка задачи
Для прогноза распространения аэрозолей в окружающей среде, нахождения количества аэрозолей в рассматриваемой области и количества выпавших аэрозолей на подстилающую поверхность, воспользуемся математической моделью [22, 24-27]:
(1)
(2)
(3)
(4)
(5)
(6)
Здесь - концентрация примеси; - мощность стационарного источника выброса вредных веществ в атмосферу; - составляющие скорости ветра; - коэффициент турбулентности; - скорость осаждения; - горизонтальный коэффициент турбулентности; - коэффициент поглощения; - коэффициент взаимодействия с подстилающей поверхностью. При имеем приподнятый источник на уровне . При наземных источниках . Задачу (1)-(6) рассмотрим в области когда источник расположен в приземном слое (рис. 1). На рис. 1 «крестиками» обозначены надземные источники загрязнения атмосферы.
Величина является функцией от и должна определяться по экспериментальным данным в зависимости от метеорологических условий, свойств подстилающей поверхности, размера и плотности пылевых частиц. В силу линейности дифференциального уравнения (1) мощности источников могут быть нормированы таким образом, чтобы максимальному значению соответствовала единица. Тогда, умножая на нормирующий множитель, получим искомые значения концентраций примесей. Соответственно коэффициент турбулентности также зависит от метеорологических условий.
Рассмотрим для определения следующие модели:
1.
2.
3.
где - высота приземного слоя, - турбулентная вязкость.
Рис. 1. Расположение источников на наземном слое в осушенной части Приаралья
Статическая обработка данных по влиянию различных погодно-климатические условий на унос частиц с поверхности земли показывает, что основным фактором, вызывающим эрозию почвы, является величина скорости набегающего потока, а главным фактором - влажность почвы. Остальные параметры хотя и играют определенную роль в развитии или предотвращении эрозии, но их влияние носит неоднозначных характер. Поэтому, при математическом моделировании процесса распространения вредных частиц и изменения их концентрации со временем, необходимо учитывать эти параметры, а физико-химические свойства почвы можно считать постоянными, входящими в расчетную формулу.
Скорость ветра и влажность почвы по-разному влияют на процесс распространения вредных частиц в атмосфере и развитие эрозии почвы. Если увеличение скорости набегающего потока ветра усиливает процесс эрозии, то повышение влажности почвы подавляет его.
Как было рассмотрено в [1] общую зависимость можно записать в виде:
(7)
Здесь - объемный расход уносимых частиц атмосферным фронтом, м3/с.
Для определения вида функции (7) переходим к анализу действующих сил, обуславливающих разрушение почвы и сопротивляющихся этому разрушению. Разрушающие силы обозначим через . Им всегда противостоят силы сопротивления , в состав которых входят влажность и другие физико-механические свойства почвы.
При превышении силы над силой начинается процесс эрозии почвы и уноса вредных частиц с поверхности земли. Сила разрушения в основном определяется величиной касательного напряжения набегающего потока воздуха. При этом, чем больше твердых частиц в потоке, тем больше суммарное касательное напряжение разрушающее почву. Для получения теоретической зависимости будем рассматривать равновесный процесс. При динамическом равновесии разность сил и должна равняться нулю, т.е.
. (8)
Составим выражение для этих сил. Связь между объемным расходом уносимых частиц и скоростью потока можно выразить в виде:
(9)
где -- касательное напряжение, кГ/м2.
Для силы сопротивления , по аналогии с возьмем выражение
(10)
где - вязкость смеси (воздух + почва), кГ*с/м2; - расстояние между отдельными частицами, м; -- константа почвы.
Подставляя (9) и (10) в (8) имеем
(11)
или
(12)
Рассмотрим отдельно взятое выражение в уравнении (12).
Если приближенно считать, что касательное напряжение определяется значением внешней скорости , т.е. , то это выражение в уравнение (12) можно представить в виде
(13)
В полученном выражении переменной величиной остается динамическая вязкость смеси , которая, кроме физико-механических свойств почвы, в основном, определяется влажностью почвы, обеспечивающей увеличение сил сцепления между отдельными частицами. На основании этого выражения (13) можем полагать зависящим от некоторой функции влажности, т.е.
. (14)
В дальнейшем функцию заменим простой зависимостью
(15)
Окончательно зависимость (14) примем в виде
(16)
где - константа почвы, с/м. Возвращаясь к выражению (12), будем иметь
(17)
Итак, получено уравнение для определения объема уносимых частиц с поверхности земли в зависимости от скорости ветра и влажности почвы. По результатам вычисленного значения - объемного расхода уносимых атмосферным фронтом частиц, можно решать задачу по процессу переноса и диффузии вредных веществ в атмосфере с помощью граничного условия (5).
Метод решения задачи
Таким образом, для рассматриваемой задачи имеем дифференциальное уравнение в частных производных первого порядка для определения объемного расхода уносимых атмосферным фронтом частиц. Интегрирование этого уравнения дает две произвольные постоянные С учетом константы почвы «», мы должны определить три постоянные Поэтому граничные условия должны быть заданы для трех значений . Чтобы использовать меньше аргументов, принимаем следующие граничные условия:
(18)
Решение (17) получено в виде
(19)
Причем, - отрицательная величина. Постоянные имеют следующие размерности
Формула (19) справедлива для случая Согласно граничным условиям получены следующие выражения для постоянных:
(20)
(21)
(22)
Для сравнения теоретически полученной зависимости (19) с экспериментами были вычислены коэффициенты для конкретного типа почвы, по данным К. Мирзажанова [3].
Например, для лугово-аллювиальной почвы, было получено:
Вычисленные значения уноса отличаются от экспериментальных в пределах 10-15%, что объясняется погрешностями эксперимента. Таким образом, на основании формул (20)-(22) по данным экспериментов можно вычислить значения для различных типов почвы. Так как опыты К. Мирзажанова [3] были проведены в области ламинарного пограничного слоя аэродинамической трубы (в начальном участке), то вычисленные коэффициенты при подстановке в формулу (19) дают значения только для уноса в пределах ламинарного пограничного слоя. Для использования формулы (19) в области турбулентного пограничного слоя следует проводить опыты в аэродинамической трубе в области турбулентного пограничного слоя и отсюда определять константы . Для этого каретка с образцом почвы должна находится на расстоянии не менее 0,8-1,0 м. от начала рабочего участка Изменение интенсивности уноса частиц в зависимости от скорости потока при различной степени увлажнения почвы .
Итак, на основе полученного решения задачи (17)-(18) можно вычислить количество выброса мелкодисперсных частиц в атмосферу в зависимости от составляющей скорости ветра на поверхности земли, влажности почвы, коэффициента шероховатости и других физико-механических параметров почвы (рис. 2-7).
Обсуждение результатов
Численные расчеты на ЭВМ были проведены при различных значениях влажности почвы (рис. 2-3). Как видно из рис. 2 и 3 - с ростом влажности почвы экспоненциально уменьшается объем выноса мелкодисперсных частиц в атмосферу. При стремлении влажности почвы к единице, уносимый объем частиц в атмосферу приравнивается нулю.