УДК 551
ОЦЕНКА ПЛОТНОСТИ ОБЪЕМНОГО ЗАРЯДА ПО РЕЗУЛЬТАТАМ ИЗМЕРЕНИЙ ВЕРТИКАЛЬНЫХ ПРОФИЛЕЙ ЭЛЕКТРИЧЕСКИХ ХАРАКТЕРИСТИК В ПРИЗЕМНОМ СЛОЕ АТМОСФЕРЫ
© 2010 г. Г.Г. Петрова
Педагогический институт Pedagogical Institute
Южного федерального университета, of Southern Federal University,
пер. Днепровский, 116/2, г. Ростов-на-Дону, 344065, Dneprovskiy Lane, 116/2, Rostov-on-Don, 344065,
geogpu@rambler.ru geogpu@rambler.ru
Вопрос о закономерностях распределения плотности объемного заряда в приземном слое актуален в силу важности данного параметра для исследований атмосферного электричества и недостатка данных о нем, обусловленного сложностью измерений этой величины. На основании результатов измерений вертикальных профилей полярных электропроводностей и характеристик электрического поля атмосферы в слое 0-3 м в летний период сделаны оценки плотности объемного заряда косвенными методами в условиях устойчиво стратифицированного приземного слоя. Несмотря на приблизительность таких оценок, удается обнаружить некоторые типичные особенности вертикального распределения этой величины, допускающие их вполне разумную физическую интерпретацию.
Ключевые слова: атмосферно-электрические измерения, электродный эффект, объемный заряд, электропроводность, турбулентная диффузия, радон-222.
The question of the space charge density distribution in the surface layer is topical owing to the importance of this parameter in the atmospheric electrical research and in view of lack of the experimental data about it due to sophistication of the measurements of this value. On the basis of measuring results of both polar conductivities and electric field characteristic vertical profiles in the lowest 3-meter layer in summer time the space charge density is estimated by the indirect route for stable surface layer stratification. In spite of the approximateness of such estimations some typical features of the vertical distribution of this value are revealed.
Keywords: atmosphere electrical measurement, electrode effect, space charge, electric conductivity, turbulent diffusion, radon-222.
Состояние проблемы
Значение приземного слоя атмосферы для глобальных электрических процессов велико, поскольку именно здесь происходит обмен зарядами между атмосферой и землей. Здесь расположен один из электрических «контактов» глобальной электрической цепи (ГЭЦ): контакт нижнего слоя атмосферы с зем-
ной поверхностью. В то же время поиск закономерностей пространственно-временных вариаций электрических характеристик приземного слоя и их интерпретации осложнен тем, что электрические процессы протекают здесь вблизи границы двух сред, каковой является земная поверхность, и многие физические характеристики претерпевают на ней разрыв.
Теоретические подходы к исследованию структуры электродного слоя атмосферы подробно рассматриваются, например, в [1, 2], а также в справочном обзоре [3]. Для описания электрических процессов в нижних слоях атмосферы приводится система уравнений:
dn±
~~дГ
+ div(n± V) + div(n±b± E) — divD'± gradn±
= v± (r) - a'±n'+ n- + K± , divE =
P
rot E = 0.
(1)
Если j = AE + pV - DT Vp, где
■ + div j = 0.
(2)
При T >>
\nA
(A0 = en0b0) будет иметь место
чают уравнения, описывающие классический электродный эффект.
Характерная толщина электродного слоя в случае классического электродного эффекта определяется Ь0 Е0
как L =
где b0 = max{b+ ,b-}; E0 - напряжен-
Л = Х (е+Ь+п+ + е'_Ь'_ п_) = Л++Л_ - удельная электро-
г
проводность атмосферы, и р = е2 (п+ _ п_) - плотность объемного заряда ионов всех подвижностей, то
Р
Здесь е0 - диэлектрическая проницаемость вакуума; п'± - объемная концентрация ионов г -й группы; Ь'± - их подвижность; V - гидродинамическая скорость;
= + БТ - коэффициент диффузии; Д± - коэффициент молекулярной диффузии; Вт - коэффициент
турбулентной диффузии; К± - член, описывающий взаимодействие ионов с ионами других групп и с аэрозольными частицами; у'+ - интенсивность ионообразо-вания ионов г -й группы; аг± - коэффициенты рекомбинации ионов; ] - плотность электрического тока.
Вместе с граничными и начальными условиями система уравнений образует полную систему уравнений для нахождения пространственного распределения электрических характеристик атмосферы вблизи поверхности земли.
Обычно задачу сводят к одномерной, исходя из того, что характерные масштабы изменения электрических величин по горизонтали много больше вертикальных, а вектор напряженности перпендикулярен к земной поверхности. При этом находят вертикальные распределения Е, X, и р в горизонтально однородном приземном слое над горизонтальным участком земной поверхности.
Ь
ность вблизи земной поверхности.
Для учета аэрозолей в правой части уравнения для п± сохраняется член К±, учитывающий взаимодействие легких ионов с аэрозольными частицами. Обычно он имеет вид К± = п± _ у№0п± (р и у -коэффициенты присоединения легких ионов к заряженным и незаряженным аэрозолям; Ы± и Ы0 - концентрация заряженных и незаряженных аэрозолей соответственно).
Для учета влияния турбулентной диффузии на электродный эффект используют метод Рейнольдса. Все входящие в основную систему уравнений величины, например, /, представляются в виде
/ = (/) + /', где /' - флуктуации величины / ; 1 Т1
(/) = — | /& - среднее значение /, полученное путем Т1 0
осреднения по времени. Здесь т1 - масштаб времени, для которого справедливо условие Тф << Т << Тср , причем т - характерный временной масштаб флук-туаций; Тср - характерный временной масштаб для средних значений; (/') = 0. Если пренебречь квадратичными членами и представить (п'± V ^ путем введения коэффициента турбулентной диффузии ВТ в
виде
in'lv ) = -Dt(r)grad(n'±)
можно получить
стационарность процессов, и стационарную систему уравнений для горизонтально-однородного слоя можно записать, исходя из условия: ёп'^/л = 0 , где Ь -характерный масштаб пространственного изменения величин; Т - характерное время электрических процессов в приземном слое атмосферы.
При решении задачи об электродном эффекте рассматриваются два предельных случая: классический и турбулентный.
Классический электродный эффект проявляется в условиях хорошей погоды, т.е. при отсутствии ветра, облачности, конвекции, загрязнения. Считая коэффициент турбулентной диффузии равным нулю, а также пренебрегая молекулярной диффузией и аэрозолями и рассматривая только одну группу легких ионов, полу-
уравнение, описывающее турбулентный электродный эффект. Для Т >> т (Т - характерное время изменения электрических характеристик, т = 4а ) система уравнений становится стационарной. Для решения системы уравнений в случае турбулентного электродного эффекта требуется задание вида коэффициента турбулентной диффузии БТ (г) [3].
Обзор теоретических подходов показывает, насколько непростой является задача математического моделирования электродного слоя атмосферы. Система уравнений, описывающих электродный слой атмосферы, аналитического решения в общем виде в настоящее время не имеет, поэтому при решении обсуждаемой задачи как аналитически, так и широко применяемыми ныне методами численного моделирования неизбежно приходится принимать ряд условий и ограничений, которые, по-видимому, не всегда соответствуют реальным состояниям электродного слоя. Принимающиеся в процессе моделирования значения констант также могут не соответствовать их реальным значениям.
Сложность математического описания электрических свойств приземного слоя атмосферы делает особенно важным экспериментальное изучение электродного слоя. Как видно из приведенных уравнений, значительную роль в приэлектродных процессах играет плот-
S
0
ность объемного заряда р, которая в приземной атмосфере очень изменчива и подвержена влиянию множества факторов. В одномерном случае вариации р определяются дивергенцией плотности вертикального тока (2), а сам заряд в соответствии с уравнением Пуассона (1) обусловливает вертикальный профиль градиента электрического потенциала вблизи земли.
Вопрос о закономерностях вариаций плотности объемного заряда в приземном слое актуален в силу важности данного параметра для исследований атмосферного электричества и недостатка данных о нем, обусловленного сложностью измерений этой величины. В поисках этих закономерностей в настоящей работе сделаны оценки плотности объемного заряда косвенными методами [4] в условиях устойчиво стратифицированного приземного слоя. Несмотря на приблизительность таких оценок, на наш взгляд, удается обнаружить некоторые типичные особенности вертикального распределения этой величины, допускающие их вполне разумную физическую интерпретацию.
Настоящую работу следует рассматривать как продолжение статьи [5]. Здесь использованы те же массивы данных об электрических характеристиках атмосферы, а метеорологические условия измерений представлены в [5]. Методы измерений и измеритель-
ный комплекс, использованный при создании этой работы, описаны в [6].
Результаты экспериментальных исследований и их анализ
Один из методов оценки средней для некоторого слоя плотности объемного заряда ионов всех подвиж-ностей основан на измерениях электрического потенциала на нескольких уровнях. Он позволяет оценить р по вертикальному профилю градиента потенциала на основании уравнения О = —рей, где О = %гай р -
градиент электрического потенциала атмосферы. Поскольку в случае горизонтальной земной поверхности характерные масштабы изменения электрических величин по горизонтали много больше вертикальных и вектор О перпендикулярен к земле, можно записать dG|dz = — ре0 , где G - проекция вектора градиента
на вертикаль. Схема подсчетов р дана в [6].
На рис. 1, 2 помещена диаграмма плотности р объемного заряда ионов всех подвижностей, оцененной по измеренному профилю потенциала.
Результаты получены посредством одного и того же измерительного комплекса в разные периоды в двух пунктах Ростовской области: на измерительной площадке Михайловка в августе 1997 г. (рис. 1) и
Электропроводность , фСм/м LP-LO, фСм/м 5 10 15 20 25 -4 -2 0 2 4
0 20 40 60
20
40
60
-150
-50
50
г-
150
250
350
450
550
Плотность, объемного заряда, пКл/м
Рис. 1. Распределение по высоте электрических характеристик приземного слоя в условиях устойчивой стратификации температуры: концентрации радона-222 (1), положительной (2) и отрицательной (3) удельных электропроводностей атмосферы, их разности Я+ — (5), градиента электрического потенциала атмосферы (4) и плотности объемного заряда ионов всех групп (серая линейчатая диаграмма). Параметр стратификации m<0,7
Электропроводность, фСм/м 10 15 20 25 30
0
20
40
60 —I
LP-LO, фСм/м '-1-1—
-2 0 0 0 2 0 4 0 Градиент потенциала, В/м
я н о
<J
3 со
40
60
80 -350
-250
-150
-50
50
150
250
Плотность, объемного заряда, пКл/м
Рис. 2. Распределение по высоте электрических характеристик приземного слоя в условиях устойчивой стратификации температуры: концентрация радона-222 (1), положительная (2) и отрицательная (3) удельные электропроводности атмосферы, их разность Л+ — Л__ (5), градиент электрического потенциала атмосферы (4) и плотность объемного заряда ионов всех групп (серая линейчатая диаграмма). Параметр стратификации m<0,5
измерительной площадке Талловеров Кашарского района в августе 2008 г. (рис. 2). Различием является то, что в Талловерове потенциал атмосферы измерялся не на трех высотах: 1, 2, 3 м, а на четырех: 0,5, 1, 2, 3 м, что позволило подробнее просчитать градиент потенциала и, соответственно, плотность объемного заряда вблизи земной поверхности.
Выборка получена для 27 ч измерений и соответствует условиям устойчивой стратификации-! для пункта наблюдений в [5].
Здесь же представлены профили концентрации ра-дона-222, удельных полярных электропроводностей и градиента потенциала, измеренные одновременно в условиях устойчиво стратифицированного приземного слоя. Планками на всех рисунках показаны значения стандартной погрешности.
Для оценки плотности объемного заряда легких ионов можно использовать разность удельных полярных электропроводностей (Л+ —Л— ). Действительно, исходя из выражения для плотности объемного заряда легких ионов рли = в(п+— п—) с учетом того, что Л+ = вЬ+п+ и Л+ = вЬ_п_, полагая Ъ+ = Ъ_ = Ъ , получим рЛии = (Л+ — Л—)/Ъ (здесь e - элементарный заряд; п+, п_ - концентрация положительных и отрицательных легких ионов; Ь - их подвижность).
Вертикальный профиль разности (Л+ — Л_) помещен на рис. 1, 2. Отрицательные значения разности полярных электропроводностей (Л+ — Л— ) на графиках соотносятся с отрицательными значениями плотности объемного заряда ионов всех подвижностей, рассчитанной по результатам измерения характеристик поля.
Следует заметить, что подвижности положительных и отрицательных легких ионов разнятся. В таком случае для подсчета плотности объемного заряда легких ионов через удельные электропроводности надо
Л. Л_
использовать соотношение рл и = —---.
' ' Ъ+ Ъ_
Если принять, что Ъ+ = 1,36 10 4 и Ъ— = 1,56 -10 —4 м2/В-с (как обычно мы полагаем в наших работах), то даже в случае отрицательных величин (Л+ — Л_) значение плотности объемного заряда легких ионов, рассчитанное по данным эксперимента, может оказываться положительным. В то же время вид профиля рл и по экспериментальным данным повторяет конфигурацию профиля (Л+ — Л—) с впадиной в сторону отрицательных значений, указывая, как отмечалось выше, на процесс накопления отрицательного объемного заряда легких ионов, который, очевидно, и компенсирует в этом слое частично или полностью положительный объемный заряд электродного эффекта. В работе [7] получена зависимость вида расчетного профиля градиента потенциала в нижнем километровом слое атмосферы от соотношения подвижностей положительных и отрицательных ионов, которая косвенно подтверждает наши выводы. По результатам этих расчетов при соотношении подвижностей, близком к 1, реализуется 2-образный профиль поля, что указывает на наличие слоя отрицательного заряда, но
уже при Ъ+/Ъ— =0,9 профиль градиента потенциала при тех же прочих условиях становится гладким и соответствует положительному объемному заряду на всем протяжении исследуемого слоя.
Плотность объемного заряда, дислоцированного на ионах всех подвижностей, в том числе и на заряженном из-за адсорбции зарядов легких ионов аэрозоле, по нашим оценкам в среднем на порядок больше, чем плотность объемного заряда легких ионов. Действительно, для (Л+ — Л— ) = 1 фСм/м при Ь = 1,5• 10-4 м2/В^с рлм = =7 пКл/м3. Как видно из рис. 1, 2, (Л+ — Л_ ) достигает 4 фСм/м, т.е. рл и имеет величину порядка десятков пКл/м3. (Аналогичн^1й результат получен в работе [8], где приводятся профили как полярных электропровод-ностей в слое до 1,5 м, так и их разности.) Плотность объемного заряда ионов всех подвижностей оценивается сотнями пКл/м3 (рис. 1, 2).
На основании анализа результатов более ранних исследований и известных представлений можно предложить следующий ход рассуждений относительно механизмов формирования объемного заряда вблизи земли. На границе атмосфера-земля вертикальный электрический ток проводимости, обусловленный глобальным атмосферно-электрическим полем и ионизацией атмосферы, претерпевает разрыв из-за различной природы проводимости граничащих сред. Вследствие этого разрыва возникает электродный эффект, который проявляется в дефиците легких ионов одной из полярностей вблизи земной поверхности, что приводит к образованию в атмосфере объёмного заряда. В свою очередь последний способствует пространственному изменению напряженности поля Е в соответствии с уравнением Пуассона. Как известно, потенциал ионосферы положителен по отношению к потенциалу земной поверхности. При этом ток положительных ионов направлен вниз, к земле, а ток отрицательных ионов - от земли. В этом случае в приземном слое образуется объемный заряд положительного знака. Плотность заряда по мере удаления от земной поверхности уменьшается и на некоторой высоте обращается в ноль. Этот уровень определяет толщину электродного слоя атмосферы.
Причиной накопления объёмного заряда легких ионов в атмосфере может также служить прохождение тока проводимости через слой с вертикальным градиентом электропроводности. Согласно Чалмерсу [9], для квазистатических условий справедливо г dЛ
рЛ=е0——, где Л- удельная электропровод-Л dh
ность атмосферы; г - плотность вертикального тока проводимости. Как показывают исследования [10-13], вблизи земной поверхности часто наблюдаются вертикальные градиенты электропроводности, вызванные уменьшением с высотой интенсивности ионооб-разования, обусловленной естественной радиоактивностью. Таким образом, при прохождении тока проводимости вблизи земли образование объемного заряда легких ионов осуществляется как в результате электродного эффекта (+рээ), так и вследствие градиента электропроводности атмосферы (+рЛ), обусловленного уменьшением интенсивности ионообра-зования с высотой. Итоговое значение плотности объ-
емного заряда легких ионов на каждом горизонте электродного слоя рли =±рээ + определится знаком и величиной каждого слагаемого, которые будут зависеть от состояния слоя [10].
Время жизни легких ионов невелико, в то время как заряженный аэрозоль может долго находиться в атмосфере. По-видимому, в процессах формирования объемного заряда в приземном слое участвуют легкие ионы, так как эти процессы связаны с током проводимости, а уже готовый заряд, оседая на частицах аэрозоля, может длительно присутствовать в электродном слое, определяя его электрическую структуру. Аэрозоль в таком случае играет роль аккумулятора объемного заряда, способного перенимать его у легких ионов и накапливать, обусловливая наблюдаемые значения плотности заряда ионов всех групп. В итоге электрическое поле в приземном слое атмосферы претерпевает сложные изменения по вертикали из-за совокупного действия электродного эффекта и наличия слоев объемного заряда иного происхождения. При этом из-за присутствия в атмосфере заряженных слоев может меняться по высоте не только численное значение напряженности поля, но и ее направление.
В исследовании структуры электродного слоя J. Law [14] одним из первых обнаружил появление отрицательного объемного заряда вблизи земли в ночные часы. Law уверенно говорит о меняющихся по высоте в нижнем метровом слое концентрациях полярных ионов для большинства условий в Кембридже, связывая это с увеличением концентрации радиоактивных эманаций вблизи земли. Им обнаружен отрицательный объемный заряд на высоте 50 см ночью, тогда как днем заряд на этой высоте положителен. Различия, прослеживающиеся как по экспериментальным, так и по расчетным данным для ночного и дневного электродного слоя, автор [14] связал со значительным изменением турбулентной диффузии в течение суток.
W. Crozier [15] получил детальное распределение характеристик поля и плотности объемного заряда в электродном слое. Измерялся потенциал атмосферы в нижнем 3-метровом слое на 4 высотах пассивной проволочной антенной и плотность объемного заряда методом Оболенского на уровне земли и на высоте 1 м. Crozier обнаружил, что тихими ночами (скорость ветра на высоте 4 м ниже 0,5 м/с) в профиле градиента потенциала имеет место минимум на высоте 15-40 см. Ниже этого минимума присутствует положительный объемный заряд с плотностью 160-320 пКл/м3 на уровне земли. Над минимумом градиента потенциала располагается отрицательный объемный заряд до высоты, по крайней мере, 3 м с плотностью порядка сотен е/см3 (100 е/см3=16 пКл/м3). Нижний положительный объемный заряд Crozier объяснил электродным эффектом, а отрицательный заряд справедливо связал со значительным перепадом по высоте скорости ио-нообразования вследствие заметной диффузии радиоактивных газов из почвы.
Из вышеизложенного следует, что значительная эксхаляция радона и накопление его у земли при устойчивой стратификации атмосферы существенно сказывается на вертикальном распределении электропроводности в электродном слое, повышая ионизацию вблизи земной поверхности. Это, в свою очередь, влияет на процессы формирования приземного объемного заряда, а значит и на электрическое поле, и на обмен зарядами между землей и атмосферой. Из-за образования в электродном слое трехслойной структуры объемного заряда в условиях устойчивой стратификации формируется Z-образный профиль градиента потенциала вблизи земли, вертикальный масштаб которого зависит от устойчивости слоя.
Литература
1. Hoppel W.A. Electrode effect: comparison of theory and measurement // Planetary Electrodinamics. N.-Y., 1969. P. 167-181.
2. Куповых Г.В., Морозов В.Н., Шварц Я.М. Теория электродного эффекта в атмосфере. Таганрог, 1998. 124 с.
3. Морозов В.Н. Атмосферное электричество // Атмосфера: справочник (справочные данные, модели). Л., 1991. С. 394-408.
4. Имянитов И.М. Приборы и методы для изучения электричества атмосферы. М., 1957. 344 с.
5. Петрова Г.Г.Вертикальные профили градиента электрического потенциала в различно стратифицированном приземном слое атмосферы // Изв. вузов. Сев.-Кавк. регион. Естеств. науки. Спецвыпуск. Физика атмосферы. 2010. С. 47-81.
6. Измерительный комплекс для исследования электричества приземного слоя атмосферы / А.И. Петров [и др.] // Изв. вузов. Сев.-Кавк. регион. Естеств. науки. Спецвыпуск. Физика атмосферы. 2010. С. 47-52.
7. Anisimov S.V., Dmitriev E.M. Aeroelectrical altitude profile in limit of regular convection // J. Atm. Electr. 1999. № 19 (1). P. 25-33.
8. Israelsson S., Knudsen E., Anisimov S.V. Vertical profiles of electrical conductivity in the lowermost part of the turbulent boundary layer over flat ground. // J. of Atm. and Terr. Phys. 1994. Vol. 56, № 12. P. 1545-1550.
9. Чалмерс Дж.А. Атмосферное электричество. Л., 1974. 422 с.
10. Петров А.И., Петрова Г.Г. Вертикальные профили полярных электропроводностей и плотности объемного заряда в электродном слое атмосферы // Тр. ГГО. 1988. Вып. 514. С. 12-16.
11. Petrov A.I., Petrova G.G., Panchishkina I.N. On factors determining the variations of the electric characteristics of a surface layer // Proc. 11th Int. Conf. Atm. Electricity. Alabama, USA, 1999. P. 547-550.
12. Petrov A.I., Petrova G.G., Panchishkina I.N. Measurements of polar conductivities in the surface layer of the atmosphere // Proc. 10th Int. Conf. Atm. Electricity. Osaka, Japan, 1996. P. 188-191.
13. Электропроводность воздуха и концентрация радона в приземном слое / А.И. Петров [и др.] // Тр. V Рос. конф. по атмосф. электр. Владимир, 2003. Т. 1. С. 124-127.
14. Law J. The ionization of the atmosphere near the ground in fair weather. // Quart. J. R. Met. Soc. 1963. № 89. P. 107-121.
15. Crozier W.D. Atmospheric electrical profiles below three meters // J. Geoph. Res. 1965. № 70. P. 2785-2792.
Поступила в редакцию_23 ноября 2009 г.