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

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

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

Графики зависимости мощности реактора и температуры зоны представлены на рис.2.18 и 2.19. Если

β2

αN

 

 

 

>

o

,

то возмущение будет затухать по

 

2

Cpτ

 

 

апериодическому закону (рис.2.20 и 2.21).

Рис.2.18. Зависимость возмущения мощности реактора от времени

при αN0 = 200 , β2 = 10,24

Cpτ 4τ

Рис.2.19. Зависимость возмущения температуры реактора от времени

при αN0 = 200 , β2 = 10,24

Cpτ 4τ

61

Рис.2.20. Зависимость возмущения мощности реактора от времени

при αN0 = 10 , β2 = 10,24

Cpτ 4τ

Рис.2.21. Зависимость возмущения температуры реактора от времени

αN

0

= 10 ,

β

2

= 10,24

при

Cpτ

 

 

 

 

62

 

 

 

Глава 3. ИЗМЕНЕНИЕ ИЗОТОПНОГО СОСТАВА ТОПЛИВА

В процессе работы реактора происходит непрерывное изменение изотопного состава топлива вследствие его взаимодействия с нейтронами. Топливо выгорает, образуются осколки деления, и изменяется реактивность реактора.

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

Изменение изотопного состава приводит к медленному изменению реактивности, тогда как при изменении температуры реактивность изменяется быстро (практически без запаздывания).

3.1. Взаимодействие нейтронов с топливом

Если в качестве топлива используется уран (смесь 235U и 238U), то имеют место следующие процессы:

1.

2.

Так как плутоний является делящимся нуклидом, то имеет место процесс деления:

3.

4. Радиационный захват

63

Радионуклид 241Pu тоже является делящимся нуклидом:

Как видно из представленных ядерных реакций, в реакторе происходит выгорание топлива 235U, накопление и расходование нового топлива 239Pu и 241Pu. При расчете изотопного состава цепочку превращений обрывают на 236U и 242Pu. Не учитываются распад изотопов урана и плутония, т.к. периоды полураспада их очень велики. В реакторе идет накопление долгоживущих продуктов нейтронных реакций. Продукты образуются как непосредственно в результате делений, так и в результате цепочек распадов. Все вещества, поглощающие нейтроны без деления, являются « вредными» и по мере накопления снижают реактивность реактора.

При использовании ториевого цикла изменение изотопного состава происходит по следующей цепочке превращений:

1.

2.

3.

4.

n + 233 U

осколки деления

 

234

,

 

 

U

5.

Накопление 235U невелико, т.к. выход 234U мал. После остановки реактора еще длительное время происходит накопление 233U, что ведет к увеличению реактивности.

64

3.2. Выгорание ядерного топлива

При составлении дифференциальных уравнений превращений сделаем следующие предположения. Первое из них основывается на том, что реактор работает на тепловых нейтронах и в качестве топлива используется слабо обогащенный уран. При расчете накопления плутония 239Р не учитываются промежуточные продукты 239U и 239Np, т.к. их периоды полураспада малы, а количество 238U остается примерно постоянным. С учетом сказанного, система уравнений имеет вид

 

 

 

dN5

 

 

 

= −σ N Ф ,

 

 

 

 

 

 

 

(3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

5

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN6

 

 

= σγ N

Ф- σ

N

Ф,

 

 

 

(3.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

5

5

 

6

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN9

= σγ N Ф + ε(1- ϕ)(ν N σ

5

+ ν N σ

9

+ ν N σ )Ф - σ N Ф ,

 

dt

8

8

 

5

5

 

 

9

9

 

1

1

1

9

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN o

 

 

 

 

= σ9 N9Ф - σo N o Ф,

 

 

 

(3.4)

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN1

= σ N Ф - σ N Ф .

 

 

 

(3.5)

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

o

 

o

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь σI – полное сечение поглощения нейтронов. Выгорание топлива идет неравномерно по объему

зоны, т.к. плотность нейтронного потока Ф является функцией координат. Кроме того, действует система компенсации и регулирования, и Ф изменяется в пространстве и времени. В таком случае систему уравнений (3.1)–-(3.5) решить можно только численно. Аналитическое решение можно найти только в предположении постоянства потока Ф.

Уравнение (3.1) не связано с другими уравнениями и может быть проинтегрировано сразу:

N5 = N5oexe(−σ5Фt) .

(3.6)

Подставим (3.6) в уравнение (3.2):

65

dN6

+ σ N

Ф = σ N

e-σ5Фt .

(3.7)

 

dt

6

6

5 5o

 

 

 

 

 

 

 

Умножим левую и правую части уравнения (3.7) на

интегрирующий множительeσ6Фt , после чего уравнение приводится к виду

 

d

(No eσ6Фt ) = σ γ 5 N5o e-(σ5 -σ6 )Фt .

(3.8)

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

Интегрирование (3.8) дает

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

N 6 = σ γ 5 N5o e-σ6Фt e −(σ5 −σ6 )Фt dt + C .

(3.9)

 

 

 

 

 

 

 

 

o

 

 

 

 

При t = 0 N6 = 0, С = 0. Вычисляя интеграл, получим

N 6 =

σ γ

5 N5o

(e

-σ

Фt

e

-σ

Фt

) .

(3.10)

 

 

6

 

5

 

σ5

− σ6

 

 

 

 

 

 

 

 

 

 

 

 

 

При решении

уравнения

 

(3.3) учтем,

что

ν5σ5N5 >> ν9σ9N9 ,

ν5σ5N5 >> ν1σ1N1

и, кроме того,

N5 за

время кампании изменяется незначительно, т.е. N5 N5o . Тогда уравнение (3.3) принимает вид

 

 

dN9

 

σ N Ф + ε(1- ϕ)(ν

5

σ

N )Ф - σ N Ф

 

 

 

 

 

dt

8

8

 

 

 

5

5o

 

 

 

9

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

dN9

+ σ N Ф = σ N Ф + ε(1- ϕ)(ν

5

σ N

)Ф . (3.11)

 

 

 

 

dt

9

9

8

8

 

 

 

 

 

5

5o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим левую часть уравнения

 

 

 

 

 

 

 

R = σ8N8Ф + ε(1- ϕ)ν5σ5N5oФ .

 

 

 

(3.12)

 

Здесь R – скорость рождения 239Pu.

 

 

 

 

 

Умножая левую

и

правую

части

 

уравнения на

интегрирующий множитель eσ9N9Фt , имеем

d

(N eσ9Фt ) = Reσ9Фt .

(3.13)

 

dt

9

 

 

 

Интегрируя (3.13), получим

N9

=

R

(1- e

−σ9Фt

) .

(3.14)

σ9Ф

 

 

 

 

 

 

 

66

Здесь учтено, что при t = 0 N9 = 0, т.е. константа интегрирования равна нулю. Аналогичным образом находятся зависимости Nо(t) и N1(t). Для Nо(t) имеем уравнение

 

 

 

 

dNo

+ σ N Ф = σ

ФN (t).

 

 

(3.15)

 

 

 

 

 

 

 

 

 

 

 

dt

o

o

9

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя сюда (3.14) и интегрируя, получим

 

N0

(t) =

К

(1 − e

− σ0Φ t

) +

 

R

 

(e

− σ0Φ t

e

− σ9Φ t

) .(3.16)

σ0Φ

 

 

(σ9 − σ0 )Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Накопление плутония в соответствии с уравнением (3.5) представим в виде

N1(t) = σ0Φe

−σ Φt

t

σ Φt

dt .

(3.17)

1

N0e 1

 

 

0

 

 

 

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

3.3. Воспроизводство ядерного топлива

Из рассмотрения ядерных реакций взаимодействия нейтронов с топливом видно, что в результате деятельности реактора идет накопление и расходование делящихся нуклидов 249Pu, 241Pu (урановый цикл) или 233U (ториевый цикл). Для количественной оценки накопления нового ядерного топлива вводится понятие коэффициента воспроизводства топлива. Далее будем иметь в виду урановый цикл:

dN9

KB =

 

dt

 

,

(3.18)

 

dN

 

 

 

 

 

 

 

 

5

 

 

 

dt

Накоплением 241Pu можно пренебречь, т.к. его количество незначительно. Разделив (3.3) на (3.1), получим

67

 

σγ N

 

 

σ9N9 ) -

σ9N9

 

KB =

8

8

+ ε(1- ϕ)(ν

5

+ ν

. (3.19)

 

 

 

σ N

 

9 σ N

σ N

 

 

5

5

 

 

 

5

5

5

5

 

Из выражения (3.19) видно, что КВ имеет максимальное значение в начале кампании реактора,

когда N9(0) = 0:

 

σγ N

 

 

KB =

8

8

+ ε(1- ϕ)ν5 .

(3.20)

 

σ5N5

 

 

 

Формула (3.20) показывает, что КВ тем больше, чем меньше исходное обогащение топлива, поэтому для накопления плутония выгодно работать на слабообогащенном уране (рис.3.1).

KB

 

 

0,74

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,66

 

 

 

 

 

 

 

 

 

 

N5

 

 

 

 

 

 

 

 

 

 

0,62

 

 

 

 

 

 

 

 

 

 

2,0 2,5 3,0 3,5

 

N8

Рис.3.1. Зависимость коэффициента воспроизводства от обогащения топлива

Второе слагаемое в формуле (3.20) говорит о том, что КВ растет с ростом резонансного поглощения. Для его увеличения можно уменьшить шаг решетки.

Чтобы понять характер зависимости КВ от времени, перепишем формулу (3.19) в виде

 

σγ N

 

σ9N9

 

 

KB =

8

8

+ ε(1- ϕ)ν5

[1 − ε(1 − ϕ)ν9 ] .

(3.21)

 

 

 

σ5N5

 

σ5N5

 

 

Первое

слагаемое

по мере выгорания

топлива

растет по экспоненте eσ5Фt , т.к. N8 за время кампании уменьшается незначительно, а последнее слагаемое

68

растет как eσ5Фt (1 − e−σ9Фt ) , что с учетом знака сначала

приводит к уменьшению КВ, а потом к незначительному росту (рис.3.2).

Рис.3.2. Качественный характер зависимости КВ от времени

Формула (3.20) получена без учета утечки нейтронов. С учетом утечки в процессе замедления она принимает вид

 

σγ N

2τ (ν

 

 

σ9N9 )

σ9N9

 

KB =

8

8

+ ε(1- ϕ)e-B

5

+ ν

. (3.22)

 

 

 

σ N

 

 

9 σ N

σ N

5

 

 

5

5

 

 

 

 

5

5

5

 

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

учетом размножения 238U или 232Th возникает ενэф

нейтронов. Без учета утечки 1 нейтрон пойдет на поддержание цепного процесса, часть будет поглощена продуктами деления и реакторными материалами q, а остальные пойдут на воспроизводство

 

 

 

ενэф = 1 + q + KB ,

(3.23)

где q =

зам

=

1 − f

– характеризует вредное поглощение.

u

 

 

 

f

 

Отсюда

 

 

KB = ενэф (1 + q) .

 

 

 

 

(3.24)

69

Из полученного результата видно, что КВ будет максимальным, когда ενэф – максимально, а q –

минимально. На рис. 3.3 приведены графики зависимости νэф от энергии для основных делящихся нуклидов.

ν

Рис.3.3. Качественная зависимость числа вторичных нейтронов νэф от энергии

Коэффициент ε заметно выше для реакторов на быстрых нейтронах, т.к. 238U и 232Th делятся нейтронами с

Е>1 МэВ. Величина q = 1 − f примерно одинакова для всех f

реакторов. Отсюда следует, что КВ для реакторов на быстрых нейтронах является самым большим и может превосходить 1. В реакторах на тепловых нейтронах предпочтительным является ториевый цикл, в котором KB ³ 1 . В реакторах на быстрых нейтронах предпочтителен уран-плутониевый цикл, для которого КВ~1,5.

3.4. Шлакование

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

70