Как рассчитывают последствия расплавления активной зоны

Как рассчитывают последствия расплавления активной зоны: методология и физика тяжелой аварии

Расплавление активной зоны ядерного реактора (core meltdown) — это наиболее тяжелый сценарий аварии в атомной энергетике. В отличие от проектных аварий, последствия такого события не ограничиваются лишь повреждением топлива. Речь идет о последовательном разрушении тепловыделяющих сборок (ТВС), плавлении топливной матрицы, взаимодействии расплава с бетоном или водой и, в худшем случае, о выходе радиоактивных продуктов за пределы гермооболочки.

Расчет последствий такого сценария — это сложная многофизическая задача. Она требует моделирования нейтронно-физических, теплогидравлических, механических и химических процессов. Ниже представлен систематический обзор того, как научное сообщество и регуляторы подходят к оценке и количественному определению последствий расплавления активной зоны.

Фазы тяжелой аварии как основа для расчета

Любой расчет последствий начинается с идентификации фаз аварии. Современная методология (основанная на документах МАГАТЭ, NUREG-1150 и SOAR) выделяет следующие стадии:

Иллюстрация к статье: Как рассчитывают последствия расплавления активной зоны

Начальная фаза (потеря теплоносителя). Давление в реакторе падает, либо падает уровень воды. Температура оболочек твэлов (цирконий) начинает расти.
Фаза паро-циркониевой реакции. При температурах выше 1000 °C цирконий вступает в экзотермическую реакцию с водяным паром. Выделяется водород и дополнительное тепло. Эта стадия катастрофически ускоряет прогрев.
Фаза расплавления и образования кориума. При температурах около 2000-2850 °C топливный материал (UO₂) и цирконий плавятся. Образуется кориум — расплав, состоящий из оксидов урана, циркония, продуктов деления и конструкционных материалов.
Фаза перемещения расплава. Расплав стекает в нижнюю часть корпуса реактора (днище). Возможно взаимодействие с остатками воды (паровой взрыв) или с сухим бетоном.
Фаза выхода из корпуса (Vessel Failure). Если расплав пробивает корпус, начинается взаимодействие кориума с бетоном фундамента (MCCI — Molten Core-Concrete Interaction) внутри гермооболочки.

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

Теплогидравлический расчет и определение источника тепла

Первый шаг — это определение остаточного тепловыделения. Даже после остановки цепной реакции деления короткоживущие продукты деления (ПД) продолжают выделять энергию. Сразу после остановки реактора остаточное тепловыделение составляет около 6-7% от номинальной мощности. Через час — около 1%, через сутки — менее 0.5%.

Для расчета последствий расплавления используется закон Вигнера-Вея (Wigner-Way formula). Однако современные коды (например, RELAP5, TRACE, MELCOR) применяют более точные табличные данные из библиотек ORIGEN. Расчетчик обязан задать корректную историю выгорания топлива, чтобы получить изотопный состав. Именно изотопный состав определяет:

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

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

Детальное фото: Как рассчитывают последствия расплавления активной зоны

Моделирование кинетики окисления и плавления

Когда активная зона оголяется, начинается газовыделение и окисление. Ключевым параметром является температура оболочек твэлов. Если в реальном реакторе нет прямой оптической диагностики, то в расчете используются корреляции Ле Бланка (Le Blanc) или Бейкера-Джаста (Baker-Just).

Эти корреляции учитывают скорость окисления циркония в зависимости от температуры и давления пара. Расчет показывает:

Температура начала деформации оболочек: 700-800 °C.
Температура начала паро-циркониевой реакции: 1000-1100 °C.
Температура плавления циркония: 1850-1900 °C.
Температура плавления топливного сердечника (UO₂): около 2850 °C.

Важно понимать, что расплавление происходит не мгновенно. В моделях используется фазовый переход. Например, в коде MELCOR для топлива UO₂ задается доля расплавленной массы как функция энтальпии. Если энтальпия топлива превышает энтальпию плавления при данной температуре, считается, что материал перешел в жидкую фазу.

Результатом этого этапа является определение времени, через which активная зона полностью расплавится, и массы кориума, которая накапливается в нижней камере реактора.

Взаимодействие расплава с водой — расчет парового взрыва

Один из самых опасных сценариев — это попадание горячего расплава в остатки воды в днище корпуса. Это приводит к взрывному парообразованию (Fuel-Coolant Interaction, FCI).

Расчет FCI основан на моделировании фрагментации расплава. Когда струя кориума попадает в воду, она разбивается на мелкие частицы из-за гидродинамической неустойчивости. Чем мельче капли, тем больше площадь контакта и тем быстрее передается тепло. Процесс может быть мгновенным — в течение миллисекунд.

Ключевые параметры расчета:

Температура и вязкость расплава.
Глубина и температура водяного слоя.
Скорость падающей струи.
Состав расплава (доля оксидов против металлов).

Для количественной оценки взрыва используются коды TEXAS-V или MC3D. Они оценивают энерговыделение и давление ударной волны. В сценарии с большим количеством воды паровой взрыв может разрушить корпус реактора и систему гермооболочки. Однако при сухой нижней камере этот этап исключается из расчета.

Взаимодействие расплава с бетоном (MCCI)

Если расплав пробивает корпус реактора, он попадает на бетонное основание гермооболочки. Это взаимодействие описывается моделью Molten Core-Concrete Interaction (MCCI).

При взаимодействии с бетоном происходят два критических процесса:

Абляция бетона. Высокотемпературный расплав (2000-2500 °C) плавит бетон. Бетон разлагается с выделением CO₂, H₂O и SiO₂. Эти газы могут барботировать через расплав, увеличивая его теплопередачу и перемешивая.
Выделение неконденсируемых газов. CO₂ и водяной пар дополнительно поднимают давление в гермооболочке. Если система фильтрации или сброса давления не справляется, гермооболочка разрушается.

Расчет MCCI требует задания состава бетона (силикатный, известковый или серпентинитовый) и температуры кориума. Коды MEDICIS или CORQUENCH вычисляют скорость эрозии бетона (обычно сантиметры в час) и время до разрушения гермооболочки.

Пример из практики: при расчете запроектной аварии для реактора ВВЭР-1000 при сухой ловушке расплава (взаимодействие с бетоном) глубина проплавления может составить 0.5-1.5 метра в зависимости от времени подачи воды в защитную оболочку.

Транспорт продуктов деления и оценка радиационных последствий

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

Расчет базируется на группировке изотопов по летучести:

Инертные газы (1-я группа). Криптон (Kr-85), ксенон (Xe-133). Покидают топливо мгновенно при его повреждении. Не осаждаются, не задерживаются фильтрами. Полностью поступают в атмосферу при разгерметизации.
Летучие элементы (2-я группа). Йод (I-131, I-133), цезий (Cs-137, Cs-134). Могут осаждаться на поверхности, сорбироваться на аэрозолях. Их химическая форма (например, CsI или I₂) критически влияет на задержку в фильтрах.
Тугоплавкие элементы (3-я группа). Стронций (Sr-90), рутений (Ru-106), лантаноиды. Выделяются только при полном расплавлении топлива до 2800 °C. В основном остаются в виде аэрозолей.

Для расчета транспорта используются модели M2 (MELCOR) или ASTEC. В них учитываются процессы нуклеации, коагуляции, осаждения на стенках (brownian diffusion, gravitational settling), а также химические реакции с паром и бетоном. Итогом является источник радиоактивности — масса каждого изотопа, выходящая в окружающую среду в единицу времени.

Вероятностный анализ безопасности и неопределенности

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

Поэтому на финальном этапе используется вероятностный подход (Probabilistic Safety Analysis, PSA уровня 2 и 3). Расчётчик задает распределения для ключевых параметров:

Время течи.
Работоспособность аварийных насосов.
Эффективность спринклерной системы.
Скорость отбора проб.

Затем методом Монте-Карло (многовариантный расчет) получают ансамбль возможных исходов. Это позволяет оценить не только средние значения выброса, но и его 95-й процентиль (худший случай). Регуляторам (NDA, ASN, Rostechnadzor) именно такие результаты предоставляются для обоснования защитных мер, например, для расчета зоны планирования защитных мероприятий (PAZ и UPZ).

Примеры расчетных кодов и их роли

Для практического понимания стоит перечислить инструменты, используемые в этой области:

MELCOR — интегральный код для анализа тяжелых аварий, разработанный Sandia National Laboratories (США). Позволяет моделировать весь сценарий — от начала течи до выброса в атмосферу. MAAP — код, аналогичный MELCOR, разработанный EPRI. Используется для быстрых параметрических оценок. ASTEC — европейский код (IRSN, GRS), который особенно силен в детальном моделировании химии йода и рутения. RELAP/SCDAPSIM — используется для детального анализа теплогидравлики первого контура и механики разрушения твэлов, хотя в последних версиях также моделирует плавление.

Выбор кода часто определяется регуляторными требованиями страны. Например, во Франции обязателен ASTEC, в США — MELCOR, в России — СОКРАТ (разработка ИБРАЭ РАН).

Заключение: от теории к практике

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

Каждый этап расчета — от кинетики окисления циркония до абляции бетона — подчиняется законам физики, верифицирован на экспериментальных стендах (Phebus FP, LOFT, ITER) и требует от специалиста глубокого понимания не только математических моделей, но и реальной конструкции реактора.

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

Сводная таблица данных

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

Фаза аварии / Параметр	Температурный диапазон / Значение	Описание / Результат расчета
Остаточное тепловыделение (сразу после остановки)	6-7% от номинальной мощности	Энергия от распада короткоживущих продуктов деления.
Остаточное тепловыделение (через 1 час)	~1%	Снижение мощности остаточного энерговыделения.
Остаточное тепловыделение (через 1 сутки)	менее 0.5%	Дальнейший спад остаточного тепловыделения.
Начало деформации оболочек твэлов	700-800 °C	Температурный порог начала механических изменений циркониевых оболочек.
Начало паро-циркониевой реакции	1000-1100 °C	Экзотермическая реакция циркония с водяным паром; выделяется водород.
Температура плавления циркония	1850-1900 °C	Переход материала оболочек в жидкую фазу.
Температура плавления топливного сердечника (UO₂)	~2850 °C	Образование кориума (расплав оксидов урана, циркония и продуктов деления).
Температура кориума при взаимодействии с бетоном (MCCI)	2000-2500 °C	Диапазон температур расплава, вызывающий абляцию бетона.
Глубина проплавления бетона (пример для ВВЭР-1000)	0.5-1.5 метра	Зависит от времени подачи воды в защитную оболочку при сухой ловушке расплава.
Скорость эрозии бетона (MCCI)	сантиметры в час	Типичная расчетная скорость абляции бетонного основания.

Частые вопросы по теме (FAQ)

Какие температуры являются пороговыми для разных фаз расплавления активной зоны?

В расчетах фиксируются следующие критические температуры: начало деформации оболочек твэлов происходит при 700–800 °C, начало паро-циркониевой реакции — при 1000–1100 °C, плавление циркония — при 1850–1900 °C, а плавление топливного сердечника (UO₂) — при температуре около 2850 °C.

Как рассчитывается источник тепла для моделирования расплавления?

После остановки реактора остаточное тепловыделение от продуктов деления составляет 6–7% от номинальной мощности, через час — около 1%, а через сутки — менее 0,5%. Расчет ведется по закону Вигнера-Вея или с использованием точных таблиц из библиотек ORIGEN, где учитывается история выгорания и изотопный состав топлива.

Что такое взаимодействие расплава с бетоном (MCCI) и как его рассчитывают?

Если расплав пробивает корпус реактора и попадает на бетонное основание гермооболочки, высокотемпературный кориум (2000–2500 °C) начинает плавить бетон с выделением CO₂ и водяного пара. Это называется Molten Core-Concrete Interaction. Расчет ведется в кодах MEDICIS или CORQUENCH, которые вычисляют скорость эрозии бетона (обычно сантиметры в час) и время до разрушения гермооболочки. Для примера, при аварии на ВВЭР-1000 глубина проплавления может составить 0.5–1.5 метра.

Как оценивается взрывное парообразование при попадании расплава в воду?

Этот сценарий смоделирован как Fuel-Coolant Interaction (FCI). Ключевые параметры: температура и вязкость расплава, глубина и температура воды, скорость струи и состав расплава. При попадании в воду струя кориума фрагментируется на мелкие частицы, что создает огромную площадь теплообмена. Процесс длится миллисекунды. Для количественной оценки используются коды TEXAS-V или MC3D, которые вычисляют энерговыделение и давление ударной волны.

Какие группы радиоактивных продуктов учитываются при расчете радиационных последствий выброса?

Расчет базируется на группировке по летучести: 1-я группа — инертные газы (Kr-85, Xe-133), покидают топливо мгновенно и полностью поступают в атмосферу; 2-я группа — летучие элементы (I-131, Cs-137), могут осаждаться и сорбироваться на аэрозолях; 3-я группа — тугоплавкие (Sr-90, Ru-106), выделяются только при полном расплавлении топлива и в основном остаются в виде аэрозолей.