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

Семенов - учебное пособие по кинетике и регулированию ЯЭУ

.pdf
Скачиваний:
317
Добавлен:
26.05.2014
Размер:
955.2 Кб
Скачать

Если в реакторе в качестве топлива используются 235U и 239Pu, то накопление шлаков определяется уравнением

dNi

= (P

σf N

+ P

σf N )Ф - N σ

Ф .

(3.25)

 

dt

5i

5 5

9i

9 9

i i

 

 

 

 

 

 

 

 

 

Здесь Рi – выход осколков в процессе деления; σi – сечение поглощения нейтронов. Полагая, что N9<<N5 и что шлаки поглощают только тепловые нейтроны, имеем

dNi

+ σ N Φ = ΦP

σf N

e-σ5Фt .

(3.26)

 

dt

i i

5i

5 5o

 

 

 

 

 

 

 

Умножая левую и правую части уравнения (3.26) на интегрирующий множитель eσiФt , имеем

d

(Nie

σiФt

) = ФP5iσ

f

e

-(σ5

-σi )Фt

.

(3.27)

dt

 

5 N5o

 

 

 

 

 

 

 

 

 

 

 

После интегрирования (3.27) с учетом начального условия Ni(0) = 0, получим

Ni (t) =

P5 N50σ5

(e

-σiФt

e

-σ5

Фt

) .

(3.28)

5

− σi )

 

 

 

 

 

 

 

 

 

 

 

Для расчета шлаков нужно знать нейтронные сечения всех осколков деления, а также продуктов их радиоактивного распада. Поэтому детальный расчет шлаков весьма сложен и громоздок. Обычно шлаки делят на три группы: с высоким, средним и малым сечениями поглощения. Поглощение нейтронов в шлаках по отношению к поглощению в ядерном топливе называют шлакованием и обозначают через q:

qi =

 

 

σiNi

 

 

.

(3.29)

 

 

 

 

 

 

 

σaN + σaN + σaN

 

5

5

8

8

9

9

 

 

При оценочных расчетах все шлаки объединяют в

одну группу с σшл = 50 бн

(барн). Кроме

того, грубо

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

71

 

 

M

c2

= M c2 + Q ,

(3.30)

 

 

т

 

 

i

 

 

 

где Q = Ef

Mт

NA – выделенная

в результате

делений

m

 

 

 

 

 

 

 

 

 

энергия.

 

 

 

 

 

 

 

 

 

Отсюда

M = M

(1 -

EfNA

) .

(3.31)

 

 

 

 

i

т

 

mc

2

 

 

 

 

 

 

 

 

 

 

Поскольку Ef NA << 1 , то количество образовавшихся mc 2

продуктов примерно равно сгоревшему топливу.

Следует заметить, что при глубоком выгорании топлива в реакторе образуется не только плутоний, но и другие трансурановые элементы: америций 242Am, 243Am, кюрий 242Cm, калифорний 250Cf. Калифорний является источником нейтронов с высокой интенсивностью I » 4,34×1012 1/с. Он находит большое применение в промышленности и научных экспериментах.

3.5. Отравление реактора ксеноном

Ксенон является продуктом деления урана со сравнительно большим выходом и большим сечением поглощения нейтронов. Он быстро накапливается и приводит к отравлению реактора. Рождение ксенона происходит по двум каналам:

Гибель ксенона идет по следующей радиоактивной цепочке:

72

Составим уравнение рождения и гибели ксенона. Концентрацию ядер йода обозначим через N1, а концентрацию ядер Xe – N2. Поскольку период полураспада теллура (Te) мал по сравнению с периодом полураспада йода, то можно воспользоваться принципом подчинения, согласно которому долгоживущие продукты

распада подчиняют себе короткоживущие, т.е. dNTe = 0 , и dt

уравнение для теллура писать не нужно. Запишем уравнения рождения и гибели для йода и ксенона:

 

 

dN1

= γ

1

f

Ф - λ N ,

 

 

 

 

(3.32)

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN2

= λ N − λ N

2

− σ N

Ф + γ

2

f

Ф .

(3.33)

 

dt

1

1

 

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь γ1 и γ2 соответственно выход йода и ксенона при делении урана.

Поскольку γ 2 << γ1 , то последним слагаемым в

уравнении (3.33) можно пренебречь. Характер решения системы уравнений (3.32) и (3.33) зависит от начальных условий, поэтому мы рассмотрим несколько типичных ситуаций.

Пуск реактора. При пуске реактора со свежей загрузкой топлива N1 = N2 = 0 при t = 0. Уравнение (3.32)

умножим на eλ1t и запишем в виде

 

 

 

 

d

(N eλ1t ) = γ

1

f

Фeλ1t .

(3.34)

 

 

 

dt

1

 

 

 

 

 

 

 

 

 

 

Так как поток быстро достигает стационарного состояния по сравнению с расходом йода, то Ф(t) можно считать постоянной величиной. Решение уравнения (3.34)

имеет вид

или

λ t

=

 

γ1 f

Φ

λ t

− 1)

N1e 1

 

 

 

 

(e 1

 

λ1

 

 

N (t) =

γ1 f Ф

(1- e−λ1t ) .

 

1

 

 

λ1

 

 

 

 

 

 

 

 

 

 

 

73

При постоянном уровне мощности концентрация йода достигает стационарного состояния:

N1ст =

γ1

f Ф

.

(3.35)

 

λ1

 

 

 

 

 

Для решения уравнения (3.33) воспользуемся тем же

приемом – умножим все слагаемые на интегрирующий

множитель e(λ 2 2Ф)t :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

(N2 e(λ2 2Ф)t ) = Φγ1Σ f (1 − e−λ1t )e(λ2 2Ф)t

. (3.36)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После интегрирования получим

 

 

γ1 f (e2 −λ12Ф)t − 1)

 

 

Ф)t

 

=

γ1 f Ф

 

 

 

Ф)t

− 1)

N 2e

 

2

2

 

 

 

 

 

 

 

(e

 

 

2

2

 

 

 

 

.

 

 

 

 

 

λ 2

 

+ σ2Ф

 

 

 

 

 

 

λ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ σ2Ф - λ1

или

 

 

 

 

 

 

 

 

γ1 f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N2

(t) =

 

Ф

(1 − e (λ 2 2Ф)t )

 

 

 

 

 

 

 

λ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

+ σ 2 Ф

 

 

 

 

 

 

.

(3.37)

 

 

 

 

 

 

 

 

Ф

 

(e −λ1t e (λ 2 + σ2Ф)t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

f

 

 

 

 

 

 

 

 

 

λ 2 + σ 2 Ф - λ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В стационарном состоянии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N2ст =

 

γ1 f Ф

.

 

 

 

 

(3.38)

 

 

 

 

 

 

 

 

 

λ2 + σ2Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует обратить внимание на то, что стационарная концентрация йода прямо пропорциональна потоку, т.е. мощности реактора, тогда как в зависимости N2(Ф) имеется насыщение (рис.3.4).

γ ∑ f

σ2

Рис.3.4. Зависимость стационарной концентрации йода (N1) и ксенона (N2) от плотности нейтронного потока

74

Отравление реактора принято определять отношением

 

 

 

 

 

 

 

 

 

 

P

= Xe

= s2N2 .

 

 

 

 

 

 

(3.39)

 

 

 

 

 

 

 

 

 

 

Xe

U

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частности, в стационарном состоянии

 

 

 

 

 

 

 

 

 

P

 

 

=

 

g f Фs2

 

=

gФs2

 

 

f

.

 

(3.40)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xe

 

(l2 + s2Ф) u

l2 + s2Ф u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь U – сечение захвата нейтронов ядрами урана. При

температуре

 

 

зоны

~700

 

К

s1

= 1,87×106

бн,

l2

=

2,1×10-5

1/с,

 

f

 

= 0,55

для естественного урана и

 

U

 

f

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0,72,

 

= 0,76 для топлива с обогащением 2 % и

 

 

 

 

 

 

U

 

U

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

3 %

соответственно.

Для

чистого

235U

 

=

0,85.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

В табл.3.1 приведены значения стационарного отравления в зависимости от потока нейтронов для естественного урана.

Таблица 3.1. Зависимость отравления реактора от потока нейтронов

Ф, 1/(см2×с)

1012

1013

1014

РХе

0,003

0,017

0,029

В табл.3.2 приведены предельные значения

отравления для больших потоков ( s2Ф >> l2 )

при

различных значениях обогащения топлива:

 

Pпр

=

g f

.

(3.41)

 

Xe

 

U

 

 

 

 

Таблица 3.2. Зависимость предельного отравления от обогащения топлива

Обогащение,

Естеств. уран,

2

3

100 чист 5U

%

0,72

 

 

 

Pпр

0,035

0,045

0,048

0,054

Xe

 

 

 

 

 

 

75

 

 

На рис.3.5 показаны зависимости PXe(t) для разных значений Ф.

Рис.3.5. Зависимость отравления реактора от времени при пуске реактора: 1 - Ф1 = 1014 1/(см2×с); 2 - Ф2 = 3×1012 1/(см2×с);

3 - Ф3 = 1012 1/(см2×с)

Отключение реактора

При отключении реактора исходная система уравнений имеет вид

 

 

 

dN1

 

= −λ N ,

 

 

 

 

 

(3.42)

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN2

 

= λ N − λ N

2

.

 

 

(3.43)

 

 

 

 

 

 

 

 

 

 

dt

 

 

1

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начальными условиями

являются стационарные

концентрации йода и ксенона:

 

 

 

 

 

 

N10

=

γ1 f

Ф

;

 

N20 =

γ ∑f Ф

.

(3.44)

 

λ1

 

 

 

λ2 + σ2Ф

 

 

 

 

 

 

 

 

 

 

 

 

Первое уравнение интегрируется сразу:

 

 

 

 

N (t) = N

e-l1t .

 

 

 

 

(3.45)

 

 

1

 

 

10

 

 

 

 

 

 

 

Интегрирование второго уравнения осуществим с

использованием интегрирующего множителя el2t . Для этого умножим все члены уравнения (3.43) на этот множитель:

d

(N

2

el2t ) = λ N e(l2 -l1)t .

(3.46)

 

dt

1

10

 

 

 

 

 

После интегрирования получим

76

 

 

 

 

 

 

λ

 

t

=

 

 

l1N10

 

(λ

 

 

−λ )t

-1) + C .

 

 

 

 

N2

e

 

 

2

 

 

 

 

 

 

 

 

( e

 

 

2

 

1

 

 

(3.47)

 

 

 

 

 

 

 

l2

- l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При t = 0 и N20 = С решение имеет вид

 

 

 

 

 

 

 

= N

 

 

 

 

 

−λ

 

t

+

 

l1N10

 

 

 

 

−λ

 

t

- e

−λ t

 

 

N

2

 

 

 

e

 

 

2

 

 

 

 

 

 

 

( e

 

2

 

 

1 ) ,

(3.48)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

l1 - l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а зависимость P(t) соответственно

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t) =

f Фs2

 

 

 

ge−λ 2t

+

 

 

g1

 

 

 

 

 

−λ

 

t

- e

−λ t

 

[

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( e

 

 

2

 

1

)] .(3.49)

U

l

2

+ s

2

Ф

l + l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения времени, когда отравление достигает максимума, нужно производную по времени от P(t) приравнять нулю:

 

dP

= -

l2 g e−λ 2t

-

g l2 e−λ 2t

+

 

g l1

 

−λ t

= 0 .(3.50)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

1

 

dt

l2

+ s2Ф

 

l1

- l2

l1

- l2

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(λ

−λ

 

)t

=

 

 

 

 

 

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

e

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(l1

- l2 )(

 

1

 

 

 

+

1

 

 

) l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

 

 

 

 

 

 

l1 - l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ s2Ф

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t =

 

 

 

 

 

 

 

 

 

 

 

 

l1 - l2

.

 

(3.51)

и

 

 

Ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1 - l 2

 

 

 

 

 

 

 

 

 

 

 

l

2

(1

+

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2 + s2 Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

 

 

 

 

данные,

 

 

 

 

 

получим

 

при

Фo = 2 ×1014 1/(см2×с), tmax @ 11 ч

Pmax = 0,51 против 0,029 в

режиме номинальной мощности. Характер зависимости приведен на рис 3.6.

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

77

Рис.3.6. Зависимость отравления ксеноном от времени после останова реактора (с=3 %): 1 – Ф = 1013 1/(см2×с); 2 – Ф = 3×1013 1/(см2×с);

3 – Ф = 1014 1/(см2×с)

Поскольку Т1/2 I < Т1/2 Хе, то реакция лимитируется ксеноном, идет его накопление и концентрация ксенона растет. Когда йода остается мало, рождения ксенона почти нет и идет распад ксенона. В свете сказанного, концентрация ксенона сначала нарастает, а затем будет уменьшаться. Точно так же будут протекать процессы при переходе реактора с высокого уровня энергии на более низкий уровень, только отравление будет меньше, чем при отключении.

3.6. Отравление ксеноном в переходных режимах

При переводе реактора с одного уровня мощности на другой начальными условиями для системы уравнений (3.32) и (3.33) являются стационарные концентрации йода и ксенона в исходном режиме. При нахождении решений будем считать, что поток меняется от Фо до Ф очень быстро по сравнению с периодами полураспада йода и ксенона.

Решение уравнения (3.32) имеет вид

N (t) =

γ ∑f Ф

(1 − e

-l t

) + N e

-l t

 

 

 

1

1 ,

(3.52)

 

 

1

λ1

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

78

 

 

 

 

где N10 = γ ∑f Ф – начальная концентрация йода.

λ1

В уравнение (3.36) также добавляется одно слагаемое:

d

(N2 e(λ2 2Ф)t ) = γ1 f Ф(1 − e−λ1t ) e(λ 2 2Ф)t +

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ N e(λ2 2Ф−λ1 )t

= γ

1

f

Ф e(λ 2 2Ф)t + (N − γ

1

f

Ф) e(λ2 −λ1 2Ф)t .

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.53)

 

 

Интегрирование (3.53) даст следующее решение:

N2 (t) =

γ1 f Ф

 

 

 

 

 

(λ

 

 

+ σ Ф)t

 

+

 

 

 

 

 

 

 

 

 

(1- e

 

 

2

 

2

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ 2 + σ2Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.54)

 

 

N10 − γ ∑ f Ф

 

 

− λ t

 

 

 

(λ

 

+ σ Ф)t

 

 

(λ

 

+ σ Ф)t

+

 

 

 

e

 

) + N20

 

 

 

 

 

(e

 

 

1

 

 

 

2

 

2

 

e

 

 

2

2

.

λ 2 − λ1 + σ2Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответственно определяется и степень отравления

реактора:

 

 

 

 

 

 

 

 

 

σ2N2 (t) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t) =

 

 

 

 

 

 

 

(3.55)

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

Следует обратить внимание на то, что характер решений уравнений (3.32) и (3.33) таков, что вначале при снижении мощности степень отравления реактора растет, а при увеличении мощности уменьшается. Это видно как из полученных решений, так и непосредственно из уравнений (3.32) и (3.33). Увеличение отравления со снижением мощности мы уже объяснили, теперь проанализируем обратный процесс. При увеличении мощности увеличивается Ф, следовательно, и слагаемое σ2N2Ф в уравнении (3.33), т.е. происходит поглощение

нейтронов ксеноном, что и приводит к снижению его концентрации. Процесс подчиняется следующей логической цепочке:

ρ

ρ

79

Отсюда видно, что по ксенону имеет место положительная обратная связь.

3.7. Связь отравления с реактивностью. Йодная яма

Установим связь отравления реактора ксеноном с его реактивностью. Отравление приводит к изменению изотопного состава зоны, что сказывается на коэффициенте размножения:

 

 

= ηε ϕ f

e

B2

τ

 

K

 

 

 

.

(3.56)

эф

 

 

 

 

1

+ B2L2

 

 

 

 

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

реакторов B2L2 << 1. Напомним, что B2 ~

1

, где R –

R2

 

 

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

 

 

=

 

U

 

f

 

 

 

.

 

(3.57)

 

U + ∑З

При отравлении

коэффициент

использования

тепловых нейтронов принимает значение

 

 

 

 

 

U

 

 

f

= U + ∑З + ∑ Xe .

(3.58)

 

Поэтому изменение

коэффициента

размножения

определяется формулой

 

 

 

 

80