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

Савандер Физическая теория ядерных реакторов ч.2 2013

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

нения (3.1.1) и (3.1.2) преобразуются в следующую систему уравнений для ядерных концентраций 135 J и 135 Xe

dρJ t

ωJ f Ф λJρJ t ;

 

 

 

 

dt

 

 

 

 

(3.2.1)

 

 

 

 

dρXe t

Xe

 

 

 

 

ωXe f Ф λJρJ t λXeρXe t σa

ρXe t Ф.

dt

 

 

 

 

 

В системе (3.2.1) Ф – средний по объему реактора поток нейтронов, f – среднее по объему реактора сечение деления. Эти вели-

чины определяют мощность реактора: W E f Ф f VА.З. , где E f -

средняя энергия на один акт деления.

В общем случае начальные условия для системы уравнений

ρJ 0 ρ0J

(3.2.1) таковы: , однако, в частном случае, при пер-

ρXe 0 ρ0Xe

воначальном запуске реактора,

когда концентрации ядер 135 J и

135 Xe равны нулю, начальные условия будут однородными

ρ

0 0;

J

.

 

ρXe 0 0.

Очевидно, что стационарное отравление имеет место в случае,

когда ядерные концентрации 135 J и 135 Xe достигают своих равновесных значений и перестают изменяться. Величину этих концентраций можно определить из условия равенства нулю производных в левых частях уравнений системы (3.2.1):

dρJ t

0

 

 

 

dt

 

 

 

 

 

 

 

t

 

dρXe

0

 

dt

 

 

 

 

 

 

ρ

ρ

JωJ f Ф ;

λJ

Xe ωJ ωXe f Ф .

λXe aXeФ

51

135 Xe

Полагая ω0 J Xe , получаем следующие выражения для рав-

новесных концентраций 135 J и 135 Xe

 

 

 

J f Ф

 

 

0 f Ф

 

 

ω

ω

 

ρJ t

 

 

 

, ρXe t

 

.

(3.2.2)

 

 

 

 

 

 

 

λJ

λXe aXeФ

 

Из соотношений (3.2.2) видно, что равновесная концентрация 135 J прямо пропорциональна потоку нейтронов. Этот факт имеет

простую физическую суть: интенсивность образования 135 J (реакция деления) прямо пропорциональна потоку нейтронов, а канал исчезновения из реактора (радиоактивный распад) от потока нейтронов не зависит и определяется только постоянной распада.

Зависимость равновесной концентрации от потока нейтронов в реакторе имеет более сложный характер. Для удобства анализа обозначим Ф0 поток, при котором скорость радиоактивно-

го

распада

135 Xe равна скорости реакции

 

поглощения:

λ

Xe

σXeФ

,

тогда

выражение для равновесной

 

концентрации

 

a

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 f Ф

 

 

 

135 Xe примет следующий вид: ρXe t

ω

 

 

. При потоке

 

λXe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф0

 

нейтронов

Ф

Ф0

равновесная концентрация 135 Xe будет про-

 

 

 

 

 

 

 

 

 

 

 

0

Σf Ф

 

 

 

порциональна

потоку нейтронов: ρXe

t

ω

,

что отражает

 

 

 

 

λXe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

факт преобладающего вклада в исчезновение

135 Xe реакции ра-

диоактивного распада. В этом случае,

 

равновесная концентрация

135Xe ведет себя аналогично равновесной концентрации 135 J . Наоборот, при потоке Ф Ф0 равновесная концентрация 135 Xe

не будет зависеть от уровня потока нейтронов : ρXe t ω0 f , что

σaXe

объясняется прямой зависимостью канала образования и канала исчезновения 135 Xe от потока нейтронов. Следовательно, в высо-

52

копоточных реакторах ядра 135 Xe не успевают испытать радиоактивный распад, поскольку за время жизни относительно распада происходит поглощение нейтрона с исчезновением по каналу радиационного захвата. Описанные зависимости равновесных концен-

траций 135 J и 135 Xe от потока нейтронов представлены на рис. 3.3.

Рис. 3.3. Зависимости равновесных концентраций 135 J и 135 Xe от потока нейтронов

Порядок величины Ф0 можно оценить таким образом: пусть

σaXe 3,5 106 барн, λXe 2,07 10 5 с-1, в результате получим

Ф0 λXe 1013 нейтр/см2·с.

σXea

Влияние стационарного отравления ксеноном на реактив-

ность. Чтобы оценить влияние стационарной концентрации ксенона на реактивность, рассмотрим величину относительного погло-

Xe

щения нейтронов на ксеноне qXe a , где a – полное макросе-

a

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

53

 

ρ

 

 

Xe

 

 

 

 

qXe

 

Xe

a

. Предельное значение эта величина достигает при

 

a

 

 

 

 

 

ФФ0 , и подставляя соответствующее выражение для равновес-

ной

концентрации ксенона

для данного

случая, а

именно

 

 

 

 

 

 

 

0 f

 

 

 

 

 

0 f

 

 

 

 

 

 

 

 

 

 

 

ρ

 

t

 

 

ω

, получим

q

ω

 

 

 

 

 

1

 

,

где

величина

 

 

 

ω

 

Xe

 

 

 

 

σXe

Xe

 

 

 

 

a

0 1

α

f

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α f

 

c

 

зависит от нуклидного состава топлива и лежит в преде-

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

лах 0,2–0,3. Таким образом, поскольку

 

0 6,0 % ,

величина по-

ω

глощения на ксеноне при равновесной концентрации будет не выше 5,0 % суммарного поглощения нейтронов в активной зоне. Покажем, что данная оценка справедлива и для отрицательной реактивности, вносимой равновесной концентрацией ксенона.

В одногрупповом приближении эффективный коэффициент

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

, где K μφθνef ,

 

2

размножения имеет вид

Keff

 

 

 

 

 

 

α0 – гео-

1 α02 L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

метрический параметр, а

L2

D

 

– длина диффузии нейтронов в

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

реакторе.

Коэффициент размножения на быстрых нейтронах μ ,

вероятность избежать резонансного поглощения φ ,

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

диффузии D практически не зависят от наличия ксенона в реакто-

ре. Произведем следующие преобразования:

 

 

 

 

 

 

 

 

 

θν

 

 

ν f f

 

ν f f

 

 

1

 

 

,

L2

 

 

D

 

 

 

1

 

 

L2

1

 

 

,

ef

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

0

1

q0

 

 

 

 

 

0

1

q0

 

 

0

1 q0

 

 

 

 

 

 

 

a

 

 

 

Xe

 

 

 

 

a

 

 

 

 

Xe

 

 

Xe

 

 

где 0

– полное макросечение поглощения нейтронов в неотрав-

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ленной

активной

зоне

(без

ксенона);

q0

 

 

Xe

поглощение

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xe

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

нейтронов ксеноном относительно неотравленной активной зоны. В результате получаем

54

 

 

μφν f f

 

 

 

 

K

0

 

 

 

 

 

Keff0

 

Keff

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

2 2

 

 

 

 

 

0

 

 

 

0

 

 

0a 1 qXe0 1

 

α0 L0

 

 

1 α02 L20

1

 

 

qXe

 

 

1

 

qXe

 

 

 

 

 

 

 

α2 L2

 

 

 

1 q0

 

α2 L2

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

Xe

 

 

 

 

 

0 0

 

 

 

 

0 0

 

Если в начальный момент времени при нулевой концентрации ксенона Kef0 1, то после установления равновесной концентрации ксенона эффективный коэффициент размножения реактора станет

 

 

 

q0

 

 

 

 

равным Keff

1

 

Xe

 

. Таким образом, в упрощенном случае,

 

α2 L2

 

1

 

 

 

 

 

 

0

0

 

α02 L20

1 , получим следую-

для физически большого реактора

щую оценку для изменения коэффициента размножения (а для случая критического реактора и оценку отрицательной реактивности,

вносимой ксеноном): Keff qXe0 , то есть также примерно 5,0 %.

Так как ксенон накапливается в топливе, а не по всей активной зоне, данная оценка является несколько завышенной. В реальных больших тепловых энергетических реакторах величина реактивности, вносимая равновесной концентрацией ксенона, не превышает 2,5–3,0 %. Данная величина реактивности не является препятствием для работы реактора и компенсируется заранее подготовленными средствами органов регулирования реактора. Однако не следует забывать, что стационарное отравление ксеноном снижает эффективность использования ядерного топлива, поскольку в этом случае нейтроны поглощаются бесполезно.

3.3. Нестационарные ксеноновые процессы, эффект «иодной ямы»

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

запуск реактора из неотравленного состояния;

переход с одного уровня мощности на другой;

55

быстрая останов реактора, работавшего на мощности (эффект «иодной ямы»).

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

нормированных концентраций ρ*J t

ρJ t

и

ρ*Xe t

 

ρXe t

си-

ρJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρXe

стема уравнений (3.2.1) преобразуется к виду

 

 

 

 

dρ* t

λJ λJρ*J t ;

 

 

 

 

 

 

 

 

J

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.3.1)

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dρ*Xe

 

 

 

 

Xe

 

 

 

 

Xe

 

 

 

 

 

 

 

ω

 

ω

*

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

λXe

 

 

 

 

 

λXeρJ t

 

λXeρXe t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

ω0

 

 

ω0

 

 

 

 

 

 

 

В системе (3.3.1) λXe λXe σaXeФ . Решением системы дифференциальных уравнений (3.3.1) с начальными условиями общего

ρ*J 0 ρ0J

типа будут следующие функции:

ρ*Xe 0 ρ0Xe

ρ*J t 1 e λJt ρ0J e λJt ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λXet (3.3.2)

 

 

ωJ

1 ρ0J

e

λJt

e

ρ*Xe t 1

e λXet ρ0Xee λXet λXe

 

 

.

 

 

λXe

 

 

 

ω0

λJ

Применим полученные соотношения при рассмотрении нестационарных процессов.

Выход реактора на мощность из неотравленного состояния.

При запуске реактора после длительной остановки, начальные кон-

центрации 135 J и 135 Xe равны нулю. Таким образом, систему уравнений (3.3.1) будем решать для случая нулевых начальных условий

56

* 0 0;J

ρ*Xe 0 0.

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

ных концентраций 135 J

и 135 Xe получаем соотношения

 

ρ*J t 1 eJt

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.3.3)

ρ*

t 1 e λXet λXe

 

 

J

 

e λJt e λXet

.

ω

 

 

 

 

 

 

Xe

 

 

 

 

 

 

λXe λ

 

 

 

 

 

ω0

 

 

 

 

 

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что из полученных решений следует

ρ*J t 1; ρ*Xe t 1,

что физически отражает факт достижения концентрациями 135 J и

135 Xe своих равновесных значений. Время установления равновесных концентраций определяется показателями соответствую-

щих экспонент вида λt ln 2

t

и составляет 3–4 периода полу-

 

 

T1 2

распада рассматриваемых изотопов, то есть порядка 1,5–2 сут. Равновесная концентрация 135 Xe устанавливается позднее, чем равно-

весная концентрация 135 J , что объясняется не только большей величиной периода полураспада, но и наличием цепочки образования

135 Xe , где 135 J является его предшественником.

Графики изменения концентраций 135 J и 135 Xe , а также графики соответствующего изменения мощности приведены на рис. 3.4.

57

Рис. 3.4. Изменение концентраций 135 J и 135 Xe , и мощности реактора при запуске реактора из неотравленного состояния. Нормировка на равновесное состояние при номинальном уровне мощности (100 %)

Данные по ксеноновым процессам, отображенные на рис. 3.4– 3.7, приблизительно соответствуют реактору типа РБМК-1000, физически большому энергетическому уран-графитовому реактору с тепловым спектром нейтронов. Номинальная мощность реактора составляет 3200 МВт.

58

Переход с одного уровня мощности на другой. Постановка за-

дачи выглядит следующим образом: до момента времени t 0 реактор работал на некотором уровне мощности W1 с потоком нейтронов Ф1 . В момент времени t 0 мощность реактора скачкообразно изменяется до значения W2 с потоком нейтронов Ф2 . Необходимо установить, как изменяются при этом концентрации

135J и 135 Xe .

Очевидно, что в начальный момент времени концентрации 135 J

и135 Xe равнялись равновесным концентрациям для уровня мощ-

ности W1 :

 

 

 

 

 

 

 

 

 

 

Σ Ф

 

 

 

 

 

0

ω

J

 

 

 

 

ρ*J

 

 

 

 

 

f 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0 f Ф1

 

ρ*

 

0

 

 

 

 

.

 

 

 

 

 

 

Xe

 

 

Xe

 

 

 

 

 

λXe

 

Ф1

 

 

 

 

 

 

 

 

 

σa

 

Однако решение системы уравнений (3.3.1) будет отвечать уровню мощности W2 , поэтому необходимо отнормировать начальные условия относительно равновесных концентраций, соответствующих потоку нейтронов Ф2 :

 

ρ*

0

ωJ f Ф1

 

ωJ f

Ф2

 

 

Ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

λJ

 

 

λJ

 

 

 

 

Ф2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0 f Ф1

 

ω0 f Ф2

 

 

λXe σaXe Ф2

 

Ф1

 

λXe2Ф1

 

ρ*

0

 

 

 

 

.

 

Xe

 

 

 

 

 

Xe

 

 

 

 

Xe

 

 

Xe

 

Ф2

 

λXe1Ф2

 

 

 

 

λXe σa Ф1

λXe σa Ф2

λXe σa Ф1

 

 

 

Итак, переходной процесс начинается со стационарных концен-

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

соответствующих потоку Ф1 . Затем кон-

центрации 135 J

и 135 Xe должны выйти на новые стационарные

значения,

соответствующие потоку

Ф2 . Поэтому главный вопрос

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

59

процесса. Для этого рассмотрим значение производной по времени от концентрации ксенона в времени t 0 :

 

 

 

 

 

dρ*Xe 0

 

 

 

 

 

 

Xe

λXe2

 

 

 

 

J

λXe2ρ*J 0 λXe2ρ*Xe 0

 

 

 

 

 

 

ω

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

ω0

 

 

 

 

 

ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xe

λXe2

 

 

 

 

 

J

 

λXe2

Ф1

λXe2

 

λXe2

 

 

Ф1

 

 

 

 

 

 

 

ω

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

0

 

 

 

 

 

 

 

 

 

 

ω

0

 

 

 

Ф2

 

 

 

 

 

 

 

 

λXe1

 

 

Ф2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф Ф

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

ω

Xe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λXe2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф

 

 

 

 

 

 

 

 

 

 

1

 

λXe

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xe

Ф1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

Поскольку

ωXe

ω0 ,

то

их

 

 

 

 

отношение

1 ,

а величина

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

λXe

 

0,5

 

при потоке нейтронов порядка 5∙1013

нейтр/см2 с, от-

 

 

 

 

 

σXeФ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

куда следует, что знак производной

 

dρ*Xe 0

 

 

 

 

в начальный момент

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

времени определяется разностью

Ф1 Ф2 .

 

Таким образом, если

происходит

 

повышение

мощности

 

 

реактора

и

Ф1 Ф2 , то

 

dρ*Xe

0

0

и в начальный момент времени концентрация ксенона

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

будет убывать. Наоборот, если происходит снижение мощности

реактора и

Ф1

Ф2

, то

dρ*Xe

0

0

и в начальный момент времени

dt

 

 

 

 

 

 

 

 

концентрация ксенона будет увеличиваться.

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

60

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