Материал: Мерзликин Г.Я. Основы теории ядерных реакторов. Курс для эксплуатационного персонала АЭС

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Тема 13. Основы кинетики подкритического реактора

241

нов в реакторе должна стабилизироваться на некотором (но не нулевом) уровне. Эту величину средней плотности нейтронов в подкритическом реакторе, которая устанавливается в нём с течением времени, условимся называть устанавливающейся подкри-

тической плотностью

Применительно к подкритическому реактору с источником нейтронов элементарное уравнение кинетики реактора имеет вид:

dn

= δkэ n(t) + s,

(13.2)

dt

l

 

отличаясь от ранее использовавшегося элементарного уравнения кинетики только включённым в правую часть положительным слагаемым - скоростью прибыли нейтронов от источника s, по смыслу представляющей собой величину удельной мощности источника нейтронов.

В подкритическом реакторе величина избыточного коэффициента размножения δkэ отрицательна, что при общей форме записи элементарного уравнения кинетики совершенно незаметно. Чтобы с первого взгляда на уравнение было понятно, что речь идёт именно о подкритическом реакторе, пустимся на маленькую хитрость: положительную величину недостатка величины эффективного коэффициента до единицы в подкритическом реакторе обозначим через δkп и назовём степенью подкритичности реактора:

δkп = 1 - kэ = - δkэ

 

 

(13.3)

С учётом (13.3) уравнение (13.2) будет иметь вид:

 

 

 

 

dn

= - δkп n(t) + s.

 

(13.4)

 

 

 

 

 

 

 

 

dt

 

 

 

l

 

 

 

 

Если вынести сомножитель перед n(t) в правой части за скобки, то:

 

 

dn

= -

δk

п

×

 

n(t) -

sl

 

 

 

 

 

 

 

.

(13.5)

 

 

 

 

 

 

dt

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δkп

 

Когда величина плотности нейтронов в подкритическом реакторе устанавливается

постоянной (то есть n(t) = idem = nу), величина производной dn/dt = 0, и если подставить нулевое её значение в левую часть (13.5), получаем выражение:

0 = - δkп × n

 

+ s, откуда следует: n

 

=

sl

.

(13.6)

 

 

 

l

у

 

у

 

δk

п

 

 

 

 

 

 

 

 

Полученная формула для устанавливающейся подкритической плотности нейтронов подтверждает первоначальные предположения:

а) Величина устанавливающейся плотности нейтронов в подкритическом реакторе тем больше, чем больше величина мощности независимого источника нейтронов.

б) Величина устанавливающейся плотности нейтронов в подкритическом реакторе будет тем большей, чем больше величина среднего времени жизни поколения нейтронов в реакторе ( l ).

в) Величина устанавливающейся плотности нейтронов в подкритическом реакторе будет тем большей, чем меньше величина степени подкритичности реактора.

В числителе выражения (13.6) стоит величина произведения sl, имеющая размерность плотности нейтронов (нейтр./см3) и смысл устанавливающейся плотности нейтронов при величине степени подкритичности δkп = 1. Эту величину называют началь-

ной устанавливающейся плотностью нейтронов в реакторе при данной мощности подкритического источника s, то есть

nу0 = sl

(13.7)

Смысл такого названия можно проиллюстрировать с помощью графика (рис.13.1) зависимости величины устанавливающейся плотности нейтронов в реакторе от величи-

Тема 13. Основы кинетики подкритического реактора

242

ны степени подкритичности δkп (или величины эффективного коэффициента размножения kэ) при различных величинах мощности подкритического источника нейтронов.

nу3 nу2 nу1

0

δkэ 1

nу

nу/nуо=ПКУ

при s3 > s2

при s2 > s1

 

 

 

 

 

 

1

 

 

 

 

при s1

kэ

0

 

 

 

kэ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

0

 

 

δkэ

 

 

1

 

0

 

 

 

Рис. 13.1. Зависимость устанавливающейся подкритической плотности нейтронов от величины степени подкритичности реактора (или эффективного коэффициента размножения) при трёх различных величинах мощности независимого источника нейтронов, и та же зависимость в безразмерном виде.

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

источника нейтронов, называется подкритическим

коэффициентом умножения

(ПКУ):

 

у (δkп )

 

 

 

 

 

ПКУ =

n

=

1

 

.

(13.8)

 

 

 

 

 

 

n

 

δk

 

 

 

уо

 

п

 

 

 

 

 

 

 

Введение величины ПКУ практически удобно тем, что она не зависит от величины мощности источника s и среднего времени жизни поколения нейтронов l, и поэтому безразмерная зависимость устанавливающейся плотности нейтронов от степени подкритичности является общей для всех реакторов (рис.13.1). И смысл самой величины

ПКУ прост: это величина, обратная степени подкритичности реактора.

Итак, самое главное, что следует из графиков рис.13.1, - возрастающий характер зависимости устанавливающейся плотности нейтронов в подкритическом реакторе при уменьшении степени подкритичности (то есть при приближении реактора к критическому состоянию из подкритического).

Из вида этой зависимости становится понятной опасность процедуры пуска: если мы приближаемся к критическому состоянию одинаковыми по величине ступенями изменения степени подкритичности, с каждой последующей ступенькой уменьшения δkп разница устанавливающихся значений плотности нейтронов возрастает всё в большей и большей степени, и значение устанавливающейся плотности нейтронов устремляется к очень большим величинам.

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

Тема 13. Основы кинетики подкритического реактора

243

троны с одинаковым временем жизни, величина устанавливающейся подкритической плотности нейтронов при δkп → 0 теоретически устремляется к бесконечности. В реальном реакторе, где есть мгновенные и запаздывающие нейтроны, величина устанавливающейся плотности нейтронов при достижении реактором критического состояния становится величиной большой, но конечной. То есть эта несуразица с бесконечно большой величиной nу при δkп = 0 - это просто одна из издержек модели первого приближения.

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

Дело в том, что штатные детекторы нейтронов (ионизационные камеры и датчики прямого заряда) в энергетическом реакторе должны быть рассчитаны на большие величины плотностей нейтронов, свойственные реальным режимам работы реактора на мощности, и при этом не должны в результате длительной работы в нейтронном потоке заметным образом изменять свои характеристики. Именно поэтому в энергетических реакторах и используются в качестве штатных детекторов ионизационные камеры, являющиеся довольно грубыми приборами. Использование более чувствительных камер деления позволяет зафиксировать нейтронный поток в процессе пуска реактора значительно раньше, то есть при гораздо более низких значениях, но они не годятся для длительной работы в предусмотренных режимах работы реактора, так как покрытия их внутренних поверхностей тонким слоем высокообогащённого UO2 или PuO2 в больших потоках нейтронов довольно быстро выгорают, из-за чего камеры деления изменяют свои характеристики и поэтому дают искажённые показания величин плотности потока нейтронов.

Частичная “ слепота” пуска технически преодолима двумя путями:

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

использованием в процессе пуска достаточно мощных независимых источников нейтронов, опускаемых в активную зону перед пуском и позволяющим поднять начальную величину подкритической плотности нейтронов в реакторе до уровня, уверенно регистрируемого штатной пусковой аппаратурой СУЗ, или, по крайней мере, существенно уменьшить “ слепой” диапазон степеней подкритичности при пуске.

Нештатная высокочувствительная пусковая аппаратура, естественно, требует места для своего размещения. Поэтому в транспортных реакторных установках, где затеснённость аппаратных выгородок над реакторами не позволяет разместить нештатную пусковую аппаратуру, пользуются вторым из упомянутых методов повышения безопасности пуска, В условиях же энергоблока АЭС, где вопрос размещения нештатной пусковой аппаратуры не представляет серьёзной проблемы, при физических пусках реакторов используется первый метод.

13.3. Переходные процессы при изменениях степени подкритичности реактора

Принципиально нам уже понятно, что переходный процесс в подкритическом реакторе при изменении степени подкритичности реактора от одного значения до другого должен быть процессом перехода величины плотности нейтронов n(t) от одного установившегося значения n1, соответствующего начальной степени подкритичности δkп1 до другого установившегося значения n2, соответствующего другому значению степени

Тема 13. Основы кинетики подкритического реактора

244

подкритичности δkп2. Поэтому единственное, что нас интересует сейчас, это характер этого переходного процесса, то есть ответ на вопрос, какой математической закономерности подчиняется переходный процесс.

Пусть вначале реактор был подкритичен при степени подкритичности δkп1, в результате чего в нём установилась плотность нейтронов nу1 = s l /δkп1. При увеличении степени подкритичности от δkп1 до δkп2, то есть на величину Δδkп = δkп2 - δkп1 (что равносильно уменьшению величины эффективного коэффициента размножения нейтронов на ту же величину: Δδkп = (1 - kэ2) - ( 1 - kэ1) = - (kэ2 - kэ1) = - k) переходный процесс n(t) будет происходить уже при степени подкритичности δkп2 = δkп1 + Δδkп = δkп1 - k, поэтому элементарное уравнение кинетики для этого переходного процесса будет

выглядеть как

 

 

 

 

 

 

 

 

 

 

dn

= - δkп2

 

sl

 

 

dn

= - δkп2 [n(t) - n

у2 ]. .

 

 

n(t) -

,

или

(13.9)

 

 

 

dt

 

dt

l

 

δkп2

 

l

 

 

Уравнение (13.9) - уравнение с разделяющимися переменными:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dn

= - δkп2 × dt.

 

 

 

 

 

 

 

 

 

 

(13.10)

 

 

 

 

 

 

 

 

 

 

 

 

n(t) - nу2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение его следует выполнять при очевидном начальном условии:

 

 

 

 

 

 

 

 

при t = 0

 

 

n(t=0) = nу1

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.11)

 

 

 

 

 

 

Интегрирование даёт следующее:

 

 

 

 

 

 

 

 

 

 

 

ln[n(t) - nу2

] = - δkп2 t,

а подстановка начального условия в это выражение:

 

 

 

n(t) nу2

l

 

 

 

 

 

 

 

 

 

n(t) nу2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δk

п2

 

 

 

 

 

 

 

 

δk

п2

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

= −

 

 

t,

 

или

 

 

 

 

= exp

 

 

t ,

 

откуда:

 

 

n у1

 

nу2

 

l

 

 

nу1 nу2

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(t) = n

 

(n

 

 

n

 

 

 

δk

л2

 

 

 

 

 

 

 

 

nу1

 

 

 

 

δk

п2

 

 

 

 

 

) exp

 

 

t , или n(t) = n

 

 

1 (1

 

 

 

) exp

 

 

t . (13.12)

у2

у2

у1

 

 

 

 

у2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

nу2

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что отношение начальной и конечной плотностей нейтронов

 

 

 

sl

 

 

 

nу1

=

 

δk

п1

 

=

δk

n у2

 

sl

 

δk

 

 

 

 

δkп2

п2

=

δkп1 + Dδkп = 1 +

Dδkп = 1 -

Dk

,

 

 

п1

 

δkп1

δkп1

δkп1

формулу (13.12) можно привести к виду, в котором в правой части фигурируют только исходные данные (nу1) и величина вносимого изменения степени подкритичности (Δδkп) или, лучше, величина изменения эффективного коэффициента размножения k = kэ2 - kэ1, соответствующая этому изменению степени подкритичности ( k = -Δδkп):

 

 

nу1

 

 

 

k

 

 

δk

п1

- Dk

 

 

n(t) =

 

 

 

 

1

+

D

exp

-

 

 

t

 

.

(13.13)

 

 

 

 

 

 

 

 

 

 

 

Dk

 

 

 

δkп1

 

 

 

 

l

 

 

 

 

1 -

 

 

 

 

 

 

 

 

 

 

 

δkп1

 

 

 

 

 

 

 

 

 

 

 

 

 

Наконец, выражение для переходного процесса станет ещё прозрачнее, если выразить величины степени подкритичности δkп1 через величину эффективного коэффициента размножения kэ1:

 

n(t)

=

 

1 - kэ1

 

+

 

Dk

k

э1 + Dk -1

 

 

 

 

 

 

 

1

 

 

exp

 

 

 

t

.

(13.14)

 

nу1

 

 

 

- kэ1

 

 

 

 

 

1 - kэ1 - Dk

1

 

l

 

 

 

Выражение (13.14) свидетельствует о том, что переходный процесс в подкритическом реакторе при изменении величины эффективного коэффициента размножения на

Тема 13. Основы кинетики подкритического реактора

245

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

щий - с уменьшении величины эффективного коэффициента размножения на

k (или

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

k) (см. рис.13.2).

 

 

kэ(t)

 

 

 

 

 

 

kэ(t)

 

 

 

kэ2

 

 

 

kэ1

 

 

 

kэ1

 

k

 

k

kэ2

 

 

0

 

 

 

t

0

 

 

t

 

n(t)

 

 

 

 

 

 

n(t)

 

 

nу2

 

 

 

 

 

 

nу1

 

 

nу1

 

nу2

 

0

 

t

0

t

Рис.13.2. Экспоненциальный характер переходных процессов n(t) в подкритическом реакторе: а) при скачкообразном уменьшении степени подкритичности реактора (увеличении эффективного коэффициента размножения) и б) при увеличении степени подкритичности реактора (уменьшении

kэ).

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

Экспонента, как известно, асимптотическая кривая: она достигает своего теоретически стационарного значения при бесконечно большом времени переходного процесса (при t → ∞). Практически же (с точностью до 1%) любая экспонента подходит к своему стационарному значения за время, равное четырём - пяти её периодов Т (на это обращалось внимание ещё при рассмотрении закона радиоактивного распада).

Период экспоненциального процесса - величина, обратная тому, что стоит в показателе перед переменной t. Следовательно, в данном случае (см. выражение (13.14)):

Т =

 

l

 

 

 

, или T =

l

.

(13.15)

 

 

 

 

 

 

 

kэ1

+

k 1

δkп2

 

Поэтому величина практического времени установления подкритической плотно-

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

 

 

 

t

 

(4 ÷5)l

.

 

 

(13.16)

 

 

у

 

 

 

 

 

 

 

 

δkп2

 

 

 

 

 

 

 

 

 

 

 

 

Время практического установления подкритической плотности нейтронов в реакторе, как и сама величина устанавливающейся плотности нейтронов, определяется величиной степени подкритичности реактора, которая задаётся ему после очередного подъёма органов компенсации реактивности.

Чем ближе реактор к критическому состоянию после очередной ступени подъёма органов компенсации реактивности, тем больше величина практического времени установления подкритической плотности нейтронов в реакторе.

13.5. Процедура ступенчатого пуска и ядерная безопасность реактора

Источник: https://studfile.net/preview/16459748/