LINEBURG


страница 1
(всего 5)

ОГЛАВЛЕНИЕ

След. стр. >>

РОССИЙСКИЙ ФЕДЕРАЛЬНЫЙ ЯДЕРНЫЙ ЦЕНТР -
ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ
ТЕХНИЧЕСКОЙ ФИЗИКИ ИМ. АКАД. Е.И.ЗАБАБАХИНА

Регистрационный №
Инвентарный №

УТВЕРЖДАЮ
Зам.директора РФЯЦ-ВНИИТФ
Б.К.Водолага




Итоговый отчет по проекту МНТЦ KR-715
(за период с 1.12.2002 по 28.02.2006)




Название проекта:

'Исследования изменений основных экологических показателей территорий, примыкающих к урановому производству в условиях его частичной конверсии'



Головная организация:
Институт физики Национальной академии наук Кыргызской Республики (ИФ НАН КР)

Руководитель проекта:
И.А.Васильев - кандидат ф.-м. наук

Организации-участники:
1. Институт физики Национальной академии наук Кыргызской Республики (ИФ НАН КР)
2. Российский Федеральный ядерный центр - Всероссийский научно-исследовательский институт технической физики (РФЯЦ - ВНИИТФ)





Руководитель субпроекта в РФЯЦ - ВНИИТФ:
А.А. Шестаков

Длительность проекта: с 1.12.2002 по 28.02.2006.

Снежинск, 2006 г.



Содержание
Введение 3
I. Математическое моделирование физических процессов ветровой эрозии с поверхности хвостохранилища 7
I.1 Перенос пыли 7
1. Модель ветровой эрозии 7
2. Математическая программа ВЕТЭРО для моделирования ветровой эрозии 18
3. Численные расчеты 19
Литература к разделу I.1. 26
I.2 Перенос радона 27
1. Модель переноса радона 28
2. Математическая программа МИФ для описания распространения радона 31
3. Численные расчеты 33
3.1 Численное исследование одномерных задач на установление стационарного распределения радона над площадью хвостохранилища 33
3.2 Численное исследование задачи на установление стационарного распределения радона над площадью хвостохранилища в безветренную погоду 40
3.3 Численное исследование распределения радона над хвостохранилищем при постоянном ветре 44
3.4 Численное исследование ширины перемешивания радона и воздуха при различных скоростях ветра 48
Литература к разделу I.2. 51
II. Создание геомиграционной модели процесса переноса загрязнений в подземных водах с учетом сорбционной модели миграции урана в многокомпонентных растворах через основание хвостохранилища 52
1. Оценка формирования загрязнений подземных вод 53
2. Модель проникновения урана в подземные воды ХВХ 57
3. Защитные свойства подстилающих пород хвостохранилища КГРК 58
4. Выбор индикатора 63
5. Обоснование выбора математической модели миграции урана в подземные воды ХВХ 63
6. Математическая программа УРАН для расчета концентрации урана в ХВХ 67
7. Численные расчеты 68
Литература к разделу II. 75
III. Создание Web-узла в сети ИНТЕРНЕТ с информацией об экологической ситуации в районах Кара-Балтинского горнорудного комбината 77
1. Создание сайта Kara-Balta в сети ИНТЕРНЕТ 77
2. Краткий обзор сайтов в ИНТЕРНЕТ по радиационной экологии 81
Заключение 84


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

ВВЕДЕНИЕ

Целью настоящего проекта, наряду с разработкой некоторых аспектов методологии реабилитационных исследований, является подготовка пакета рекомендаций по использованию Кара-Балтинского горнорудного комбината (КГРК) для создания предприятий, не связанных с военно-промышленным комплексом. Выполнение проекта будет способствовать решению крупной проблемы - реабилитации территорий Кыргызской Республики, на которых были расположены предприятия по добыче и переработке урановых руд; подготовке и созданию рабочих мест на реабилитированных территориях и вновь образованных предприятиях.
Ранее город, в котором находится КГРК, назывался Кош-Тегирмен, позднее (9 сентября 1975 г), когда его объединили с селом Калининское, он получил название - город Кара-Балта. Он расположен на высоте 650-750 метров над уровнем моря, в 62 км на запад от Бишкека. В 1993 году в городе проживало 54,2 тысячи человек. Здесь живут в основном кыргызы, русские, украинцы, казахи и др. В городе имеется сахарный завод, 2 хлебных, консервный и молочные заводы, пивзавод, завод электротехнической игрушки, быткомбинат, 6 строительных и 4 транспортных организации, ковровый комбинат и др. Также имеется в городе 11 средних школ, медучилище, 2 профессионально-технических училища, пять клубов.


Фото 1. На переднем плане дамба Карабалтинского хвостохранилища. Здесь в 1959 году после многодневного ливневого дождя произошел прорыв дамбы и селевидная радиоактивная масса залила на много километров окружающие ее поля.

Гидрометаллургический завод (ГМЗ) Кара-Балтинского горнорудного комбината был создан в 1950 г. как комплекс горнорудодобывающих и перерабатывающих предприятий атомной отрасли бывшего Советского Союза. Он функционировал на базе Каджи-Сайского и Кавакского урановых месторождений. К 1955 году, с завершением строительства в г. Кара-Балта гидрометаллургического завода, ТЭЦ, вспомогательных производств, ГМЗ становится крупным предприятием отрасли по выпуску уранового концентрата. Был освоен выпуск попутно извлекаемого из руд молибдена в виде парамолибдата аммония. Наряду с выполнением основного своего предназначения завод освоил переработку молибденсодержащих отходов. Весь выпуск молибдена обеспечивался за счет переработки отходов и промпродуктов, которые ранее поступали более чем с 20 предприятий СССР. Наряду с традиционно выпускаемым аммонием молибденокислым также интенсивно осваивалась технология выпуска других продуктов. На заводе было создано производство по извлечению рения из промпродуктов других отраслей промышленности.
В 1989 году в Кыргызстане была прекращена добыча урановой руды подземным и открытым способами. Основным направлением на горнорудном комбинате стала переработка олововольфрамового концентрата и золотосодержащих руд и концентратов. Был создан цех по выпуску золота высокой чистоты из продукции комбината "Макмалзолото", а в последнее время там же проводится аффинаж Кумторского золота. Подразделения вспомогательного производства, кроме того, ранее выпускали соединения редких и радиоактивных металлов, комбайны для проходки восстающих горных выработок, буровые станки, резинотехнические, железобетонные и столярные изделия, респираторы, электродвигатели, антифрикционную молибденовую композицию "Моликом" для обогащения моторных и трансмиссионных масел, товары народного потребления. Если основным поставщиком урановой руды или полуфабрикатов являлись в 50 годы "Восьмерка" и Кавак, то позднее таковым становится Казахстан. С распадом СССР на Карабалтинском горнорудном комбинате резко сократилось производство урана. В 1997 году за год было произведено урана не более 500 тонн. Потребность в уране в таких странах, как США, Германия, Франция и Япония, достаточно велика. Россия также нуждается в центральноазиатском стратегическом сырье, так как после образования СНГ Россия фактически осталась без запасов урана. А он, в основном, необходим сейчас для деятельности АЭС, поэтому предполагалось, что в предприятиях, желающих купить продукт нового АО, недостатка не будет. На базе ГМЗ, совместно с Казахстаном создано АОЗТ "Уран". При этом хозяином акций стала казахская атомная компания Казатомпром, что позволяет ей диктовать условия. Уран на мировой рынок приходит под маркой Казахстана, а производство переносится на территорию Кыргызстана. Согласно существующему плану предполагается переработать 1500 тонн уранового концентрата за три планируемых года. Около половины от этого объема останется на территории Кыргызстана в виде вредных отходов.


Фото 2. Здесь, в заболоченной впадине, в 1950 году в результате аварийного сброса радиоактивных отходов с Карабалтинского гидрометаллургического завода образовалась техногенная радиоактивная аномалия (до 500 мкр/час).

Рядом с ГМЗ находится хвостохранилище, а от г. Кара-Балта - в 1,5 км. Там ведется регулярный контроль за состоянием подземных вод. В последние годы наблюдается тенденция снижения, как площади ореола загрязнения, так и концентраций загрязняющих веществ (это в основном сульфаты, нитраты). У экологов есть опасения с точки зрения безопасности хранилища, так как в процессе эксплуатации не всегда соблюдались проектные решения и хотя построена дополнительная заградительная дамба, хранилище не обезвожено, т.к. не действует система полного водооборота. В 1959 году на Кара-Балтинском хвостохранилище после многодневных ливневых дождей произошла авария - размыло дамбу. Радиоактивные отходы затопили поля, попали в ирригационную сеть. В том месте, где прошли радиоактивные отходы, выросла гораздо более высокая пшеница и дольше не осыпалась. Оказалось, что в радиоактивных отходах было много азота, и он то повлиял на рост пшеницы. Пшеница не содержала в себе никаких радиоактивных элементов, но урожай ее весь сожгли, так как при уборке за счет пыли она получала радиоактивное загрязнение. После этой аварии в целях исключения из обихода у населения колодезной воды в селе Калининское был построен водопровод. Кроме того, в Кара-Балте была выявлена радиоактивная аномалия в результате аэрогаммасъемки с вертолета (в 1990 году). Она представляла собой заболоченную впадину, заполненную радиоактивными отходами ГМЗ. Радиоактивная аномалия была сформирована еще в 50 годы в начальный период работы гидрометаллургического завода в результате аварийного сброса по канализационным трубам радиоактивных отходов. После ее выявления радиоактивный слой был снят и вывезен на карты хвостохранилища ГМЗ. Уровень гамма-радиации здесь ранее превышал 500 мкр/час. Сейчас здесь отмечается повышенная радиоактивность, но не более 60 мкр/час. Взятая здесь болотная вода показала суммарную альфа-активность в 3,8*10-10 Ки/л, суммарную бета-активность в 5,5*10-11 Ки/л (ПДК питьевой воды составляет по альфа и бета активности - 3*10-11 Ки/л).
Хвостохранилище ГМЗ действует более 50 лет. Оно намывное, равнинного типа, площадью порядка 240 га, состоит из отдельных карт. Карты недостаточно экранированы, ореол загрязнения хвостохранилища имеет грушевидную форму, вытянутую по направлению потока с севера на юг, на расстоянии около 14 км от заградительной дамбы. Площадь загрязнения 40-50 км2. Глубина загрязнения не превышает 110-120 м. Промстоки не оказывают влияния на качество подземных вод глубже 150 м. Несмотря на весь комплекс природоохранных мероприятий, проводимых ГМЗ, загрязнение подземных вод остается стабильным по площади и во времени. Максимальное содержание загрязняющих компонентов отмечается в районе заградительной дамбы хвостохранилища, в скважинах водоперехвата.

Фото 3. Аккумулирующие емкости (на заднем плане - населенные пункты, находящиеся вблизи хвостохранилища).

В свое время все руководители горнорудной и перерабатывающей урановой промышленности, начиная от министра Средмаша Славского Е.П., и, особенно начальника Первого Главка Средмаша Карпова Н.Б., нацеливали руководителей урановых предприятий на местах на то, чтобы все урановые подразделения в отношении санитарии содержались на порядок выше, чем в химической промышленности. Но хвостохранилища выщелаченных руд до сих пор представляют из себя радиоактивную опасность. Эти отвалы истощены в смысле наличия природного урана, но в них имеются другие радиоактивные элементы. Они являются источником эманации радона и выщелачивания тяжелых металлов, кислотных растворов и токсичных элементов. Кроме того, в хвостохранилищах размещены частично непереработанные руды - это отвалы бедных руд. Хотя опасность от эманации радона и выщелачивания радиоактивных элементов относительно низка, но при использовании этих руд в качестве строительного материала для дорожного покрытия или домов, внешняя радиации может быть высока. Отвалы должны быть закрыты и просачивающиеся воды должны контролироваться и обеззараживаться в отстойниках, которые также представляют определенную опасность. Они защищены дамбами, которые могут прорваться и просачивающиеся воды могут загрязнить как поверхностные, так и подземные воды. Обедненные руды подземного выщелачивания, обработанные серной кислотой или щелочным раствором, остаются основным источником загрязнения грунтовых вод, даже если откачка химических растворов закончена. Наиболее существенным являются загрязнения поверхностных и грунтовых вод, используемых в качестве источника питьевой воды, а также для сельскохозяйственных и строительных целей.
Исследования изменений основных экологических показателей территорий вокруг хранилища радиоактивных отходов в районе г. Кара-Балта можно разделить на процессы переноса загрязнений в атмосфере и процессы переноса загрязнений в подземных водах (рис.1). Процессы переноса в атмосфере в свою очередь можно разделить на процессы ветровой эрозии с поверхности хвостохранилища (перенос пыли) и процессы переноса радиоактивного газа радона (в проекте это описано в задаче 12). Процессы переноса в подземных водах можно разделить на процессы переноса радиоактивных и химических загрязнений в зоне аэрации и процессы переноса радиоактивных и химических загрязнений в зоне насыщения (в проекте это описано в задаче 13). Загрязнение еще может переноситься водным путем через реки, ручьи, родники, сели, но в данном проекте эти пути не рассматриваются, т.к. главная ожидаемая опасность распространения загрязнений от Карабалтинского хвостохранилища состоит в инфильтрации радионуклидов в подземные воды и пылеуносе радионуклидов на жилые и сельскохозяйственные площади.

Геомиграционная модель процесса распространения загрязнений вокруг хранилища



Процессы переноса в атмосфере Процессы переноса Процессы переноса в подземных водах
в наземных водах




перенос пыли перенос радона в зоне аэрации в зоне насыщения


Рис.1. Схема распространения загрязнений вокруг Карабалтинского хвостохранилища.




I. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ ПРОЦЕССОВ ВЕТРОВОЙ ЭРОЗИИ С ПОВЕРХНОСТИ ХВОСТОХРАНИЛИЩА

I.1 Перенос пыли

1. Модель ветровой эрозии
Изучение явления ветровой эрозии с поверхности ХВХ является актуальным, т.к. ветровая эрозия - второй после водной фильтрации фактор распространения загрязнений. Под ветровой эрозией (дефляцией) понимается совокупность взаимосвязанных процессов отрыва, переноса и отложения частиц ветром. Но большая или меньшая подверженность почв ветровой эрозии определяется не одним лишь ветром, а также совокупностью таких факторов как режим атмосферных осадков, температура, рельеф, гранулометрический, агрегатный состав почвы и другие. Таким образом, хорошо проанализировав все эти факторы можно выделить места наиболее подверженные дефляции. Работы Р.А. Багнольда, осуществленные в 30-ых годах и опубликованные в его выдающейся монографии "Физика выдуваемого песка и дюны пустынь", ознаменовали собой кардинальный сдвиг в понимании системы ветровой эрозии. Багнольдом была рассмотрена проблема зависимости "ветер/песок" как один из видов аэродинамики, поддающийся непосредственному измерению. Он апробировал свои гипотезы в лабораторных экспериментах с использованием аэродинамической трубы, а также в полевых испытаниях в Ливийской пустыне. Используя аэродинамическую трубу, Багнольд принимал, что крупные завихрения не играют значительной роли в поддержании песка в воздухе и что перемещение песка происходит примерно в пределах метра от поверхности земли. Многие идеи и результаты, полученные Багнольдом, создали основу для дальнейших исследований проблемы ветровой эрозии.
Сегодня существует много различных моделей ветровой эрозии. Но все они сходятся в одном, необходимо как можно более точно решить проблему, связанную с прогнозированием ветровой эрозии и разработать надежные и экономически выгодные методы защиты от эрозии. Для решения этих вопросов необходимо глубокое понимание этого на первый взгляд простого процесса. Изучение ветровой эрозии осложняется разномасштабностью составляющих ее процессов, в основе которых лежат разные механизмы. Ученые разных стран по-разному пытаются решить эту проблему, стремясь отыскать наиболее правильный подход. Разрабатываются все новые и новые модели ветровой эрозии, которые имеют дальнейшее применение на практике в противоэрозионных мероприятиях. И от того, насколько модель точна и проста в расчетах, можно судить об ее пригодности к применению на практике.
При описании основных физических процессов ветровой эрозии существует необходимость в проведении сравнительного анализа различных моделей, для того, чтобы выявить наиболее полное и комплексное решение этой задачи. Это обозначит правильные направления в дальнейших исследованиях основных закономерностей ветровой эрозии и позволит использовать полученную модель в практических целях на Карабалтинском ХВХ. При проведении работы по предварительному описанию модели ветровой эрозии с поверхности ХВХ был проведен анализ двух моделей. Первая модель принадлежит американским ученым Роберту С. Андерсону и Бернарду Халлету. Их работа называется "Общая модель переноса частиц почвы ветром" и является обобщением общепринятых представлений и подходов в ветровой эрозии [1]. Вторая модель разработана российскими учеными из МГУ Г.П. Глазуновым и В.М. Гендуговым и изложена в работе "Механизмы ветровой эрозии" [2]. Она является достижением в том плане, что в ее разработке успешно используются методы, ранее не применявшиеся и, в итоге, предлагается физически содержательная модель. Обе модели являются фундаментальными научными работами известных ученых в области эрозиоведения и поэтому заслуживают пристального внимания.
Эти две модели, имея некоторые общие выводы, различаются кардинально в фундаментальных понятиях и подходах к изучению механизмов ветровой эрозии почв. Так, например, российские исследователи давно используют представление о вихревой природе сил, отрывающих частицу от поверхности. Исходя из этого представления, строятся все дальнейшие рассуждения. Американские же ученые, принимая во внимание этот факт, в основном пренебрегают им при построении своей теории.
Работа Г.П. Глазунова и В.М. Гендугова посвящена выявлению закономерностей выдувания на микро- и макроуровнях, исследованию структуры почво-воздушного потока и обобщению полученных результатов в виде физически содержательной модели потерь почвы от ветровой эрозии. В данной работе выводится уравнение выдувания почвы, выявляется физический смысл параметра массообмена в уравнении выдувания, решается задача о скорости вылета почвенной частицы с поверхности под действием ветра в поверхностном слое почвы, выводится теоретическое уравнение траектории почвенных частиц в воздухе, определяется понятие и способ нахождения критической скорости ветра, при которой начинается горизонтальный полет почвенной частицы, обсуждается и выясняется структура почво-воздушного потока, решается задача о возможных потерях почвы от ветровой эрозии с бесконечного поля.
Работа Андерсона и Халлета отличается от работы российских ученых изначально другим подходом к изучению механизмов ветровой эрозии. Американцы проводят четкую границу между процессами скачкообразного передвижения частиц и движением частиц в подвешенном состоянии. Соответственно, они и ищут отдельные решения для этих двух процессов, в отличие от работы российских ученых, в которой найдено общее решение для обоих процессов. Их работа посвящена выявлению закономерностей механизмов ветровой эрозии на микро- и макроуровнях, исследованию структуры почво-воздушного потока, исследованию начальных условий, при которых происходит отрыв частицы от поверхности. Обсуждается природа сил действующих на частицу, в результате чего приводятся уравнения зависимости концентрации и массового потока от определенной высоты для частиц передвигающихся скачкообразно. Для этого же процесса авторами приводятся различные эмпирические закономерности, требующие проверки. Также предлагается уравнение для плотности вероятности скорости вылета, распределение которой имеет прямую зависимость с распределением длины скачка. Для отдельно рассматриваемого процесса переноса частиц в подвешенном состоянии дается вывод уравнения для общей концентрации частиц в определенном объеме.
Российские авторы при решении проблемы математического описания процесса ветровой эрозии использовали метод перехода от микроуровня, при котором анализируется движение отдельной частицы, к макроуровню, при котором процесс описывается в рамках законов механики многофазных сред. Авторы исходили из представления о том, что существует некоторая пороговая скорость ветра Uкр. (м/с), называемая критической, при превышении которой начинается интенсивное выдувание почвы, характеризуемое величиной q (кг/м2/с). Интенсивность выдувания при данной скорости Ue на границе пограничного слоя зависит от касательного напряжения трения ? (Н/м2). Плотность энергии Е (Дж/кг), необходимая для выдувания почвы, может быть выражена через кинетическую энергию ветра при его значении, равном критическому для почвы: Е =0.5(Uкр.)2. Величины Ue, q, ?, E достаточно полно описывают процесс ветровой эрозии. Из них можно составить только две безразмерные комбинации: и , которые являются независимыми в том смысле, что любые другие безразмерные комбинации могут быть выражены через них. Величина В представляет собой параметр массообмена, аналогичный по физическому смыслу параметру массообмена в физической теории испарения, привлекаемому для описания газодинамических величин на границе раздела пограничного турбулентного слоя и граничной поверхности.
Согласно ?-теореме (Седов, 1972), физически процесс ветровой эрозии почв может быть определен с помощью функции, связывающей эти параметры: [3]. Вид этой зависимости найден российскими авторами экспериментально и потом подтвержден математическим выводом. Строился график, на оси ординат, которого откладывались значения ln B, а на оси абсцисс . Во всех случаях на графиках ln наблюдается одна и та же зависимость: при значениях скорости, превышающих Uкр., указанная зависимость имеет вид прямой линии. При скоростях, не превышающих Uкр., закономерная связь между переменными не прослеживается. Следовательно, при скоростях больших, чем Uкр., искомая зависимость имеет вид: . Здесь ? - эмпирический коэффициент, характеризующий свойства почвы, которые определяют ее устойчивость к выдуванию. Его находят, измеряя тангенс угла наклона прямой. Полученные из опытов с монофракциями значения этого коэффициента свидетельствуют о большей "скорости" нарастания выдувания почвы с увеличением скорости потока. Для реальных почв следует ожидать меньших по абсолютной величине значений ?. Величины коэффициента ? имеют тенденцию к увеличению с увеличением размера частиц. Коэффициент ? = ln Bкр. находят по пересечению прямой с осью у. Точка пересечения и будет ? = ln Вкр.
Эта зависимость выполняется при скоростях потока, превышающих Uкр., поэтому ее преобразовывают к виду, содержащему параметр Вкр., отражающий это обстоятельство. В итоге получается уравнение выдувания почвы: . Уравнение выдувания описывает поток почвенных частиц, направленный от поверхности в атмосферу. Это уравнение позволяет утверждать, при превышении критической скорости начинается закономерное выдувание почв, которое можно прогнозировать на основе выведенного уравнения. Для того, чтобы данное уравнение можно было применять для почв с реальным распределением частиц по размерам, необходимо экспериментальное определение соответствующих значений коэффициента ? и параметра Вкр.
Традиционный подход к моделированию сил, действующих на частицу, использовали и российские авторы, и американские. Он основан на рассмотрении баланса сил, приложенных к единичной частице поверхностного слоя почвы, и последующей экстраполяции результатов на все остальные частицы.
В обеих работах показателем поведения частицы при действии данных сил будет служить траектория частицы, о которой будет сказано позже.
В русской теории принимается существенным действие только потока на частицу, влиянием же частиц на поток пренебрегают, т.к. масса частиц потока определенного объема составляет ничтожную часть в масштабах ветровой эрозии от массы воздуха рассматриваемого объема. В американской теории говорится о важности действия частиц на поток, но все равно это нигде не учитывается в силу сложности подобных расчетов.
Для рассмотрения сил, действующих на частицу, как и принято, вводится система координат так, что ее начало совпадает с местом вылета частицы. Ось х совпадает с направлением ветра, а ось у - с нормалью к почвенной поверхности, направленной в сторону атмосферы. По мнению российских авторов, в направлении оси х на частицу радиуса ri, действует только сила лобового сопротивления , где К* - коэффициент лобового сопротивления, ?в - плотность воздуха, ui - составляющая скорости частицы в направлении оси х, U - скорость ветра. Эта сила придает частице ускорение в направлении этой оси. Поэтому согласно второму закону Ньютона можно записать:
,
где m - масса почвенной частицы плотностью ?п, равная (4/3)?ri2?п, а t - время.
В направлении оси у на частицу действуют? подъемная сила Жуковского , направленная вертикально вверх, сила Архимеда , результирующая которой направлена вниз, и сила сопротивления Стокса , где - коэффициент подъемной силы, g - ускорение свободного падения, ? - коэффициент формы частицы в законе Стокса, ? - вязкость воздуха, vi - составляющая скорости движения частицы в проекции на ось у. Уравнение движения почвенной частицы в проекции на ось у имеет вид:
.
Конечно же, на частицу действуют и другие силы: упругие силы, возникающие при взаимном соударении частиц или при ударе их о подстилающую поверхность, электрические силы, сила Магнуса. Считается, что силы электрической природы не вносят существенного вклада в перемещение почвенных частиц ветром. А сила Магнуса, возникающая при вращательном движении частиц в потоке, оказалась недостаточной для подъема частиц вследствие сравнительно малой скорости их вращения. Поэтому российские ученые пренебрегают ею.
Российские ученые считают, что главную роль играют подъемная сила и сила лобового давления. Подъемную силу связывают с возникновением различия в давлениях между верхней и нижней поверхностями обтекаемой потоком частицы. Подъемная сила, определяемая таким образом, быстро убывает по мере удаления почвенной частицы от поверхности и на высоте в несколько диаметров частицы она стремится к нулю, что не позволяет объяснить подъем частиц на большие высоты. Прямое измерение сил сцепления между частицами затруднено, поэтому ее значение совместно со значением коэффициента подъемной силы ищут по уравнению , где FС является проекцией суммы сил сцепления на ось ординат. Это уравнение представляет собой проекцию сил, действующих на частицу в момент, предшествующий отрыву. В момент, предшествующий взлету частицы вертикальная составляющая скорости равна нулю, поэтому уравнение и имеет такой вид. С подстановкой оно будет выглядеть следующим образом: = mg +FC . Неизвестными в этом уравнении являются FC и . Данное уравнение линейное и поэтому по тангенсу угла наклона находим , а по пересечению прямой с осью ординат находим точку, которая и будет FC .
Сильно отличаются представления американских ученых о наиболее значимых силах, которые следует учитывать при расчетах траектории частицы. Рассматриваемые ими силы: сила тяжести (Fg), касательное напряжение (Fd), подъемная сила (Fl) и подъемная сила Магнуса (Fm), возникающая в результате вращения частицы.
Ввиду того, что плотность частиц несоизмеримо больше плотности воздуха, сила тяжести становится просто равной массе частицы. Касательное напряжение действует на частицу в том же направлении, в котором действует Urel = U - Up, где U - средняя скорость ветра, Up - скорость частицы в потоке. При подъеме вектор направлен под таким же углом, как и вторая половина траектории, где частица опускается. При падении картина обратная.
Для касательного напряжения , где А - сечение частицы, перпендикулярное вектору силы касательного напряжения (A= для частиц сферической формы с диаметром D), ?а - плотность воздуха, Сd - коэффициент касательного напряжения. Для сферических частиц Сd является функцией мгновенного значения числа Рейнольдса, , где v кинематическая вязкость воздуха. Андерсон и Халлет, используя данные Уайта и Шульца принимали, что Сd = Сd(Re).
Подъемная аэродинамическая сила, по мнению авторов, возникает в результате образования разности давлений на противоположных концах частицы, достаточного для подъема ее в направлении возрастающей скорости. Уравнение, используемое авторами для подъемной силы, выглядит следующим образом: , где и - скорости наверху и у основания частицы, - коэффициент подъемной силы, значение которого авторы приняли 0.85 (Чепил, 1958).
Авторы полагают, что подъемная сила имеет особое значение только у самой поверхности, где создается максимальная разница в скоростях сверху и снизу частицы. Например, для частиц с радиусом 0.25 мм на высоте, равной 10 диаметрам частицы, значение подъемной силы уменьшается более, чем в 10 раз от первоначального значения. Подъемной силой можно пренебречь, считают авторы, в случае отрыва частиц вследствие удара другой частицы и передачи энергии первой. Но для частиц, которые не поднимаются высоко над поверхностью, эта сила оказывается очень важна. Подъемная сила вносит свой вклад в начало скачкообразного движения частиц, которое в общем-то требует совместного влияния касательного напряжения и подъемной силы.
Сила Магнуса рассматривается американскими учеными как основная подъемная сила. Уравнение для нее выглядит следующим образом: , где - угловая скорость частицы, значение которой принимают положительным для верхнего вращения. Ускорение вращения частицы происходит в результате разницы скоростей ветра на противоположных концах частицы. Однако, принимается во внимание и действие вязкости воздуха, замедляющего скорость вращения частицы.
Сила Магнуса будет давать частице ускорение при подъеме до тех пор, пока горизонтальная скорость частицы не станет равной скорости ветра. Если рассматривать горизонтальную составляющую силы Магнуса, то будет очевидно, считают американские ученые, что при взлете сила будет давать частице ускорение, а при падении, наоборот, тормозить.
Все эти силы можно разложить по осям х и у и составить суммарные по осям уравнения, из которых можно найти вертикальную, горизонтальную составляющие начальной скорости частицы, а также начальную угловую скорость. Но для этого необходимо знать угол подъема, скорость подъема частицы. Эти данные авторы получают, решая два уравнения с подстановкой известных данных о диаметре, плотности частицы и профиле скорости.
Следующий важный факт, который следует рассмотреть, это то, как описываются траектории частиц в американской и российской моделях. Американские ученые считают, и подтверждают свое мнение фотографиями, что отдельные траектории частиц имеют асимметричную форму, где угол вылета больше угла падения частицы. Российские ученые вполне согласны с ними, но уточняют, что поскольку, как говорилось выше, частица летит по хаотичной траектории, форма этой траектории для отдельной частицы может быть какой угодно различной. В российской модели рассматривается некая усредненная траектория для определенных сортов частиц, имеющая симметричную выпуклую форму.
Халлет и Андерсон не предлагают нам никакого конкретного уравнения для нахождения траектории частицы. Они предлагают методом итерации, т.е. небольшими шажками, рассчитывая для каждого этапа мгновенное значение числа Рейнольдса, коэффициент касательного напряжения, скорость вращения, искать траекторию отдельной частицы с заданными параметрами, где учитывается их мгновенное изменение. Ясно, что встает необходимость проведения громоздких расчетов, которые осуществляются машинами ЭВМ. Когда траектория частицы определена, то не составляет труда найти длину и высоту скачка.
Российские ученые дают в свой фундаментальной работе вывод уравнения траектории, основанный на балансе сил, приложенных к единичной частице. При выводе используется уравнение движения почвенной частицы . С помощью несложных математических преобразований этих выражений получается следующее уравнение? , где для расчетов необходимо иметь данные лишь о радиусе ri , скорости ветра и начальной скорости подъема частицы vi0. Значения также можно найти, используя уравнение vio = kU, выведенное авторами, где в первом приближении принимают . Параметры bi и ai - константы, и . Начальная скорость подъема принимается одинаковой для всех частиц данной почвы, что объясняется неразборчивостью вихрей, действующих на частицу.
Если скорость ветра U меньше скорости при которой начинается перенос частиц в подвешенном состоянии, то наступит момент, когда частица снова окажется на поверхности, т.е. у = 0, как и в момент взлета. Этот факт используется для нахождения длины скачка х методом подстановки у = 0 в уравнение траектории. Получаемое уравнение решается графически. Это возможно, поскольку его левая часть представляет собой уравнение прямой, а правая уравнение экспоненты. Искомая длина хс соответствует абсциссе точки пересечения графиков этих двух функций.
Из уравнения траектории можно найти абсциссу и ординату наивысшей точки траектории скачка, которая является важным показателем при обсуждении структуры почво-воздушного потока.
Если рассмотреть случай, когда скорость потока равна Uкр.2, т.е. начинается горизонтальный полет частиц, уh = Н, т.е. высота скачка достигает предельной величины, которая является по сути высотой горизонтального полета. Радиус частицы, совершающий при скорости Ue горизонтальный полет, назовем критическим радиусом, rкр. Величины H и rкр легко выводятся из уравнения траектории и в итоге получаются уравнения: H = , rкр =.
В американской статье делается несколько выводов о различных зависимостях траектории частицы, которые являются результатом обобщения эмпирических данных. Так, например, ученые говорят, что поскольку вертикальное торможение частицы, пропорционально квадрату вертикальной компоненты скорости отрыва частицы, то получается, что частица, вылетевшая под большим углом и с высокой скоростью, упадет быстрее, чем та, которая вылетит под меньшим углом к поверхности и с меньшей скоростью. Это иллюстрируется графиком зависимости высоты скачка от скорости отрыва частицы. По нему выходит, что чем больше скорость отрыва частицы, тем ниже подскакивает частица, и, соответственно, чем меньше скорость отрыва частицы, тем выше она подлетает.
Длина траектории также сильно зависит от вертикальной компоненты подъемной скорости. Эта зависимость иллюстрируется графиком, который говорит о том, что чем больше скорость отрыва частицы, тем больше длина скачка. Это никак не противоречит теории русских ученых. Такая зависимость отражена в уравнении траектории частицы, где начальная скорость подъема частицы прямо пропорциональна и высоте, и длине траектории.
Халлет и Андерсон считают, что поскольку торможение обратно пропорционально размеру частицы, то маленькие частицы не поднимаются высоко, но получают большее горизонтальное ускорение, в отличие от больших частиц. Это утверждение требует более детального рассмотрения.
Американские ученые утверждают также, что угол снижения частицы является величиной, не зависящей от размера частицы, а лишь немного зависит от скорости отрыва частицы. По данным российских ученых он также является величиной варьирующейся в очень узких пределах. Но этот вопрос еще требует соответствующего изучения.
С разных позиций подходят ученые к такому важному вопросу как структура почво-воздушного потока. Российские ученые подошли к изучению этого вопроса, используя методы механики сплошной среды. Они ввели понятия многоскоростного континуума и взаимопроникающего движения составляющих почво-воздушного потока [3,4]. Составляющими потока является множество m компонент. Воздух характеризуется средней скоростью, направленной вдоль абсциссы, а также постоянной плотностью ?в. Каждый i-тый почвенный компонент представляет собой совокупность всех передвигающихся в воздушном потоке почвенных частиц i-того сорта, заполняющие тот же объем, что и другие почвенные континуумы и воздух, и именно поэтому движение континуумов является взаимопроникающим. Каждая частица i-того компонента имеет характерный радиус ri (i = 1, ..., m), плотность ?п и скорость vi.
Механика смесей строится на основе физических законов сохранения массы, импульса, энергии, поэтому можно записать уравнение сохранения массы потока почвенных частиц i-того сорта и использовать его для вывода. Для этого выделяется элементарный объем пространства в виде кубика и в нем рассматривается изменение массы, которое равно разности входящих в элементарный объем и выходящих из него потоков частиц. Не сложным математическим путем в статье выводится доказательство того, что вдоль любой i-той траектории концентрация сi частиц i-того сорта остается постоянной. Этот факт позволяет анализировать структуру почво-воздушного потока при том условии, что имеются данные о составе и свойствах частиц почвы.
Анализ структуры почво-воздушного потока российские ученые начинают со случая выдувания идеальной почвы, состоящей из частиц одного сорта, ri, в условиях "бесконечного" поля. Под "бесконечным" понимается та часть поля, где можно пренебречь его границами. Наиболее важным свойством структурного потока является пространственное распределение концентрации почвенной фазы в потоке.
При скорости потока меньшей той критической скорости, когда начинается перенос частиц в подвешенном состоянии, но большей той, при которой начинается скачкообразное движение, частицы i-того сорта будут передвигаться по предположительно одинаковым выпуклым траекториям. Поскольку частицы предполагаются абсолютно одинаковыми, их траектории также будут одинаковыми. Если представить это на графике у(х), то получится изображение пересекающихся парабол, концы которых являются местом вылета и падения частицы. В этих точках концентрация почвенных частиц будет оставаться постоянной, равной сi, вследствие равенства количества взлетевших и упавших частиц. Концентрация почвенных частиц на максимальной высоте подъема также будет оставаться постоянной, вследствие равенства количества поступающих и покидающих ее частиц.
Каждая траектория состоит из участков взлета и падения частицы. Из геометрического равенства траекторий одинаковых частиц следует, что произвольную восходящую ветвь пересечет бесконечное множество нисходящих ветвей траекторий, а нисходящую ветвь - бесконечное множество восходящих ветвей траекторий. Отсюда следует, что в любой точке потока, через которую проходит отдельная траектория, за исключением точек ее начала, конца и максимума высоты подъема, концентрация почвенных частиц равна 2сi. Этот вывод справедлив и применим ко всему слою сальтации одной фракции.
Если скорость достигает Uкр.2, все траектории частиц i-того сорта будут стремиться к одной горизонтальной плоскости, расположенной на определенной высоте скачка.
Если число сортов скачущих частиц равно двум, то образуются два взаимопроникающих слоя сальтации. Форма траектории зависит от размера частицы, чем мельче (легче) частица, тем выше и дальше она подскакивает при данной скорости потока, тем, соответственно, больше толщина слоя сальтации. В каждой точке меньшего слоя сальтации концентрация равна, где и - концентрации частиц первого и второго сорта. На верхней границе меньшего слоя сальтации концентрация твердой фазы при суммировании будет равна . Суммарная концентрация частиц на поверхности почвы равна .
Если частицы обеих фракций при данной скорости ветра способны лететь, то концентрация почвенной фазы в каждой точке потока будет равна сумме концентраций этих фракций . Таким образом, концентрация летящих почвенных частиц в любой точке потока над бесконечным эродируемым полем равна сумме их концентраций, , которая не зависит от высоты.
Все вышесказанное говорит о том, что почво-воздушный поток над эродируемым ветром полем имеет слоистую структуру. Почво-воздушный поток, состоящий только из двух сортов скачущих частичек, формирует шесть параллельно расположенных слоев с различными концентрациями почвенной фазы, естественным образом уменьшающимися с высотой. При большем количестве сортов частиц увеличивается количество слоев, но качественно картина не меняется.
В отличие от полной картины структуры почво-воздушного потока, которую дают нам российские ученые в своей модели, Халлет и Андерсон приводят два уравнения для расчета концентрации? одно для расчетов при скачкообразном движении частиц, другое для случая переноса частиц в подвешенном состоянии. Первое уравнение выведено на основе анализа отдельной траектории частицы. Оно дает зависимость общей концентрации от высоты при известных вертикальной скорости, количестве частиц, улетающих с единицы поверхности за единицу времени, известной зависимости концентрации от высоты частиц определенного сорта и рассчитанной плотности вероятности скорости вылета частицы. Плотность вероятности скорости вылета частицы рассчитывается по уравнению, в котором необходимо знать динамическую скорость и вертикальную компоненту начальной скорости вылета.
Уравнение для нахождения полной концентрации твердой фазы при условии, что все частички летят, выводится на основе использования созданной Шмидтом диффузной теории. Существенным ее недостатком является отсутствие представления о силе, перемещающей твердые примеси в потоках жидкостей и газов. Она не объясняет, каким образом частицы, имеющие плотность гораздо большую, чем у воздуха, поднимаются высоко над поверхностью. Ученые выводят свое уравнение для общей концентрации твердой фазы летящих в горизонтальном направлении частиц сложением предварительно найденных концентраций отдельных слоев с частицами разного размера.
Американские ученые говорят о том, что структура почво-воздушного профиля состоит из слоев с различной концентрацией преобладающих в них частиц, и общая концентрация уменьшается с высотой по степенному закону.
В обеих разработанных моделях очевидна большая заинтересованность авторов в решении задачи о возможных потерях почвы. Выведенное американскими учеными уравнение для массового потока как функции от высоты предлагается использовать как показатель потерь почвы. Уравнение это выведено также на основе анализа отдельных траекторий частиц. Для решения этого уравнения необходимо знать вертикальную составляющую начальной скорости вылета частицы, плотность вероятности скорости вылета, количество частиц, улетающих с единицы поверхности за единицу времени и зависимость массового потока от высоты для частиц одного размера, т.е. частиц, имеющих одинаковую траекторию.
Российские ученые в своей физически-содержательной модели предложили другое решение задачи о возможных потерях почвы. Эта задача идет заключительным этапом и является еще не полностью решенной, по мнению, авторов. Во-первых, для решения этой задачи рассматривается ограниченный участок бесконечного поля, где наблюдается ветровая эрозия при постоянной скорости ветра. Очевидна необходимость дальнейшего детального исследования влияния с последующим учетом таких факторов как: размеры эродируемого поля, рельеф, растительность, пожнивные остатки, абразия, динамика агрегатного состава и другие. Во-вторых, для вывода нужной нам формулы необходимо рассматривать закон сохранения массы потока почвенных частиц на поверхности участка бесконечного поля. Если Ue > Uкр., то формируются два потока частиц, один из которых движется от почвенной поверхности в атмосферу, а другой - в обратном направлении. Если для всех частиц, вылетающих с поверхности, выполняется условие Ue > Uкр.2, то поток из атмосферы в сторону почвы не формируется, поскольку все вылетающие с поверхности частицы уносятся ветром "безвозвратно". Однако такой случай маловероятен, и наряду с безвозвратно уносимыми частицами всегда имеются и перемещаемые скачками, поскольку распределение частиц по размерам в составе агрегатов, слагающих поверхность почвы, непрерывно. В таком случае, поток почвенных частиц, покидающих поверхность, состоит из двух потоков: потока частиц, уносимых безвозвратно, и потока частиц, которые переносятся скачками. Зная концентрации частиц в этих двух потоках, можно рассчитать потери почвы от эрозии. Общее уравнение имеет вид:
,
где Q - потеря почвы от ветровой эрозии (кг) за время T (с) с площади S (м2) при скорости ветра за пределами слоя шероховатости Uе (м/с), p1 и p2 - концентрации почвенных частиц (кг/кг), ? - касательное напряжение на почвенной поверхности (Н/м2), возникающее вследствие воздействия на нее ветра, имеющего скорость Uе, Bк - параметр массообмена (кг/кг), характеризующий данную почву при Ue = Uкр. , где Uкр. - критическая скорость ветра для данной почвы (м/с), по физическому смыслу аналогичная размывающей скорости водного потока, взятая на той же высоте над поверхностью, что и Uе; ? - безразмерная почвенная константа.
В уравнении потерь почвы от ветровой эрозии используется основополагающее в данной модели уравнение выдувания почв:
B = Bk,
выведенное Г.П. Глазуновым и В.М. Гендуговым, где В - параметр массообмена. Следует заметить, что не существует аналогов применения подобного параметра ни в одной из областей науки почвоведения. Огромное преимущество его заключается в том, что он, как величина безразмерная, является физически-содержательным параметром, т.е. при подстановке размерностей они сокращаются. В российской модели параметр массообмена В используется для вывода уравнения выдувания почв, являющегося физически-содержательным. Его физический смысл стал понятен после анализа этого процесса выдувания. Параметр массообмена представляет собой концентрацию частиц на поверхности почвы, потерявших межагрегатное сцепление под действием воздушного потока. Фактически, это частицы, принадлежащие не поверхностному слою почвы, а почво-воздушному потоку. Параметр массообмена можно найти имея данные о критической скорости ветра для данной почвы - Uкр. и о безразмерной почвенной константе - ?. Уравнение выдувания почвы имеет основополагающее значение для целей прогнозирования возможных потерь почвы от ветровой эрозии. Впервые это уравнение для почв в водных потоках было получено методами теории подобия и анализа размерностей Кузнецовым и Гендуговым.
Расходятся мнения авторов в представлениях о природе сил, отрывающих частицу от поверхности. Основанием для применения российскими учеными к почвенным частицам, являющимися плохо обтекаемыми телами, выражения для подъемной силы Жуковского послужило представление о вихревой природе силы, отрывающей частичку от поверхности и переносящей их по воздуху. Частица отрывается от поверхности замыкающимся на нее вихрем, а в процессе переноса передается от одного вихря к другому. Важное значение имеет то, что вихри, порожденные взаимодействием воздушных масс с поверхностью, имеют особенность диффундировать по направлению к атмосфере, т.е. концентрация их максимальна у поверхности почвы и постепенно убывает с высотой. Такое представление о подъемной силе объясняет как частицы с массой гораздо большей, чем у воздуха, поднимаются на значительные высоты.
Американцы применили для подъемной силы уравнение Жуковского для крыла самолета. В нем используется разность скоростей на верхнем и нижнем концах частицы. Для того, чтобы образовалась эта разность скоростей, необходимо, чтобы частицы имели обтекаемую форму как у крыла самолета. Даже если частица будет иметь такую форму, что уже мало вероятно, необходимо, чтобы выпуклой частью она была обращена к встречному движению ветра. Возможно, подобные частицы и окажутся на поверхности, расположенные определенным правильным образом, но, очевидно, что количество их будет ничтожно. Американцы, утверждают, что разность скоростей возникает в результате очень быстрого возрастания скорости ветра по профилю у поверхности почвы. И если рассматривать участок с максимально резким изменением скорости, то по отношению к размерам частицы, действительно окажется, что сверху и снизу скорость различна. Это так, но существенное значение этот факт может иметь только в самой близи поверхности, т.к. скорость достаточно быстро приобретает постоянное значение с высотой. Такое представление подъемной силы объясняет скачкообразные перемещения частиц у самой поверхности, но не объясняет всей масштабности явления ветровой эрозии.
Второй факт, вызывающий подозрение, это применение силы Магнуса в качестве силы отрывающей частицу и дающей ей ускорение при взлете и торможение при падении. Вообще, сила Магнуса является отклоняющей силой, т. е. она отклоняет частицу в различные стороны при полете. Сила Магнуса не может выступать в качестве подъемной силы по тем простым причинам, что, во-первых, она может возникнуть только когда частица уже находится в полете, но никак не на поверхности, а, во-вторых, доказано, что то количество вращений в единицу времени, которое наблюдается у частицы при отрыве от поверхности, не достаточно для подъема. Таким образом, американская теория не может объяснить подъема частиц на большие расстояния и заводит рассуждения в тупик. К тому же расчеты, предлагаемые на основе такого представления механизма подъема частицы, получаются очень громоздкими. А скорости вверху и у основания частицы, например, вообще невозможно измерить имеющимися на сегодняшний день приборами. То есть формулу для подъемной силы нельзя использовать при расчетах. Сомнительным оказывается применение формул для начальных вертикальной, горизонтальной и угловой скоростей, где используется формула для подъемной силы.
Что касается траектории частицы, то российская модель представляет общее уравнение траектории, достаточно простое в расчетах. Оно не может описать траекторию каждой частицы в отдельности, но, с другой стороны, является справедливым в отношении усредненной почвенной частицы и характеризует движение всех частиц определенного сорта. То же обстоятельство характеризует и формулы для высоты и длины скачка. Такой подход оказался очень удобным и правильным в том плане, что это позволяет вполне детально анализировать структуру почво-воздушного потока. Знание того, что частицы определенного размера (массы) ведут себя очень похожим образом и позволяет разделить почво-воздушный поток на континуумы.
Таким образом, уравнение траектории частицы, выведенное российскими учеными на основании рассмотрения сил, действующих на частицу, более содержательно и удобно в применении. В соответствии с этим уравнением при скоростях Uкр.< Uе< Uкр.2 почвенная частица движется скачкообразно; при скорости Uе = Uкр.2 частица в своем движении стремится к прямой горизонтальной линии; при Uе > Uкр.2 частица летит по вогнутой восходящей траектории. Из этого, становится понятен смысл критической скорости Uкр.2 для индивидуальных почвенных частиц. А поскольку в реальных условиях на поверхности почвы представлена не одна фракция частиц, а широкий спектр частиц различного размера, то полученное уравнение траектории автоматически делит почвенную фазу почво-воздушного потока на континуумы с собственными траекториями. Для баланса ветровой эрозии существенным является то, что они по типу движения частиц делятся на два потока - скачущих и летящих частиц.
Необходимым условием обобщения всех результатов, полученных для индивидуальных почвенных частиц, стало использование методов механики сплошной среды и представлений о взаимопроникающих движениях многоскоростных континуумов. Это позволило российским ученым доказать постоянство концентрации - , частиц i-того сорта вдоль i-той траектории, а это обстоятельство, в свою очередь, позволило успешно исследовать структуру почво-воздушного потока. В результате установлено, что реальный почво-воздушный поток, возникающий вследствие ветровой эрозии над бесконечным, однородным горизонтальным полем, имеет слоистую структуру, что обеспечивается различными типами движения частиц в потоке. Концентрация почвенной фазы - с, при переходе от почвенной поверхности к потоку убывает резко от единицы до малой величины , численно равной параметру массообмена (= B). Затем, в непосредственной близости от поверхности, она так же, скачком, увеличивается, вследствие сальтации, после чего, вплоть до границы слоя сальтации убывает с высотой. Выше этого слоя она остается постоянной. Чем больше сортов частиц вовлечено в скачкообразное движение, тем плавнее убывает концентрация с высотой. Детали изменения концентрации с высотой определяются особенностями распределения по размерам частиц в эродируемой почве и должны быть исследованы отдельно.
Помимо общепринятых методов исследования российские ученые использовали методы подобия и анализа размерностей при решении задачи о параметре массообмена, а при анализе структуры почво-воздушного потока использовали методы механики сплошной среды и представления о движении многоскоростных взаимопроникающих континуумов. Все это позволило российским ученым установить закономерности движения компонентов почво-воздушного поток на микроуровне, соответствующем размерам отдельных почвенных частиц, и на макроуровне, соответствующем масштабам всего почво-воздушного потока, которые в совокупности и составляют основы физически-содержательной теории ветровой эрозии почв.
Общепризнанной сейчас является вихревая структура потока, но до сих пор мало кто принимает во внимание столь важный факт. Например, объяснение американских ученых подъемной силы базируется на применении формулы для подъемной силы Жуковского для крыла самолета и силе Магнуса. Российские ученые в свое работе подразумевают неспособность подобных представлений объяснить механизмы ветровой эрозии. Они успешно строят свою модель на представлении о вихревой структуре потока и о вихревой природе сил отрывающих частицу от поверхности.
Вся механика сплошных сред строится на усреднении поведения отдельных индивидуумов. Поэтому российские ученые выводят уравнение траектории частицы при пыльных бурях не для конкретной частицы, как это делают американские ученые методом итерации, а для какой-то усредненной частицы определенного сорта.
В заключение надо сказать также, что американцы, предложив свои формулы для расчета начальных условий вылета частицы, зависимости концентрации и почво-воздушного потока от высоты, не предложили методов с помощью которых можно было бы измерить или посчитать такие, например, величины как: количество частиц, улетающих с единицы почвенной поверхности за единицу времени или разницу скоростей наверху и у основания частицы. Они предложили профиль скорости на высоте менее одного миллиметра, но не объяснили, как его измерили. Обе модели, тем не менее, интересны и, очевидно, необходимы дальнейшие исследования в этой области, поскольку хорошая, полная, удобная в расчетах и получении необходимых данных модель может послужить инструментом прогнозирования потерь почвы в составе общей геомиграционной модели.

2. Математическая программа ВЕТЭРО для моделирования ветровой эрозии

Для моделирования процессов ветровой эрозии почвы была написана программа "ВЕТЭРО", которая решает следующие задачи:
1) определение ореола загрязнения (площади покрытия радиоактивной пылью);
2) определение потерь ХВХ от ветровой эрозии;
3) определение концентрации пылесодержания на площади ореола загрязнения.
При решении каждой из этих задач возникают вопросы с выбором математического алгоритма реализации физической модели, который адекватно описывал бы каждый из этих процессов.
При решении первой задачи возникает вопрос о выборе закона распространения пыли. Если для определения расстояния пользоваться формулой прямолинейного переноса с учетом розы ветров, т.е. S = ?UT, где S - расстояние, U - скорость ветра, V= ?U - скорость частиц, ? - коэффициент замедления скорости частиц, T - время, то, как показывают расчеты, за сравнительно небольшой промежуток времени пыль переместится на значительное расстояние, что на самом деле вряд ли имеет место, т.к. существуют различного вида естественные преграды (дома, деревья, водоемы и т.д.). Также возникают вопросы об учете скачкообразного переноса пыли за счет порывов ветра, о неравномерном распределении пыли по поверхности и вторичном переносе за счет других факторов. Поэтому для решения первой задачи был выбран простейший алгоритм определения максимального ореола загрязнения без учета естественных преград, порывов ветра, неравномерного распределения и вторичного переноса (линейная модель).
Для решения второй задачи необходимы специальные исследования состава поверхностного слоя ХВХ и климатических условий. Такие исследования не проводились, поэтому использовалось уравнение определения потери почвы от ветровой эрозии Г.П. Глазунова, В.М. Гендугова с параметрами, полученными из различных литературных источников:
,
где Q - потеря почвы от ветровой эрозии, кг;
T - время прогноза (потерь), с;
S - площадь эрозии, м2;
Uе - скорость ветра за пределами слоя шероховатости, м/с;
p1 и p2 - концентрации почвенных частиц, перемещаемые безвозвратно и скачкообразно, кг/кг;
? - касательное напряжение на почвенной поверхности, возникающее вследствие воздействия на нее ветра, Н/м2;
Bк - параметр массообмена, характеризующий почву при Ue = Uкр., кг/кг;
Uкр. - критическая (размывающая) скорость ветра для данной почвы, м/с;
? - безразмерная почвенная константа.
При уточнении оценочных расчетов был принят ряд допущений:
* критическая скорость ветра (Uкр.) была принята равной 8 м/с;
* на основании усредненных метеоданных в расчете учитывали только дни, когда средняя скорость ветра (Uе) была выше Uкр. ;
* скорость ветра (Uе) в расчетные дни была принята равной среднему значению скорости ветра в эти дни;
* параметр массообмена (Bк) был принят равным промежуточной величине (средней величине для различных видов почвы) - 0.0000005 кг/кг;
* величина безразмерной почвенной константы (?) была принята равной - 8 (средней величиной из литературных данных);
* величины p1 и p2 были приняты 0,5 и 0,1 кг/кг, соответственно (0,4 кг/кг - вклад крупных непереносимых частиц);
* климатические условия: климат умеренно-континентальный, с небольшим количеством осадков. В летние месяцы, июль-август, приходится период максимума суточных температур, самые низкие температуры в январе-феврале. Атмосферные осадки распределены очень неравномерно. Наибольшее их количество выпадает в виде дождя и приходится на весну и осень, минимальное на летние месяцы, сентябрь. Среднегодовое количество 380 мм. Снег выпадает в декабре - феврале, толщина покрова до 15-20 см. Снеговой покров весьма неустойчив, держится 50-70 дней;
* в расчетах оценочно учитывали влияние снежного и травяного покрова на эрозионной площади, а также увлажнение почвы в результате дождевых осадков, при этом введен поправочный коэффициент -? в уравнение (1):
Qпопр.=?Q,
где
? = 0 - в случае наличия снежного покрова;
? = 1 - в случае "голой", "сухой" почвы;
? = 0,5 - в случае выпадения осадков и установления травяного покрова.
Расчеты проводили для трех "критических" вариантов:
1). без учета климатических условий, но для дней с превышением критической скорости ветра Uкр.;
2). для эрозии почвы (медленное расползание зеркала ХВХ) с оценочным учетом климатических условий;
3). для пылевой бури (распыление ХВХ).
3. Численные расчеты
Исходные данные и результаты оценки потерь почвы приведены в табл. 1-3.


Таблица 1. Результаты расчетов размывания ХВХ без учета климатических условий.
Месяц
p1
p2
p1/(p1+p2)
?, Н/м2
Кол-во дней с Uе ? 8 м/с
Uе в учетные дни, м/с
Потеря Q, кг/м2
Потеря Q, мм
Январь
0,5
0,1
0,833
0,320
2,5
8,1
0,004
0,003
Февраль
0,5
0,1
0,833
0,590
5
11,0
0,418
0,321
Март
0,5
0,1
0,833
0,568
17,5
10,8
1,226
0,943
Апрель
0,5
0,1
0,833
0,527
15
10,4
0,717
0,552
Май
0,5
0,1
0,833
0,487
15
10,0
0,469
0,361
Июнь
0,5
0,1
0,833
0,600
20
11,1
1,820
1,400
Июль
0,5
0,1
0,833
0,547
15
10,6
0,873
0,671
Август
0,5
0,1
0,833
0,855
10
13,25
3,750
2,885
Сентябрь
0,5
0,1
0,833
0,395
5
9,0
0,042
0,033
Октябрь
0,5
0,1
0,833
0,955
2,5
14,0
1,343
1,033
Ноябрь
0,5
0,1
0,833
0,395
7,5
9,0
0,063
0,049
Декабрь
0,5
0,1
0,833
0,702
5
12,0
0,896
0,689
Суммарные потери за год
11,621
8,939






Таблица 2. Результаты расчетов размывания ХВХ с учетом климатических условий.
Месяц
p1
p2
p1/(p1+p2)
?, Н/м2
Кол-во дней с Uе ? 8 м/с
Uе в учетные дни, м/с
?
Условия
Потеря Q, кг/м2
Потеря Q, мм
Январь
0,5
0,1
0,833
0,320
2,5
8,1
0
снежный покров
0
0
Февраль
0,5
0,1
0,833
0,590
5
11,0
0
снежный покров
0

страница 1
(всего 5)

ОГЛАВЛЕНИЕ

След. стр. >>

Copyright © Design by: Sunlight webdesign