Статья: Анализ аварий с плавлением активной зоны РУ ВВЭР по кодам СОКРАТ/В1 И RELAP/SCDAPSIM

При выполнении анализа использован реалистический подход:

- начальные условия АЭС соответствуют режиму нормальной эксплуатации без учета возможных отклонений и неопределенностей в параметрах, величинах уставок и т.д.;

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

- не учитываются отказы оборудования и ошибки персонала сверх установленных по сценарию;

- анализ выполняется с использованием современных компьютерных кодов “наилучшей оценки” (“best estimate”).

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

Вследствие интенсивного истечения теплоносителя давление внутри корпуса реактора резко снижается до давления насыщения при температуре теплоносителя в СКР. При достижении давления в реакторе до 5,9 МПа, происходит срабатыванию гидроемкостей САОЗ, и начинается подача воды в напорную и сборную камеры реактора. Работа ёмкостей ограничивает величину опорожнения реактора и к моменту их отключения обеспечивает существенное заполнение активной зоны и реактора охлаждающей водой. Активная зона охлаждается водой, поступившей из гидроемкостей САОЗ. По мере поступления охлаждающей воды САОЗ температура материалов активной зоны, в том числе топлива, снижается.

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

При снижении уровня теплоносителя в активной зоне ниже ее верха, за счет ухудшения теплоотвода от оболочек твэлов к парогазовому потоку, начинается постепенный рост их температуры. По мере роста температуры твэлов начинает расти давление газа под оболочкой, что приводит к их распуханию и разрыву, и авария переходит в стадию разрушения активной зоны. Разрыв оболочек обеспечивает доступ пара внутрь твэлов. Как следствие, происходит двустороннее окисление оболочек вследствие пароциркониевой реакции. В результате интенсивного окисления материалов активной зоны происходит генерация и поступление в объём ЗО водорода. Вследствие экзотермичности пароциркониевой реакции окисление оболочек сопровождается быстрым ростом температуры оболочек.

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

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

Изменение параметров РУ в данной аварии, определенное по кодам СОКРАТ/В1 и RELAP/SCDAPSIM показано на рисунках 9 - 18. В целом, можно отметить удовлетворительное качественное и, по отдельным параметрам, количественное совпадение. Важными параметрами являются скорость генерации и суммарный выход водорода в ЗО на внутрикорпусной стадии аварии, поскольку они определяют требования к системам водородной безопасности. Отличия, полученные для суммарного количества водорода, обусловлены учетом генерации водорода в коде СОКРАТ/В1 на всех стадиях тяжелой аварии.

Рисунки 9 и 10 позволяют заключить, что, несмотря на различные модели истечения (рис. 12, 15) полная потеря теплоносителя описывается по этим кодам взаимосолгасованно. Параметры истечения соответствуют выполненным ранее оценкам [4].

Уровень в активной зоне рис.11, по результатам, полученным по различным кодам, качественно согласуется, что отражается на поведении температуры оболочек твэл (рис. 13). Более раннее обрушение (рис. 14) элементов активной зоны по коду RELAP/SCDAPSIM вызвано присутствием “горячих” твэлов в расчетной схеме активной зоны, тогда как в коде СОКРАТ/В1 присутствуют только твэлы со средним энерговыделением по активной зоне.

Отличия в данных, характеризующих повреждения активной зоны, более значительные (рис. 16, 17), т.к. модели окисления, обрушения элементов активной зоны существенно различаются. Наличие бассейна расплава, сформировавшегося в зоне с более высоким энерговыделением, приводит к более раннему разрушению активной зоны. Отказ от использования SCDAP моделей ВКУ (БЗТ, шахта), также внесло свой вклад в данное различие, т.к. не учитывается окисление стали, нагрев излучением и плавление элементов ВКУ в использованной модели кода RELAP/SCDAPSIM.

Для оценки возможности залива активной зоны после достижения высокой температуры теплоносителя на выходе из активной зоны, а также для исследования эффективности залива был рассмотрен сценарий (по коду RELAP/SCDAPSIM) с подачей от одного канала активной части САОЗ через 700 с от начала аварии. Экспериментальные исследования [7,8] показали, что даже при превышении максимальной температуры оболочек твэл 1700К можно охладить сборку, используя расход, эквивалентный расходу от одного канала САОЗ.

Максимальная температура оболочек твэл для сценариев с учетом и без учета подключения одного канала активной части САОЗ приведена на рис. 18. Видно, что подача борного раствора резко ускоряет пароциркониевую реакцию и при таких температурах подача воды ускоряет процесс обрушения, тем не менее дальнейшее разрушение активной зоны удается предотвратить.

Однако, времени для подключения САОЗ у оператора очень мало. Вариантные расчеты показали, что уже через 750 с после начала аварии, т.е. спустя всего 60 с после превышения максимального проектного предела, работы одного канала САОЗ недостаточно для предотвращения разрушения активной зоны. Согласно рисунку 18, включение одного канала САОЗ через 700 с после начала аварии позволяет предотвратить дальнейшее разрушение активной зоны, т.е. задержка в подаче борного раствора всего на 50 с уже не может остановить плавление.

Анализ аварии “Полное обесточивание станции”

Тяжелая авария в результате полного обесточивания станции характеризуется протеканием процессов осушения, разогрева и разрушения активной зоны при высоком давлении. Хронология событий аварии приведена в таблице 3. Длительность различных стадий аварии представлена в таблице 4.

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

Начинается выбег ГЦН, отвод остаточных тепловыделений осуществляется за счет работы БРУ-А. После снижения уровня воды в ПГ, уменьшается эффективность отвода тепла, давление в первом контуре начинает расти до открытия ИПУ КД. Потеря теплоносителя первого контура приводит к ухудшению охлаждения, нагреву и разрушению активной зоны.

Рассмотрены два сценария данной аварии: с учетом действий оператора и без учета действий оператора. Предполагается, что оператор открывает ИПУ КД на 5000 с с целью снижения давления в первом контуре, чтобы обеспечить срабатывание гидроемкостей САОЗ, и подачу воды в напорную и сборную камеры реактора. Работа гидроемкостей САОЗ ограничивает величину опорожнения реактора и к моменту отключения гидроемкостей обеспечивает существенное заполнение активной зоны и реактора охлаждающей водой.

Суммарный расход от САОЗ, приведенный на рисунке 30, отличается, что влияет на поведение параметров РУ, т.к. процесс развития аварии во многом определяется количеством теплоносителя в первом конуре. В расчетах по коду СОКРАТ/В1 гидроемкости опорожняются довольно быстро, тогда как по коду RELAP/SCDAPSIM подача борного раствора из гидроемкостей происходит периодически, что приводит к увеличению длительности стадии разогрева активной зоны и сдвигает время начала тяжелого повреждения активной зоны примерно на 6000 с.

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

В отличие от предыдущей аварии, в данном режиме процесс разрушения активной зоны протекает медленней, обрушение топлива происходит периодически, порциями (рис. 19).

Таблица 3

Хронология аварии “Полное обесточивание станции”

* - вар2 - управление аварией - открытие оператором контрольного ИПУ КД (RELAP/SCDAPSIM) и контрольного и рабочего ИПУ КД (СОКРАТ/В1) на 5000 с.

Таблица 4

Длительность стадий разрушения

Таблица 4 показывает, что длительность стадии осушения, оцененная по кодам СОКРАТ/В1и RELAP/SCDAPSIM, в аварии большая течь совпадают, тогда как длительность фаз разрушения активной зоны и формирования расплава существенно различаются. В коде RELAP/SCDAPSIM расплав формируется в горячей части активной зоны (начиная со стержней СУЗ) и постепенно распространяется, как на участки с более низким энерговыделением в этом же горячем канале, так и на соседние каналы с твэлами, имеющими меньшее энерговыделение. Соответственно, и разрушение активной зоны происходит постепенно (рис. 14, 17). В коде СОКРАТ/В1не использовалось разбиение зоны на твэлы разной мощности и процесс обрушения активной зоны протекает существенно быстрее, что существенно снижает время расчета.

Более детально проработанная модель ПГ в коде RELAP/SCDAPSIM приводит к существенно более позднему осушению ПГ (примерно на 2500 с, табл. 3). Примерно эта разница во времени между моментами наступления различных процессов и сохраняется на протяжении всей аварии обесточивания без действий оператора (вар1).

В коде СОКРАТ/В1идет более быстрая потеря воды второго контура, что связано с моделированием парогенератора. Расход естественной циркуляции по петле показывает хорошее согласование по перепадам давлений на различных участках первого контура по двум кодам.