Атмосферное электричество и качество воздуха

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

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

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

1.Распределение зарядов частиц для стационарной аэродисперсной системы при симметричной биполярной зарядке описывается больцмановской функцией вида [ 1 ]

f (и) = ^■e%p[—(ne)2/87T£or/cT], где – зарядэлектрона; – электрическая постоянная; –радиус частицы; -постоянная Больцмана; -абсолютная температура;

-положительные и отрицательные целые числа или нуль; -суммарное число частиц, получающееся как результат суммирования по всем значениям п.

70÷ 90% частиц в тропосфере имеют размеры от 10 -7 до 5 - 10 -6 см . [ 2 ] . Частицы радиусом 0,01 мкм большей частью (90%) нейтральны и только 10% частиц имеют заряд ± . Частицы радиусом 0,05 мкм на 43% нейтральны, 48% частиц имеют заряд ± , 8,6% - ±2 и 0,4% - ±3 . Частицы так называемого комнатного аэрозоля крупнее (около 0,1 мкм ), а доля заряженных частиц среди них больше. Из практики давно известно, что заряд субмикронных частиц обычно составляет 1 e . Отметим также, что напряженность электрического поля у поверхности частицы, несущей заряд 1 е и имеющей радиус 0,1 мкм , равна приблизительно 1,6 кВ-см "\ то есть значительно меньше границы ионизации (около 5 кВ-см ^-1) [ 3 ] , что и подтверждает наблюдаемый единичный заряд частицы.

Кинетика зарядки хорошо исследована аналитически во многих работах. Кирш и За-гнитько [4] обосновали следующее соотноше ние: ()=4    3 ̅      ()+

(

1+2     ⧸( +2)       ⧸(4   +) где: -средний радиус частицы;   -средний квадрат радиуса частицы; - напряженность электрического поля у поверхности частицы;

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

Уравнение стационарного состояния ионного баланса при условии отсутствия адвективного (посредством ветра), электрического и турбулентного переноса ионов записывается в виде [5] :v = a±n+п - + ±±n±Nа, где ν –скоростьобразования ионов. Первое слагаемое справа описывает рекомбинацию ионов (  ,   -концентрации ионов, ± - коэффици енты рекомбинации). Второе слагаемое справа описывает процесс конверсии ионов в аэрозольные частицы или адсорбции ионов на аэрозольных частицах, включая составляющие (адсорбция положительных ионов на нейтральной частице),         (адсорбция от рицательных ионов на нейтральной частице), (адсорбция положительных ионов на отрицательно заряженной частице) и другие слагаемые. Использовались известные [6] экспериментальные значения конверсионных сла- гаемых, в частности,     =1,7∙10 см ∙ с ,    =2,34∙10 см ∙с ,    =5,22∙

10 см ∙с , а также и всех других слагаемых. Среднее значение для всех составляющих получается равным =3,4∙ 10 см ∙с . Запишем далее уравнение стационарного состояния при условии     =   = в следую щем виде:

+2    - =0.

Элементарный анализ этого уравнения выражается в следующих оценках. Низкая концентрация частиц:

≫ 1, откуда =( ⁄ ) ^⁄ .

Подставляя фоновые значения [ 7 ] ν= 10см ∙с и =10 см ∙с , получаем =10 см .

Высокая концентрация частиц: ⁄      ≪ 1, откуда = ⁄2   .

Подставляя значение =3,4∙10 см ∙ с , получаем условие >10 см . Полученные оценки имеют следующий физический смысл. Значение, полученное для концентрации ионов, то есть =10 см при   =0, типично для свободной атмосферы и соответствует эволюционно сложившейся физиологической норме. Это значение примерно соответствует ионизации над поверхностью океана. Это же значение соответствует принятой в РФ гигиенической норме (биологическая единица аэроионизации в естественных условиях, БЕ [8]). Для взрослого человека это означает вдыхание в течение суток 8∙10 легких ионов, а суточной нормой для здорового взрослого человека считается (2-3)БЕ. Хорошо известно, что снижение концентрации легких ионов до 50см –3 вызывает головокружение, тошноту и общую слабость.

При высокой, но все же умеренной концентрации частиц    =10 см (при этом массовая концентрация частиц радиуса 0,1мкм равна примерно 0,004мг-м-3, а концентрация ионов =1,37∙10 см ) присутствие частиц заметно не изменяет естественную ионную компоненту.

Постоянная времени =1⁄2 определяет скорость ионной конверсии при заданных значениях концентрации ионов и производительности ионного источника. Принимая =10 см    и  ν=10см  ∙с , получаем

=100 с.

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

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

2.Для измерения спектра подвижностей высокоподвижных частиц (ионов) используются аспирационные конденсаторы, имеющие чаще всего цилиндрическую конструкцию. Пример такой конструкции [ 9 ] , используемой для атмосферных исследований, показан на рис.1.

Анализируемый воздух (sampleair) через кольцевую щель, образованную направляющей поток острой кромкой (knifeedge), течет в тонком слое вдоль внутренней поверхности внешнего цилиндрического электрода, имеющего нулевой (ground, GND) потенциал. Этот поток стабилизирован потоком чистого, сухого воздуха (sheathair), подаваемого через сетчатый фильтр (gauze), служащий для получения ламинарного течения.

Центральный положительно заряженный собирающий электрод (centralrod) имеет высокий (до 12000 В) потенциал (hightension, HT). Отклоняемые в электрическом поле отрицательно заряженные частицы, имеющие при заданном потенциале центрального электрода определенную подвижность, осаждаются на собирающем электроде. Неосажденные частицы удаляются через выходную щель (exitslit). Устройство заключено во внешний кожух (outercasing), электрически хорошо изолированный от нулевого электрода. Соотношение производительностей 1∶ 10 анализируемого потока (1 л/мин) и стабилизирующего потока (10 л/мин) необходимо для задания направле- ния дрейфа частиц к собирающему электроду уже на входе в конденсатор.

Sample out

Рисунок 1 – Цилиндрическийаспираци-онный конденсатор

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

( )≥   ( ⁄ ) 2    (   -   ), где: ;  - радиусы соответственно внутренне го (центрального) и внешнего электродов конденсатора; -время пребывания частицы в конденсаторе, определяемое скоростью невозмущенного (ламинарного) потока и длиной в направлении движения; -приложенное напряжение.

Тогда уравнение, связывающее вольт-амперную характеристику конденсатора с функцией плотности вероятности распределения подвижностей ( ) может быть записано для безразмерного тока в виде ( )

()1

=1-   [ ( )]-( )

н                             ()

Vо где ( )- функция распределения.

Значение н определяет максимальный ток или ток насыщения. Нетрудно показать, что полученное соотношение удовлетворяет в интервале напряжений 0≤  ≤    очевидным условиям (0)=0 и (    )= н.

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

В качестве наиболее общей аналитической формы может быть выбрано гамма – распределение: ( )=  ⧸             ,где а»;в>0;k>0, причем форма кривой определяется параметром а, а параметр в носит характер масштабного множителя; ( )- гамма-функция Эйлера. В зависимости от величины а функция ( ) видоизменяется от экспоненциального распределения ( =1) до почти сим метричного (нормального) распределения с ростом а.

Математическое ожидание и дисперсия подвижности даются в следующей простой форме:

т- аа

=  ; ()=.

Функция распределения ( ) дается выражением

()= ()

где ( )- неполная гамма-функция.

Используем теперь для определения параметров распределения начальный ( конечный (асимптотический) (  →     ) участки вольтамперной характеристики.  При

→0 и, соответственно, при ( )→     (то есть для наибольшего измеряемого значения подвижности) получаем приближенное равенство

• ⁽⁰⁾ _н ≅  ⁄ (  )=¹   _{(  )}      .

Введем обозначение (0)⁄ н = -произвольное, но малое значение безразмерного тока; для значения ( )=      перепишем предыдущее равенство в виде

51 — 1/ктах i I^ , из которого сразу следует возможность оценки величина математического ожидания подвижности на основании соотношения:

=

н у=о

Коэффициент определяется выражением =1 2  (  -  ), содержащим только параметры аспирационного конденсатора.

Аналогично, для конечного участка вольтамперной характеристики при

(( )→ н), то есть для значения ()= к0 получаем соотношение

( )

• 1-        н=,

которое после некоторых преобразований приводит к соотношению

()

≅      ⧸ (+2)

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

Приведенные формулы справедливы для часто рассматриваемого диапазона подвижностей 0,01 ≤ k ≤2 (в единицах см² / В ∙ с ).

Следовательно, необходимо использовать значения = 0,01 и      =2. Таким образом, если имеется экспериментально полученная вольтамперная характеристика и, следовательно, значения и , то можно определить параметры распределения и , а затем математическое ожидание и дисперсию подвижности, а также и все остальные параметры распределения (коэффициент вариации, коэффициент асимметрии и другие).

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