Елохин Автоматизированные системы контроля радиационной обстановки окружаюсчей среды 2012
.pdf
Рис. 8.6. Иллюстрация к определению мощности выброса радиоактивной примеси по коэффициенту отражения электромагнитной волны ||R||
Таким образом, из (8.19) действительно следует | R | Pв , что
и позволяет по измеренному коэффициенту отражения определять радиационные характеристики радиоактивного выброса.
8.2.1. Учет метеорологического состояния реальной атмосферы при решении задачи определения мощности источника радиоактивных выбросов
Рассматривая ионизационные образования в атмосфере, т.е. в смеси газов, в которой кислород составляет 21 %, еще раз подчеркнем, что они, в основном, должны представлять ион-ионные, но не электрон-ионные образования. В основе этого утверждения лежат следующие соображения.
Как показано в п. 8.1.2, время жизни свободного электрона в воздухе до его захвата молекулой (атомом) кислорода составляет
τз 1 (kзNO2 ) ≈ 2 10-9 с, а время электрон-ионной рекомбинации
(при условии n = N + = |
GI k |
ei |
≈ 2,4 106 |
см-3) составляет τp ~ |
e |
|
|
|
1/keine ≈ 4 с. Из сравнения τз, τp видно, что электрон после «рождения» в результате акта ионизации атома или молекулы и сброса своей энергии до тепловых значений через 2 нс захватывается нейтральным атомом или молекулой кислорода с образованием отрицательного иона. Система уравнений (8.2), как показано в п. 8.1.2, при равных нулю производных сводится к уравнению третьего порядка:
251
|
|
kei N+3 + kзN0kei N+2 −keiGIN+ −GIkзN0 kei |
|
kii |
= 0, |
|
|
|
|||||||||
решение |
|
которого |
|
|
и |
определяет |
|
|
концентрации |
||||||||
n |
= GI −k |
ei |
N 2 |
|
k |
N |
0 |
1−k |
ei |
k |
=1,036.10-3, N |
– |
= N |
– n |
e |
= |
|
e |
|
+ |
|
з |
|
( |
|
ii ) |
|
+ |
|
|
|||||
= 6,021 105, |
|
N+= 6,021 105см-3. Такое существенное различие в кон- |
|||||||||||||||
центрации ионов и электронов обусловлено именно тем, что кислорода в 1013 раз больше, чем положительных ионов даже при высокой мощности дозы (I ~ 10-2 Зв/ч). Поэтому при длительности локации τ ~ 2−3 мкс в реальной атмосфере и тем более в ее пограничном слое, где и формируются ионизированные образования – плазмоиды, необходимо рассматривать их как ион-ионные образования, но не как электронные. В земных условиях рассматривать электрон-ионные плазменные образования возможно в том случае, если увеличить время жизни свободного электрона, например, путем уменьшения концентрации кислорода в воздухе за счет уменьшения давления атмосферы, что достигается в ее верхних слоях, но не рассматривается автором. С учетом найденных значений концентраций электронов и ионов оценим изменения диэлектрической проницаемости воздушной среды в области радиоактивного загрязнения подстилающей поверхности [17]:
|
4πn e2 |
|
|
4πN |
+ |
e2 |
|
|
|
|
|
4πN |
− |
e2 |
|||||
ε =1− |
e |
|
− |
|
|
|
|
|
|
− |
|
|
|
|
≈ |
||||
me (ω2 +ν2 ) |
M+ (ω2 +ν2 ) |
M− (ω2 +ν2 ) |
|||||||||||||||||
|
≈1− |
|
4πe2 N |
i |
|
|
|
2 + |
N |
M |
i |
|
|
|
|
||||
|
|
|
|
|
|
|
e |
|
|
, |
|
|
|
||||||
|
|
2 |
2 |
) |
|
|
|
|
|
|
|
|
|
||||||
|
|
Mi (ω +ν |
|
|
|
me Ni |
|
|
|
||||||||||
где ν – частота соударений, ω – частота электромагнитной волны,
М+ = M- = Mi – масса ионов, те – масса электрона и N+ = N- = Ni. Принимая в качестве массы иона молекулу кислорода, из послед-
него выражения находим, что влиянием электронов на изменение диэлектрической проницаемости можно пренебречь, поскольку при найденных значениях ne, Ni отношение Mine/теNi ≈10-4 << 1. Полученное неравенство остается справедливым и при массе тяжелых ионов равных 103 Mi.
Поскольку формирование плазмоидов осуществляется в реальной атмосфере, характеризующейся определенной влажностью и дисперсным составом ионов [21], постольку эти особенности среды обязательно должны быть учтены при описании их физической
252
структуры, и необходимо внести определенную ясность в условия их устойчивости. Действительно, наличие определенной влажности
врегионе, в котором наблюдаются эти явления (предполагается, что обязательно имеет место радиоактивное загрязнение подстилающей поверхности или происходит выброс радиоактивной газоаэрозольной примеси, или выход радиоактивного газа, например радона, из разломов земной коры как предвестник землетрясений [22]) должно приводить к образованию тяжелых ионов, подвижность и коэффициент диффузии которых значительно меньше легких (см. табл. 7.2), что обязательно должно сказываться на условиях переноса этих ионов и поперечной устойчивости плазмоида
вцелом. Кроме того, на тяжелые ионы, в отличие от легких аэроионов, действуют дополнительные силы – силы тяжести, которые приводят к пространственному разделению зарядов.
Идея разделения тяжелых ионов, из которых большей массой
обладали отрицательные, под действием силы тяжести позволила Я.И. Френкелю объяснить причину электрической поляризации облаков и их механическую устойчивость при наличии восходящих воздушных потоков [23]. Также следуя результатам этой работы полагаем, что в поле действия источника γ-β-излучения, т.е. в приземной области, загрязненной при выпадении аэрозолей после прохождения радиоактивного облака, в результате ионизации воздуха происходит образование положительных и отрицательных ионов, которые при наличии влажности образуют гидратированные комплексные ионы того и другого знаков [21].
Рассмотрим движение водных капель массой m1, m2 с радиусами rl, r2 (r2 >> r1), имеющих электрический заряд соответственно q+ и q-. На движущуюся каплю действует сила тяжести mg, сила Стокса Fs= 6πνr1w и сила электростатического взаимодействия Fe = Eq, где g – ускорение силы тяжести, ν – динамическая вязкость воздуха, w – скорость относительного движения, E – напряженность электрического поля, возникающего в результате поляризации ионизированного слоя. Движение каждой капли в поле действия этих сил иллюстрируется рис. 8.7 и описывается уравнениями:
m |
dw |
= m g −6πνr w − Eq− ; |
m |
dw |
= m g −6πνr w + Eq+ , |
||
2 |
dt |
2 |
2 |
1 |
dt |
1 |
1 |
253
Рис. 8.7. Иллюстрация к вопросу о поляризации водного облака, образовавшегося в атмосфере при захвате ионов каплями воды или микроскопическими частицами
пыли
решение которых при начальной скорости капли, равной нулю, имеет вид:
w |
= |
m g − Eq− |
|
−exp |
|
− |
6πνr |
t |
|
; |
2 |
1 |
|
2 |
|
||||||
2 |
|
6πνr2 |
|
|
|
m2 |
|
|
||
|
|
|
|
|
|
|
|
|
||
w |
= |
m g + Eq+ |
|
−exp |
|
− |
6πνr |
|
|
. |
1 |
1 |
|
1 t |
|
||||||
1 |
|
6πνr1 |
|
|
|
m1 |
|
|
||
|
|
|
|
|
|
|
|
|
||
|
Для установившегося движения |
|||||||||
при |
r2 = 10-3см и |
|
|
t >> τmax = |
||||||
= m 6πνr = 2ρr 9ν ≈10−3 c, |
|
т.е. |
||||||||
|
2 |
2 |
2 |
|
|
|
|
|
|
|
когда ω1 = ω2, величина электрического поля при m = (4/3)πρ r3 будет равна
E = g m2r1 −m1r2 = |
4 πρgr r |
r2 −r1 |
, |
||
|
|||||
q(r2 + r1 ) |
3 |
1 2 |
q |
|
|
|
|
|
|||
которое при r1 = 10-4 см и r2 = 10-3 см, ρ = 1 г/см-3 и q = e = 1,6.10-19 K, составляет величину ~2,5.105 В/см, тогда как напряженность электрического поля Земли составляет ~1,3 В/см. Заметим, что с увеличением влажности атмосферы величина напряженности электрического поля пробоя также будет расти [24]. Скорость коллективного
установившегося |
движения |
|
капель равна |
w = |
2ρg |
r2 |
−r r |
+ r2 |
|
|||
|
|
2ρg |
|
|
|
|
9ν |
2 |
1 2 |
1 |
|
|
или при r2 >> r1 |
w = |
r2 |
, |
которая при заданных значениях па- |
||||||||
|
||||||||||||
|
|
9ν 2 |
|
|
|
|
|
|
|
|
||
раметров и ν = 1,94.10-5 кг/мс составляет 1 см/с.
Таким образом, в радиационном поле излучения источника образующиеся аэроионы при наличии влажности или пыли переходят в тяжелые ионы, формирующие водные или пылевые облака, которые поляризуются под действием силы тяжести этих ионов. Убедительным аргументом в пользу существования такого механизма может служить экспериментальное измерение разности электропроводности δλ = (λ+ – λ–) как функция высоты, проведенное в работе [25] и приведенное на рис. 8.8.
254
Рис. 8.8. Высотный профиль разности электропроводности воздуха δλ= (λ+ – λ–). Полигон КЭ ИЭМа в Заполярье в 7 км к югу от аварийного блока ЧАЭС [25]
При наличии восходящих воздушных потоков, скорость которых может составлять ~ 1 м/с, согласно [23], аэроионы будут увлекаться этим потоком, захватываться каплями воды, образуя тяжелые преимущественно отрицательные ионы, которые, разделяясь под действием силы тяжести, наряду с положительным зарядом мелких капель и микроскопическими частицами пыли, адсорбирующими положительный заряд из воздуха, будут создавать слоистую структуру, продольная устойчивость которой будет обусловлена как электростатическим взаимодействием между слоями, так и постоянством восходящего воздушного потока.
Оценим величину силы электростатического взаимодействия между слоями, рассматривая подобную систему как плоский конденсатор с поверхностным зарядом σ = εε0Е, где Е – найденное ранее значение электрического поля и ε0 = 8,85.10-12 Ф/м. В условиях весомых «пластин» конденсатора энергия их электростатического взаимодействия CU2/2, где C – емкость конденсатора, U – разность потенциалов «пластин», должна уравновешиваться работой внешних сил Fed, где d – расстояние между»пластинами», не позволяющими им «схлопнуться». Из равенства CU2/2 = Fed находим Fe = CU2/2d.
255
Для плоского конденсатора C = εε0S/d, a U = Q/C, где Q = σS. Таким образом, Fe = εε0SE2/2 и при ε = 1 (для воздуха) и S = πr02 (r0 = 5 м, r0 – радиус плазмоида) Fе = 2,17 105 Н. Очевидно, что с ростом размеров капель эта сила будет возрастать, поскольку напряженность электрического поля будет расти как квадрат радиуса
капли. При скорости воздушного потока U0(h) найдем составляющую силы Fc, разрывающую плазмоид, при воздействии на него воздушным потоком (рис. 8.9).
При давлении воздушного потока P = ρaU02(h)/2, где ρa – плотность воздуха, сила равна
Рис. 8.9. Иллюстрация к вопросу о продольной устойчивости плазмоидов
m = 0,5 [26], находим
Fc = |
ρa ∫h U02 (h)ds, ds = 2r0dh. |
|
2 0 |
Аппроксимируя U0(h) выражени-
h m
ем вида uz0 z0 при z0 = 10 м и
|
ρa |
|
2 |
|
|
h |
2m |
|
h |
|
|
|
|
|
1 |
H |
|||
|
|
|
|
|
|
|
|
|
∫ Fc (h)dh, |
||||||||||
Fc = |
|
uz0 |
2r0 |
|
|
|
|
|
|
, Fc |
= |
|
|
||||||
2 |
|
|
|
2m +1 |
|
H |
|||||||||||||
|
|
|
|
z 0 |
|
|
|
|
|
|
|
|
0 |
||||||
что при заданных параметрах и Н = 500 м дает |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
H |
|
2m |
|
|||
|
|
|
|
Fc |
= |
|
|
ρauz0 r0 H |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
. |
|
|||||||
|
|
|
|
|
(2m +1)(2m |
|
|
||||||||||||
|
|
|
|
|
|
|
|
+ 2) z0 |
|
|
|
|
|||||||
Последнее при u |
z0 |
= 3−5 м/с составляет (2,25−6,25) 105 Н. При от- |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
клонении плазмоида на угол β << 10° составляющая Fc, действующая на разрыв плазмоида, равна Fcsin(10°) ≈ 3,9 04 −1,09 105 Н, т.е. представляет собой величину, либо значительно меньшую, чем Fс, либо сравнимую с ней. Но с ростом скорости ветра будет, очевид-
но, увеличиваться и угол наклона, так что при β = 15° Fc0 sin (β) ≈
≈ 1,62 105 Н.
Таким образом, если скорость ветра такова, что угол, на который отклоняется плазмоид от вертикального положения не превы-
256
шает некоторое предельное значение βпр ≈ arcsin (Fe
Fc ), то можно
говорить о продольной устойчивости плазмоида в данных условиях и о реальности существования таких объектов техногенного происхождения, учитывая, что формирование плазмоидов осуществляется при наличии восходящих воздушных потоков.
Уяснив основные механизмы, характеризующие устойчивость плазмоидов, сформулируем систему уравнений, определяющих концентрацию аэроионов, тяжелых ионов, концентрацию капельной влаги, способствующей их образованию, и приведем краткое изложение методики измерений основных параметров, необходимых для определения указанных величин.
Относительная влажность D |
|
|
Таблица 8.1 |
|||||
при заданной температуре окру- |
|
|
||||||
|
|
|
|
|||||
жающего воздуха Т определяется |
Зависимость давления |
|||||||
отношением D = P/Pнас(Т), где P – |
насыщенных паров воды |
|||||||
парциальное давление паров воды |
от температуры Pнас(Т) [27] |
|||||||
в атмосфере |
(искомая величина) |
|
|
|
|
|||
t, °C |
P, кПа |
t, °C |
P, кПa |
|||||
при температуре Т; Pнас – парци- |
||||||||
- 5 |
0,401 |
12 |
1,401 |
|||||
альное давление паров воды при |
|
|||||||
|
- 4 |
0,437 |
13 |
1,497 |
||||
термодинамическом |
равновесии |
|
- 3 |
0,463 |
14 |
1,597 |
||
при заданной температуре (задан- |
|
- 2 |
0,517 |
15 |
1,704 |
|||
ное таблично как функция темпе- |
|
-1 |
0,563 |
16 |
1,817 |
|||
ратуры, табл. 8.1) [27]; D – изме- |
|
0 |
0,611 |
17 |
1,937 |
|||
ряемая величина и P = DPнас. |
|
1 |
0,656 |
18 |
2,062 |
|||
|
2 |
0,705 |
19 |
2,196 |
||||
Используя |
закон |
Авогадро |
|
|||||
|
3 |
0,757 |
20 |
2,337 |
||||
PV = (m/μ)RT, где V – рассматри- |
|
|||||||
|
4 |
0,813 |
21 |
2,248 |
||||
ваемый объем (Vед = 1 см3); т – |
|
5 |
0,872 |
22 |
2,642 |
|||
масса пара (воды); μ – вес грамм- |
|
6 |
0,935 |
23 |
2,809 |
|||
молекулы воды ( μH2O = 18); R – |
|
7 |
1,005 |
24 |
2,984 |
|||
|
8 |
1,072 |
25 |
3,168 |
||||
газовая постоянная, находим m = |
|
|||||||
|
9 |
1,148 |
26 |
3,361 |
||||
= PVед μ/RT или m = DРнасVедμ/RT. |
|
10 |
1,227 |
27 |
3,565 |
|||
Полагая, что |
капля |
воды пред- |
|
11 |
1,312 |
28 |
3,780 |
|
ставляет собой шар радиуса r, величину которого можно найти как усредненное значение радиуса r по спектру капель вида
∞
f (V ) = Aexp(−ν
ν0 ), где A определяют из условия A∫ f (V )dV = Nk ,
0
257
где Nк – количество капель в единице объема, т.е. их концентрация
[23], для |
|
r , |
переходя |
к |
радиусу |
|
капли |
|
V = |
4 |
πr3 |
|
, |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
f (r ) = B exp −(r r0 )3 r2dr , находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
= |
∞ exp |
− |
(r r |
)3 r3dr |
|
∞ exp −(r r |
)3 r2dr . |
|
|
|
|
|||||||||||
|
|
∫ |
|
0 |
|
|
∫ |
|
|
|
|
0 |
|
|
|
|
|
|
|
||||
|
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
Вычисляя интеграл в числителе, находим |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
∞ exp −(r |
|
)3 r3dr = |
1 |
|
|
1 |
|
−4 |
3 |
|
4 |
|
= r4 1 Γ |
|
1 , |
|
|
|
|
||||
r |
|
Γ |
|
|
|
|
|||||||||||||||||
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|||||||||||
∫ |
|
|
0 |
|
|
3 |
|
|
|
0 |
9 |
|
|
|
|
|
|||||||
0 |
|
|
|
r0 |
|
|
|
|
3 |
|
|
|
3 |
|
|
|
|
||||||
где Г(z) – гамма функция и Г(z + 1) = zГ(z). Вычисление в знамена-
теле дает: r03/3, так, что для |
r |
получаем: r = 0,893r . Найденное |
r |
||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
позволяет определить массу воды Vед как |
|
(4 3)πr 3ρNk , |
где ρ = |
||||||||||||
= ρ(T) – плотность воды (табл. 8.2), и, в конечном итоге, вычислить |
|||||||||||||||
|
|
|
Таблица 8.2 |
значение Nk как |
|
|
|
|
|
||||||
|
|
|
|
|
|
N |
|
= |
|
3DPнас (T )μ |
V . |
|
|
|
|
Зависимость плотности воды |
k |
4πr 3ρ(T )RT |
|
|
|
||||||||||
|
|
|
ед |
|
|
|
|||||||||
ρ(Т) от температуры [27] |
При относительной влажности |
||||||||||||||
|
|
|
|
|
|
||||||||||
t, °C |
ρ, г/см3 |
t, °C |
ρ, г/см3 |
D ≈ 50−100 % и t = 20 °С (Т = 293 К) |
|||||||||||
0 |
0,999841 |
16 |
|
0,998943 |
R = 461 Дж/(кг К) для r0 |
= 10-3 см |
|||||||||
1 |
0,999900 |
17 |
|
0,998775 |
значение Nk |
составляет величину |
|||||||||
2 |
0,999941 |
18 |
|
0,998596 |
~105 в 1 см3. Задаваясь значением |
||||||||||
3 |
0,999965 |
19 |
|
0,998406 |
скорости |
конвективного |
подъема |
||||||||
4 |
0,999973 |
20 |
|
0,998205 |
воздушного потока |U0| = 10 |
2 |
см/с |
||||||||
5 |
0,999965 |
21 |
|
0,997994 |
|
||||||||||
6 |
0,999941 |
22 |
|
0,997772 |
и учитывая установившуюся ско- |
||||||||||
7 |
0,999902 |
23 |
|
0,997540 |
рость падения капель относитель- |
||||||||||
8 |
0,999849 |
24 |
|
0,997299 |
но воздуха W = 1 см/с, сформули- |
||||||||||
9 |
0,999782 |
25 |
|
0,997047 |
руем систему уравнений, описы- |
||||||||||
10 |
0,999701 |
26 |
|
0,996785 |
вающих |
|
формирование |
|
рассмат- |
||||||
11 |
0,999606 |
27 |
|
0,996515 |
риваемого ионизированного обра- |
||||||||||
12 |
0,999498 |
28 |
|
0,996235 |
зования при радиоактивном за- |
||||||||||
13 |
0,999377 |
29 |
|
0,995946 |
грязнении подстилающей поверх- |
||||||||||
14 |
0,999244 |
30 |
|
0,995649 |
ности радиусом R0, с мощностью |
||||||||||
15 |
0,999099 |
|
|
|
|
дозы γ-,β-излучения равной H Р/с. |
|||||||||
|
|
|
|
|
|
||||||||||
258

