Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Елохин Автоматизированные системы контроля радиационной обстановки окружаюсчей среды 2012

.pdf
Скачиваний:
65
Добавлен:
12.11.2022
Размер:
4.84 Mб
Скачать

Рис. 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

1k

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ν ≈103 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 = насVедμ/RT.

 

10

1,227

27

3,565

Полагая, что

капля

воды пред-

 

11

1,312

28

3,780

ставляет собой шар радиуса r, величину которого можно найти как усредненное значение радиуса r по спектру капель вида

f (V ) = Aexp(−νν0 ), где A определяют из условия Af (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

Полагая, что период полураспада радиоактивных продуктов загрязнения много больше времени подъема аэроионов, увлекаемых воздушным потоком на высоту 1 км, будем исходить из условий стационарности, при которых можно положить равными нулю члены, содержащие производные по времени от искомых функций. Учитывая, что основную плотность ионизации создают высокоэнергетические электроны, пробег которых в воздухе составляет 1,5−2,0 м, источник ионизирующего излучения представим в виде

цилиндра радиуса R0 и высотой z0: H(r, z) = Н0[η(r) η(r-R0)][η(z) –

η(z z0)]. Систему уравнений, в рамах диффузионно-дрейфовой модели, полагая тяжелые ионы однозарядными , а аэроионы пред-

ставляя как N0+ = N0= N0± = G Hkr , запишем следующим обра-

зом:

kr N + N kr N +Nk−μ+E grad(N + )−μ+ N +div(E)

U0 grad(N + )+ D+div grad(N + ) = 0;

kr N + N kr N Nk+ E grad(N )N div(E) + +U0 grad (N )+ Ddiv grad(N ) = 0;

kcp N + Nks kr Nk+ N kr Nk+ Nk−μW+ E grad (Nk+ )−μW+ Nk+div(E)

U1 grad(Nk+ )+ D+W div grad (Nk+ ) ;

kcn N Nkg kr NkN + kr Nk+NkWE grad(Nk)WNk+div(E) +

+U1 grad(Nk)+ DW div grad(Nk) ;

div

E

)

=

e

N + + N + N = N

;

 

(

 

 

 

 

k

k

 

 

 

 

 

εε0

 

 

 

при условии квазинейтральности:

N + + Nk+ N Nk= 0,

(8.20)

(8.21)

(8.22)

(8.23)

(8.24)

(8.25)

где kr – коэффициент рекомбинации; kcp , kcn – коэффициенты захвата положительных и отрицательных ионов мелкими Nks и крупными Nkg нейтральными каплями ( Nks , Nkg их концентрации); μ-; μ+,

В работе [25] показано, что тяжелые ионы могут быть и многозарядными.

259

μW, μW+ ; D-, D+; DW , D+W – подвижности и коэффициенты диффу-

зии аэроионов и тяжелых ионов (с верхним индексом W) соответственно. Направляя ось Z вверх, имеем: U1 = (U0 W)ez. В цилиндрической геометрии с учетом азимутальной однородности задачи имеем:

 

 

 

 

 

 

E = Ezez

+ Erer ; grad( N )

=

N er + N ez ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

z

 

 

 

 

 

 

 

 

 

divE =

 

 

1

(rE

) +

(E

 

);

div grad( N ) =

1

 

r

 

N

+ N .

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

r r

r

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

r

 

Граничные условия представим следующим образом:

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

=

N

+

 

 

=

N

 

 

=

N

+

 

 

= 0;

(8.26)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

k

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r=0

 

 

 

r=0

 

 

 

 

 

r

 

 

r=0

 

 

r

 

 

 

r=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N (r, z)

 

r→∞

 

= N +

(r, z)

 

 

r

→∞

= Nk(r, z)

 

r→∞

 

= Nk+

(r, z)

 

r→∞

= 0;

(8.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N (r, z)

 

z→∞

 

= N + (r, z)

 

z→∞

= Nk(r, z)

 

z→∞

= Nk+ (r, z)

 

z

→∞

= 0.

(8.28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При радиационном выходе пар ионов на 1 Р поглощенной энергии в воздухе в 1 см3 равном G = 2,08.109 Р-1 см-3 концентрации аэроионов и тяжелых ионов при z = 0 определятся следующими выражениями:

 

(

 

 

) z=0

 

 

(

 

 

) z=0

 

 

 

 

r

 

(

 

)

 

 

 

(

 

0 )

N +

 

r, z

 

 

 

= N

 

r, z

 

 

 

=

G H k

 

η

 

r

 

−η

 

r R

;

N +

(

r, z

 

 

 

= N

(

r, z

 

 

 

= N

 

kc η

(

r

)

−η

(

r

R

.

 

 

 

 

 

 

 

 

 

 

k

 

 

)

z=0

 

k

 

 

)

z=0

 

k

 

 

 

 

 

 

 

 

0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициенты захвата ионов каплями воды имеют вид [25]:

(1+

 

ln (μ+ μ)

 

)

 

 

 

 

kc = 2πDa D± 1+ 2

 

D± (Dac± )

 

 

,

(8.29)

(8.30)

(8.31)

где D±, c± – средние коэффициенты диффузии (0,05 см2/Вс) и тепловая скорость (3,0 104 см/с) легкого аэроиона. При диаметре тя-

желых капель Da 2 10 см-3 kcn 0,63.10-3см3/с, а время захвата τс ~1 (kcn Nkg ) =102 c. Аналогично и для мелких капель при Da ~ 2 10

см-4 kcp 0,63 10-4см3/с.

260

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]