Статьи Соросовского Образовательного журнала в текстовом формате
Статья посвящена физике подкритического реактора, работающего как усилитель энергии ускоренных протонов или более тяжелых ядер. Этот вариант в настоящее время рассматривается как основа безопасной ядерной энергетики. Фундаментальные исследования элементарных процессов с участием релятивистских ядер очень важны для оптимального выбора ускоренного пучка.
ФИЗИКА ПОДКРИТИЧЕСКОГО ЯДЕРНОГО РЕАКТОРАБ. А. АРБУЗОВ
Московский государственный университет
им. М.В. Ломоносова
ВВЕДЕНИЕ
Основным прикладным результатом фундаментальных исследований в ядерной физике явилось становление атомной энергетики. Производимая в ядерных реакторах энергия составляет около 6% всего мирового производства энергии. В некоторых странах (Франция, Швеция) атомные электростанции дают более половины всей электроэнергии. Однако развитие атомной энергетики породило и общественные проблемы, которые наиболее ярко проявились в трагической Чернобыльской катастрофе. После Чернобыля опасность для здоровья людей и окружающей среды, связанная с ядерной энергетикой, вызвала обоснованную негативную реакцию общественного мнения. Возникшие при этом вопросы относились не только к промышленникам и политикам, но и к научному сообществу физиков, работающих в области ядерной физики и физики элементарных частиц. В конце концов выяснилось, что физики разработали ядерный реактор, который, как оказалось, может выйти из-под контроля. Поэтому задача развития безопасной ядерной энергетики, проведение фундаментальных исследований по этой тематике в последние годы привлекают повышенное внимание.
Уже давно известен возможный вариант безопасной ядерной энергетики - освоение управляемого термоядерного синтеза. Однако, несмотря на принципиальную осуществимость этой программы, до сих пор перед исследователями стоят еще не преодоленные технологические трудности. Для завершения программы исследований по управляемому термоядерному синтезу необходимы большие материальные вложения и значительное время. В то же время также достаточно давно известен и другой вариант безопасной энергетики, основанный на работе ядерного реактора в подкритическом режиме, для чего требуется облучение реактора потоком нейтронов. Эти нейтроны могут быть получены с помощью интенсивных пучков протонов или более тяжелых ядер. В последние годы работа в этом направлении значительно активизировалась как в области фундаментальных исследований, так и в разработке конкретных проектов установок, производящих энергию.
В настоящей статье описываются физические основы работы подкритических реакторов и представлены некоторые направления фундаментальных исследований, связанных с этой проблемой.
АТОМНЫЙ РЕАКТОР
Источником энергии реактора служит процесс деления тяжелых ядер (см., например, [1]). Напомним, что ядра состоят из нуклонов, то есть протонов и нейтронов. При этом количество протонов Z определяет заряд ядра Ze: оно равно номеру элемента в Периодической системе Менделеева, а атомный вес ядра A - суммарному количеству протонов и нейтронов. Ядра, имеющие одинаковое число протонов, но различное число нейтронов, являются различными изотопами одного и того же элемента и обозначаются символом элемента с атомным весом слева вверху. Например, существуют следующие изотопы урана: 238U, 235U, 233U, _
Масса ядра M не просто равна сумме масс составляющих его протонов и нейтронов, а меньше ее на величину DM, определяющую энергию связи DE = D(A )Ac 2 (в соответствии с соотношением DE = DMc 2)
M = Zmp + (A - Z )mn - D(A )A,
где D(A )c 2 - энергия связи, приходящаяся на один нуклон. Величина D(A ) зависит от деталей строения соответствующего ядра. Однако наблюдается общая тенденция зависимости ее от атомного веса. А именно, пренебрегая мелкими деталями, можно описать эту зависимость плавной кривой, возрастающей при малых A, достигающей максимума в середине таблицы Менделеева и убывающей после максимума к большим значениям A. Представим себе, что тяжелое ядро с атомным весом A и массой M разделилось на два ядра A1 и A2 с массами соответственно M1 и M2 , причем A1 + A2 равно A либо несколько меньше его, так как в процессе деления могут вылететь несколько нейтронов. Возьмем для наглядности случай A1 + A2 = A. Рассмотрим величину разности масс начального ядра и двух конечных ядер, причем будем считать, что A1 © A2 , так что D(A1) = D(A2),
DM = M - M1 - M2 = - D(A )A + D(A1)(A1 + A2) =
= A(D(A1) - D(A )).
Если A соответствует тяжелому ядру в конце Периодической системы, то A1 находится в середине и имеет максимальное значение D(A1). Значит, DM > 0 и, следовательно, в процессе деления выделяется энергия
Eд = DMc 2.
Для тяжелых ядер, например для ядер урана, (D(A1) - - D(A ))c 2 © 1 МэВ. Так что при A © 200 имеем оценку Eд © 200 МэВ. Напомним, что электрон-вольт (эВ) - внесистемная единица энергии, равная энергии, приобретаемой элементарным зарядом под действием разности потенциалов 1 В (1 эВ = 1,6 " 10-19 Дж). Например, средняя энергия, выделяемая при делении ядра 235U:
Eд = 180 МэВ = 180 " 106 эВ.
Таким образом, тяжелые ядра являются потенциальными источниками энергии. Однако самопроизвольное деление ядер происходит исключительно редко и практического значения не имеет. Если же в тяжелое ядро попадает нейтрон, то процесс деления может резко убыстриться. Это явление происходит с различной интенсивностью для различных ядер, и мерой его служит эффективное поперечное сечение процесса. Напомним, как определяются эффективные сечения и как они связаны с вероятностями тех или иных процессов. Представим себе пучок частиц (например, нейтронов), падающих на мишень, состоящую из определенных объектов, скажем ядер. Пусть N0 - число нейтронов в пучке, n - плотность ядер, приходящаяся на единицу объема (1 см3). Пусть нас интересуют события определенного сорта, например деление ядер мишени. Тогда число таких событий N будет определяться формулой
N = N0nlsэф ,
где l - длина мишени и sэф называется эффективным поперечным сечением процесса деления (или любого другого процесса) при заданной энергии E, соответствующей энергии налетающих нейтронов. Как видно из формулы (5), эффективное сечение имеет размерность площади (см2). Оно имеет вполне понятный геометрический смысл: это площадка, при попадании в которую происходит интересующий нас процесс. Очевидно, если сечение большое, процесс идет интенсивно, а маленькое сечение соответствует малой вероятности попадания в эту площадку, следовательно, в этом случае процесс происходит редко.
Итак, пусть для некоторого ядра мы имеем достаточно большое эффективное сечение процесса деления. При этом, как мы уже отмечали, при делении наряду с двумя большими осколками A1 и A2 могут вылететь несколько нейтронов. Среднее число этих дополнительных нейтронов называется коэффициентом размножения и обозначается символом k? . Тогда реакция идет по схеме
n + A A1 + A2 + k?n.
Родившиеся в процессе (6) нейтроны, в свою очередь, реагируют с ядрами A, что дает новые реакции деления и новое, еще большее число нейтронов. Если k? > 1, такой цепной процесс происходит с нарастающей интенсивностью и, вообще говоря, приводит к взрыву с выделением огромного количества энергии. Однако процесс этот можно контролировать. Во-первых, не все нейтроны обязательно попадут в ядро A: они могут выйти наружу через внешнюю границу реактора, могут поглотиться в веществах, которые специально вводятся в реактор. Таким образом, величину k? , характеризующую процесс (6), можно уменьшить до некоторой kэф , которая равна единице или лишь незначительно ее превышает. Тогда можно успевать отводить производимую энергию и работа реактора становится устойчивой. Тем не менее в этом случае реактор работает в критическом режиме. Неполадки с отводом энергии привели бы к нарастающей цепной реакции и катастрофе. Разумеется, во всех действующих системах предусмотрены меры безопасности, направленные на предотвращение этой возможности, однако аварии, пусть с очень малой вероятностью, могут происходить и, к сожалению, происходят.
Как выбирается рабочее вещество для атомного реактора? Необходимо, чтобы в топливных элементах присутствовали ядра изотопа с большим эффективным сечением деления. В табл. 1 приведены значения этого сечения для некоторых изотопов при реакциях с достаточно быстрыми нейтронами (с энергиями порядка мегаэлектрон-вольт, которые характерны для условий, обсуждаемых ниже). Единица измерения сечения 1 барн = 10- 24 см2. Мы видим две группы значений сечений: большие (233U, 235U, 239Pu) и малые (232Th, 238U). Для того чтобы представить себе разницу, вычислим, какое расстояние должен пролететь нейтрон, чтобы произошло событие деления. Воспользуемся для этого формулой (5). Для N = N0 = 1 имеем
Здесь n - плотность ядер,
где r - обычная плотность и m = 1,66 " 10- 24 г - атомная единица массы. Для урана и тория n = = 4,8 " 1022 см- 3. Тогда для 235U имеем l = 10 см, а для 232Th l = 35 м. Таким образом, для реального осуществления процесса деления следует использовать такие изотопы, как 233U, 235U, 239Pu. Изотоп 235U в небольшом количестве содержится в природном уране, состоящем в основном из 238U, поэтому в качестве ядерного топлива обычно используют уран, обогащенный изотопом 235U. При этом в процессе работы реактора вырабатывается значительное количество еще одного расщепляющегося изотопа - 239Pu. Плутоний получается в результате цепочки реакций
238U + n (g)239U (b)239Np (b)239Pu,
где g означает излучение фотона, а b - b-распад по схеме
Z (Z + 1) + e + .
Здесь Z определяет заряд ядра, так что при распаде (9) происходит переход к следующему элементу таблицы Менделеева с тем же A, e - электрон и - электронное антинейтрино.
Необходимо отметить также, что изотопы A1 и A2 , получающиеся в процессе деления, как правило, являются радиоактивными с временами полураспада от года до сотен тысяч лет, так что отходы атомных электростанций, представляющие собой выгоревшее топливо, очень опасны и требуют специальных мер для хранения. Здесь возникает проблема геологического хранения, которое должно обеспечивать надежность на миллионы лет вперед.
Таким образом, несмотря на очевидную пользу современной атомной энергетики, основанной на работе ядерных реакторов в критическом режиме, она имеет и серьезные недостатки. Это, во-первых, риск аварий, аналогичных Чернобыльской, и, во-вторых, проблема радиоактивных отходов.
Предложение использовать для атомной энергетики реакторы, работающие в подкритическом режиме, полностью разрешает первую проблему и в значительной степени облегчает решение второй.
ЯДЕРНЫЙ РЕАКТОР В ПОДКРИТИЧЕСКОМ РЕЖИМЕ КАК УСИЛИТЕЛЬ ЭНЕРГИИ
Представим себе, что мы собрали атомный реактор, имеющий эффективный коэффициент размножения нейтронов kэф немного меньше единицы. Облучим это устройство постоянным внешним потоком нейтронов N0 . Тогда каждый нейтрон (за вычетом вылетевших наружу и поглощенных, что учтено в kэф) вызовет деление, которое даст дополнительный поток N0kэф . Каждый нейтрон из этого числа снова произведет в среднем kэф нейтронов, что даст дополнительный поток и т.д. Таким образом, суммарный поток нейтронов, дающих процессы деления, оказывается равным
Если kэф $ 1, ряд в (10) расходится, что и является отражением критического поведения процесса в этом случае. Если же kэф < 1, ряд благополучно сходится и по формуле суммы геометрической прогрессии мы имеем
Выделение энергии в единицу времени (мощность) тогда определяется выделением энергии (4) в процессе деления
где k < 1 - коэффициент, равный отношению числа нейтронов, вызвавших деление, к полному их числу. Этот коэффициент зависит от конструкции установки, используемых материалов и т.д. Он надежно вычисляется. В примерах, которые мы будем рассматривать,
k = 0,6.
Осталось выяснить, как можно получить первоначальный поток нейтронов N0 . Для этого можно использовать ускоритель, дающий достаточно интенсивный поток протонов или других частиц, которые, реагируя с некоторой мишенью, порождают большое количество нейтронов. Действительно, например, при столкновении с массивной свинцовой мишенью каждый протон, ускоренный до энергии 1 ГэВ (109 эВ), производит в результате развития ядерного каскада в среднем = 22 нейтрона. Энергии их составляют несколько мегаэлектрон-вольт, что как раз соответствует работе реактора на быстрых нейтронах. Удобно представить поток нейтронов через ток ускорителя
где e - заряд протона, равный элементарному электрическому заряду.
Когда мы выражаем энергию в электрон-вольтах, как в (4), это значит, что мы берем представление E = eV, где V - соответствующий этой энергии потенциал, содержащий столько вольт, сколько электрон-вольт содержит энергия. Это значит, что с учетом (14) можно переписать (12) в виде
Наконец удобно представить мощность установки в виде
где V - потенциал, соответствующий энергии ускорителя, так что VI по известной формуле есть мощность пучка ускорителя: P0 = VI, а R0 в (15) есть коэффициент усиления этой мощности. Оценим этот коэффициент для kэф = 0,98, что обеспечивает надежный запас подкритичности. Все остальные величины известны, и для энергии протонного ускорителя 1 ГэВ имеем
Мы получили коэффициент усиления 120, что, разумеется, очень хорошо. Однако коэффициент (16) соответствует идеальному случаю, когда полностью отсутствуют потери энергии и в ускорителе, и при производстве электроэнергии. Для получения реального коэффициента нужно умножить (16) на эффективность ускорителя rу и КПД тепловой электростанции rэ . Тогда R = rуrэR0 . Эффективность ускорения может быть достаточно высокой, например в реальном проекте сильноточного циклотрона на энергию 1 ГэВ rу = 0,43. Эффективность производства электроэнергии может составлять 0,42. Окончательно реальный коэффициент усиления
R = rуrэR0 = 21,8,
что по-прежнему вполне хорошо, потому что всего 4,6% производимой установкой энергии нужно возвращать для поддержания работы ускорителя. При этом отметим еще раз, что реактор работает только при включенном ускорителе и никакой опасности неконтролируемой цепной реакции не существует.
ВОСПРОИЗВОДСТВО ТОПЛИВА
Для производства энергии в подкритическом режиме требуется хорошо делящийся изотоп. Обычно рассматриваются три возможности: 239Pu, 235U, 233U. Очень интересным оказывается последний вариант, связанный с 233U. Этот изотоп может воспроизводиться в реакторе при облучении интенсивным потоком нейтронов, а это и есть непременное условие работы реактора в подкритическом режиме. Действительно, представим себе, что реактор заполнен смесью природного тория 232Th и 233U. Тогда при облучении реактора нейтронами, полученными с помощью ускорителя, как описано в предыдущем разделе, идут два основных процесса: во-первых, при попадании нейтрона в 233U происходит деление, которое и является источником энергии, и, во-вторых, при захвате нейтрона ядром 232Th идет цепочка реакций
232Th + n (g)233Th (b)233Pa (b)233U.
Каждая реакция деления приводит к убыли одного ядра 233U, а каждая реакция (18) приводит к появлению такого ядра. Если сравниваются вероятности процесса деления и процесса (18), то количество 233U при работе реактора остается постоянным, то есть топливо воспроизводится автоматически. Вероятности процессов определяются их эффективными сечениями согласно формуле (5). Из этой формулы мы получаем условие стабильной работы реактора с постоянным содержанием 233U:
n(232Th)s(232Th) = n(233U)s(233U),
где n( " ) - плотность ядер соответствующего изотопа. Сечение деления s(233U) = 2,784 барн приведено выше, а сечение захвата нейтрона торием при тех же энергиях s(232Th) = 0,387 барн. Отсюда получаем отношение концентраций 233U и 232Th
Таким образом, если мы в качестве рабочего вещества выберем смесь из 88% природного тория и 12% изотопа 233U, то такой состав будет длительное время сохраняться при работе реактора. Положение изменится после того, как будет выработано достаточно большое количество тория. После этого нужно производить смену рабочего вещества, но 233U следует выделить из отработанного вещества и использовать в следующей загрузке. Оценим время, которое может проработать реактор при одной загрузке. Возьмем в качестве примера параметры установки, предлагаемые группой проф. К. Руббиа (сотрудничество нескольких европейских центров). Здесь ток ускорителя 12,5 мА при энергии 1 ГэВ и исходная масса топлива 28,41 т. Топливо состоит из окислов ThO2 и 233UO2 . Исходное количество ядер 232Th 5,58 " 1028. При приведенном значении тока производится 1,72 " 1018 нейтронов в секунду. В силу соотношений (19) и (5) половина нейтронов захватывается торием, это соответствует 2,7 " 1025 захватов в год. Отсюда делается заключение, что при времени работы на одной загрузке порядка нескольких лет будет выработано менее 1% всего количества тория. В проекте принята периодичность замены топлива 5 лет.
Необходимо отметить, что продукты деления 233U, представляющие большую радиационную опасность, с большой вероятностью участвуют в реакциях с нейтронами, в результате которых наиболее опасные продукты деления со средним временем жизни пережигаются, то есть либо переходят в устойчивые изотопы, либо, наоборот, в очень нестабильные, которые быстро распадаются. Таким образом, отпадает необходимость геологического хранения отходов работы атомной электростанции. Это еще одно несомненное преимущество подкритического режима работы ядерного реактора. При этом, разумеется, часть потока нейтронов расходуется на пережигание отходов, что несколько понижает коэффициент усиления (17). Однако эти затраты, вне всякого сомнения, оправданны.
О ВЫБОРЕ СОРТА ЧАСТИЦ В УСКОРИТЕЛЕ
В проекте, разрабатываемом группой К. Руббиа, а также в ряде других проектов для получения пучка нейтронов предлагается использовать ускоритель протонов. Действительно, технология сооружения сильноточных ускорителей протонов хорошо разработана, изучены процессы рождения нейтронов при взаимодействии пучка протонов с массивными мишенями.
Однако отметим, что в последние годы развиваются исследования с использованием пучков более тяжелых ядер высоких энергий, в том числе и в применении к проблеме создания интенсивных пучков нейтронов. В этом случае при столкновении ускоренного ядра с ядром мишени рождается некоторое количество нейтронов и ядерные фрагменты, которые, будучи достаточно энергичными, сами вступают в реакции, порождающие нейтроны и новые ядерные фрагменты, вновь вступающие в реакции, и т.д. Такой процесс называется ядерным каскадом. В результате развития ядерного каскада рождается значительное число нейтронов. Проблема заключается в выборе исходной частицы, дающей максимальное число нейтронов на единицу затраченной на ее ускорение энергии.
Для анализа процессов, вызываемых ускоренными ядрами, удобно ввести удельную энергию, то есть энергию, приходящуюся на один нуклон. Это величина E * = E / A. В первом приближении ядро, летящее в пучке с энергией E, можно рассматривать как совокупность A нуклонов с энергией E * каждый. Тогда действие пучка ядер представляется эквивалентным действию пучка протонов, в A раз более интенсивного и в A раз менее энергичного, что даст то же число нейтронов на единицу затраченной на ускорение энергии (при этом ускорение ядер - процесс технологически несколько более сложный, чем ускорение протонов).
Однако этот вывод справедлив лишь в первом приближении. Величина в выражении (14), вообще говоря, является функцией двух переменных: E и A, а не только их отношения A *. С одной стороны, эту зависимость можно рассчитать из теоретической модели, а с другой - изучить на опыте. Теоретический расчет дает максимальное число нейтронов на единицу затраченной энергии для пучка дейтронов 2H, а далее с ростом A эффективность ядерного пучка медленно убывает. В эксперименте проявился неожиданный эффект. Эти экспериментальные результаты были получены двумя группами физиков в опытах на синхрофазотроне Объединенного института ядерных исследований в Дубне (Россия), который в последние годы работает в режиме ускорения пучков ядер. Одна группа представляла физиков ОИЯИ, другая объединяла в рамках сотрудничества физиков из ОИЯИ, Германии (Марбург), Франции (Страсбур), Греции (Салоники). Обе группы получили согласующиеся между собой результаты: измеренный поток нейтронов, порожденный пучком ядер 12C с полной энергией 44 ГэВ (E * = = 3,65 ГэВ), в полтора раза превышает расчетный, теоретический (см., например, работы [2]). При этом отклонение результатов наблюдений от расчетных предсказаний начинается при достаточно большом значении энергии E, превышающей, согласно данным второй группы, энергию 22 ГэВ.
С большой степенью вероятности причиной такого рассогласования можно считать коллективные эффекты в ядрах. Дело в том, что при столкновении двух ядер наряду с взаимодействием отдельных составляющих их нуклонов между собой может происходить обмен энергией между взаимодействующими ядрами как целыми, то есть в игру вступают сразу все 44 ГэВ, запасенные ядром 12C. В результате образуется сильно возбужденное ядерное состояние, дающее при развале большое количество также возбужденных ядерных фрагментов. По-видимому, эти процессы с заметной интенсивностью происходят при энергиях E порядка 40 ГэВ и более. Например, для ядер аргона 40Ar это происходит уже при удельной энергии E * = 1 ГэВ.
Элементарные процессы, приводящие к увеличению интенсивности ядерных каскадов, составляют предмет исследования автора статьи. Для коллективных эффектов в ядре важным является действие вязкости ядерной материи, что приводит к эффективному трению при движении частиц в ядре. Трение приводит к тому, что область взаимодействия налетающего ядра с ядром мишени как бы расширяется. Вследствие этого увеличивается вероятность вылета возбужденных ядерных фрагментов, что ведет к увеличению выхода нейтронов. Справедливость такой интерпретации составляет предмет теоретических и экспериментальных исследований.
Изучение этой проблемы даст возможность выбрать оптимальный пучок для поддержания работы подкритического реактора. В самом деле, усиление ядерных каскадов при реакциях тяжелых ядер с достаточно высокой энергией может привести к выводу о преимуществе использования тяжелых ядер вместо протонов для работы установок, которые мы описали выше. Таким образом, технический, казалось бы, вопрос о выборе пучка для генерации потока нейтронов оказывается связанным с фундаментальными проблемами физики ядра и элементарных частиц. В этой статье мы постарались подчеркнуть важность фундаментальных исследований, в том числе и для таких жизненно важных проблем, как энергетика.
ЛИТЕРАТУРА
1. Вальтер А.К., Залюбовский И.И. Ядерная физика. Харьков: Основа, 1991.
2. Воронько В.А. и др. // Атом. энергия. 1990. Т. 68. С. 449; 1991. Т. 71. С. 563.
* * *
Борис Андреевич Арбузов, доктор физико-математических наук, профессор физического факультета Московского государственного университета им. М.В. Ломоносова, главный научный сотрудник Института физики высоких энергий (Протвино Московской обл.). Специалист по теоретической физике элементарных частиц. Область научных исследований: квантовая теория поля, взаимодействия частиц и ядер при высоких энергиях. Автор курса лекций по элементарным частицам и более чем 120 статей по квантовой теории поля, физике частиц и ядер.