6-я МНТК «Обеспечение безопасности АЭС с ВВЭР»
ОКБ «ГИДРОПРЕСС», Подольск, Россия
26-29 мая 2009 г.
анализ АВАРИЙ с плавлением активной зоны РУ ВВЭР
по кодам СОКРАТ/В1 и RELAP/SCDAPSIM
В.В. Щеколдин, С. И. Пантюшин, Ю.С. Сорокин
ОКБ «Гидропресс», г. Подольск
Введение
В докладе представлен анализ внутрикорпускной стадии тяжелой аварии, включая исследование поведения кориума в напорной камере реактора и оценку параметров кориума, выходящего из корпуса реактора после его разрушения с использованием зарубежного кода RELAP/SСDAPSIM и российского тяжелоаварийного кода СОКРАТ/В1.
В качестве характерных для РУ ВВЭР-1000 рассмотрены исходные события:
Большая течь теплоносителя с отказом активной части САОЗ;
Полное обесточивание станции.
При выполнении анализа предполагается, что начальное состояние соответствует работе РУ на номинальной мощности. Рассматриваются процессы, начиная с исходного события, и далее, включая потерю теплоносителя, разогрев и деградацию активной зоны, разрушение внутрикорпусной шахты и корпуса реактора, выход кориума из корпуса реактора.
Для ОКБ «ГИДРОПРЕСС» наибольший интерес представляет анализ внутрикорпусной стадии тяжелой аварии с оценкой времени наступления характерных событий, поведения основных параметров РУ. Несмотря на постоянное совершенствование тяжелоаварийных кодов (включая рассматриваемые в данном докладе), до сих пор существуют трудности с обоснованием достоверности результатов расчетов. Все еще существуют значительные неопределенности в понимании и описании явлений, имеющих место в тяжелых авариях. В связи с этим необходимо настроить параметры кодов и согласовать поведение параметров РУ до перехода в тяжелую стадию аварии по аттестованным теплогидравлическим кодам.
Коды RELAP/SCDAPSIM/MOD3.2 и СОКРАТ/В1 включены в состав библиотеки программ предприятия ОКБ «ГИДРОПРЕСС», наряду с другими тяжелоаварийными кодами (например MELSIM) введены в опытную эксплуатацию и использовались для выполнения проектных расчетов.
Системный код СОКРАТ/B1 является современным отечественным средством для рассмотрения тяжелой стадии ЗПА, включая процессы разрушения и плавления топлива твэлов активной зоны, образование и перенос неконденсирующихся газов, разрушения шахты и корпуса реактора. В настоящее время код СОКРАТ/B1 проходит процедуру аттестации в Ростехнадзоре.
Отличительной особенностью кода RELAP/SCDAPSIM/MOD3.2 является его идентичность в теплогидравлической части с известным кодом наилучшей оценки RELAP5/MOD3.2. Вместе с тем имеется возможность анализа тяжелых аварий для сценариев с плавлением и разрушением активной зоны. Таким образом, можно провести оценку влияния точности моделирования начальной стадии тяжелой аварии на явления, происходящие при разрушении активной зоны.
Проведенная работа позволила разработать файлы исходных данных для расчетных кодов RELAP/SCDAPSIM и СОКРАТ/В1 для выполнения анализов тяжелых аварий РУ ВВЭР-1000. Сравнительный расчет позволил провести проверку и отладку файлов исходных данных для обоих расчетных кодов, как в теплогидравлической, так и в тяжелоаварийной части.
При условии настройки начальной стадии аварии по кодам улучшенной оценки (такие как ТРАП-97 и КОРСАР/ГП), учитывая достигнутую степень адаптации к РУ ВВЭР, можно применять коды СОКРАТ/В1и RELAP/SCDAPSIM для выполнения анализов тяжелых аварий при обосновании безопасности АЭС с РУ ВВЭР.
Описание кодов и методика расчета
Реальные физические системы моделируются в коде RELAP/SCDAPSIM/MOD3.2 [1] при помощи гидродинамических и термодинамических моделей. Гибкая схема нодализации предоставляет широкие возможности для моделирования компонентов РУ. Для моделирования гидроемкостей используется специальный элемент, позволяющий учитывать процессы теплообмена между газом, стенками сосуда и водой при истечении борного раствора из гидроемкости.
Расчетная схема реактора представлена на рисунке 1. Напорная камера реактора моделируется опускным участком из четырех элементов, который соединяется с объемом, моделирующем нижнюю часть НКР.
В реакторе внутрикорпусные устройства моделируются как тепловые структуры и компоненты модуля SCDAP. Поведение кориума на дне корпуса реактора описывается с помощью модуля COUPLE. Активная зона представлена двумя каналами: первый – с твэлами средней и максимальной мощности, второй – тапловыделяющая сборка максимальной мощности. Обогреваемая часть активной зоны разделена на два параллельных канала, каждый из которых включает 10 объемов. С точки зрения энерговыделения в активной зоне выделено семь отдельных компонентов, описывающих поведение твэлов, стержней СУЗ, направляющих каналов. Моделирование энерговыделения в топливе в результате реакции деления основано на модели точечной кинетики. Нейтронно-физические характеристики активной зоны (коэффициенты реактивности по плотности теплоносителя, температуре топлива и др.) соответствуют значениям для номинальной мощности. Мощность остаточных энерговыделений соответствует ANS 79 1.
Первый контур представлен четырьмя петлями, трубчатка парогенераторов описывается несколькими теплогидравлическими каналами (рис. 2, рис. 3). По второму контуру моделируются паропроводы, питательные насосы с трактами подачи питательной воды, отсечные и сбросные устройства.
Элементы конструкций над активной зоной представлены тепловыми структурами, которые не могут разрушаться. Специальные компоненты модуля SCDAP для построения моделей этих конструкций (БЗТ и т.п.) не использовались, т.к. их применение приводило к аварийным остановам расчетов.
Рис. 1 Расчетная схема реактора
Рис. 2 Расчетная схема первой петли
Рис.3 Расчетная схема парогенератора по второму контуру
Описание кода СОКРАТ
Для анализа тяжелых аварий в данной работе применен программный комплекс СОКРАТ/B1 [2], в котором соединены преимущества системного подхода к анализу РУ и современные достижения в области моделирования отдельных элементов и процессов. СОКРАТ/В1позволяет моделировать физические процессы тяжелой аварии РУ ВВЭР на внутрикорпусной и внекорпусной фазах после выхода расплава активной зоны в ЗО. Выполненные верификационные исследования кода СОКРАТ/В1[3] на отечественных и зарубежных экспериментальных данных подтверждают способность кода адекватно описывать совокупность процессов и явлений, определяющих протекание тяжелой ЗПА в РУ ВВЭР.
Функциональное наполнение СОКРАТ/В1 имеет блочно-модульную структуру. Основными модулями являются РАТЕГ, СВЕЧА и ГЕФЕСТ. Системный двухжидкостный теплогидравлический модуль РАТЕГ позволяет создавать полные расчетные модели РУ, включая системы управления и регулирования, и моделировать их работу с учетом следующих процессов: течение двухфазного теплоносителя с примесью неконденсируемых газов, перенос тепла в элементах РУ, теплообмен теплоноситель–стенка, теплоперенос излучением и др. Модуль СВЕЧА предназначен для моделирования физико-химических процессов, существенных для описания явлений деградации активной зоны в ходе тяжелой аварии: окисление, плавление и перемещение материалов, генерация водорода, теплообмен излучением и др. Модуль ГЕФЕСТ служит для моделирования взаимодействия корпуса реактора с расплавом активной зоны. В нём моделируются теплопередача, перемещение и плавление материала, перемешивание и расслоение расплава, тепловая эрозия корпуса, пограничное взаимодействие горячего материала с теплоносителем, теплообмен излучением и другие явления и процессы.
Моделирование тяжелой аварии с осушением, разогревом и плавлением активной зоны является сложной задачей и требует учета широкого круга взаимосвязанных явлений. При подготовке исходных данных и проведении анализа с использованием кода СОКРАТ/В1 учтен опыт анализов тяжелых аварий для РУ ВВЭР-1000 с использованием различных тяжелоаварийных кодов, а также результаты работ [4] - [6].
На рисунках 4-8 представлены нодализационные схемы главных циркуляционных петель, реактора, парогенератора, системы компенсации давления, а также области разрыва ГЦТ.
Рис. 4 Нодализационная схема расчетной петли с компенсатором давления
|
|
Рис. 5 Нодализационная схема реактора
|
Рис. 6 Нодализационная схема парогенератора
|
|
|
Рис. 7 Нодализационная схема системы компенсации давления первого контура
|
Рис. 8 Моделирование разрыва ГЦТ
|
Стационарное состояние
Параметры РУ в стационарном состоянии при работе на номинальной мощности определяются методом установления, при этом имитируется работа систем нормальной эксплуатации. При установлении давления в первом контуре в расчетах по коду RELAP/SCDAPSIM используется граничное условие, состоящее из искусственно введенного в нодализационную схему парового объема с постоянными параметрами. Этот объем исключается из нодализационной схемы первого контура при переходе к расчету аварии.
В таблице 1 приведены параметры РУ в стационарном состоянии, определенные по кодам СОКРАТ/В1и RELAP/SCDAPSIM.
Таблица 1
Характеристики стационарного режима
Параметр
|
СОКРАТ/В1
|
RELAP/SCDAPSIM
|
ВВЭР-1000*
|
Мощность реактора, МВт
|
3000
|
3000
|
3000
|
Давление на выходе из активной зоны, МПа
|
15,69
|
15,78
|
15,7
|
Расход теплоносителя через активную зону, кг/c
|
17575
|
17540
|
17678
|
Максимальная температура в среднем канале топлива, К
|
1419
|
1362
|
-
|
Максимальная температура в среднем канале оболочек твэл, К
|
609
|
606
|
-
|
Давление в ГПК
|
6,11
|
6,08
|
6,08
|
Масса теплоносителя в реакторе, т
|
79,2
|
75,5
|
74
|
Масса теплоносителя в первом контуре, т
|
259,9
|
249,2
|
245
|
Масса теплоносителя в ПГ, т
|
45,9
|
47,1
|
38-47
|
Объем ПГ, м3
|
127,5
|
-
|
127,5
|
Уровень в КД, м
|
7,65
|
8,1
|
8,17
|
Расход пара на турбину, кг/c
|
-
|
1630
|
1633
|
Температура на входе в активную зону, К
|
565
|
563,7
|
564,2
|
Температура на выходе из активной зоны, К
|
595
|
593,0
|
594,2
|
* Примечание. Указаны номинальные значения
Можно отметить удовлетворительное совпадение расчетных значений с номинальными значениями параметров РУ ВВЭР-1000. Таким образом, сравнение показало, что созданные пакеты (файлы) исходных данных позволяют корректно определить теплогидравлические параметры как отдельного оборудования, так и РУ в целом. Различия между сраниваемыми значениями, приведенными в таблице 1, могут быть признаны допустимыми в пределах погрешности их определения и не оказывают заметного влияния на результаты анализа аварии.
Анализ аварии «Большая течь с отказом активной части САОЗ»
В качестве исходного события рассмотрен мгновенный гильотинный разрыв ГЦТ при наложении полного обесточивания станции. Рассматриваемая авария характеризуется быстрым осушением реактора и активной зоны, что приводит к ее разогреву и плавлению. Хронология событий при аварии приведена в таблице 2.
Таблица 2
Хронология событий при аварии разрыва ГЦТ
с отказом активной части САОЗ
Событие
|
Время, с
|
СОКРАТ/В1
|
RELAP/SCDAPSIM
|
Исходное событие
|
0,0
|
0,0
|
Срабатывание АЗ
|
1,54
|
2,0
|
Срабатывание емкостей САОЗ
|
7,7
|
7,3
|
Достижение первого максимума температуры оболочек твэл в среднем канале
|
6,7 (1100 К)
|
3,0 (838К)
|
Начало разогрева среднего канала активной зоны (>880К)
|
558
|
312
|
Температура теплоносителя на выходе из активной зоны превышает 880К
|
1370
|
730
|
Температура максимального теплонапряженного участка твэлов в среднем канале (>1473К)
|
765
|
688
|
Начало плавления стержней СУЗ
|
884
|
1071
|
Начало пароциркониевой реакции
|
728
|
660
|
Начало плавления оболочек твэл
|
2284
|
1094(1941)
|
Формирование бассейна расплава в аз
|
-
|
1506
|
Формирование бассейна расплава в НКР
|
2846
|
-
|
Капли расплава попадают в НКР
|
2043
|
1507
|
Массовая деградация активной зоны и выход расплава в НКР
|
2657
|
2020
|
Разрушение шахты реактора
|
2912
|
-
|
Разрушение корпуса реактора
|
4472
|
6100
|
При выполнении анализа использован реалистический подход:
- начальные условия АЭС соответствуют режиму нормальной эксплуатации без учета возможных отклонений и неопределенностей в параметрах, величинах уставок и т.д.;
- характеристики активной зоны (коэффициенты реактивности, коэффициенты неравномерности и т.д.) принимаются в соответствии с нейтронно-физическими расчетами без учета неопределенностей и погрешности в расчетах;
- не учитываются отказы оборудования и ошибки персонала сверх установленных по сценарию;
- анализ выполняется с использованием современных компьютерных кодов «наилучшей оценки» («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 с уже не может остановить плавление.
Суммарный расход от САОЗ
Рис. 9 Большая течь с отказом САОЗ
|
Давление на выходе из активной зоны
Рис. 10 Большая течь с отказом САОЗ
|
Уровень в активной зоне
Рис. 11 Большая течь с отказом САОЗ
|
Суммарный расход теплоносителя в течь
Рис. 12 Большая течь с отказом САОЗ
|
Максимальная температура в среднем канале активной зоны
Рис. 13 Большая течь с отказом САОЗ
|
Мощность энерговыделений в активной зоне
Рис. 14 Большая течь с отказом САОЗ
|
Масса теплоносителя выброшенного в течь
Рис. 15 Большая течь с отказом САОЗ
|
Масса водорода
Рис. 16 Большая течь с отказом САОЗ
|
Масса расплавленного топлива в НКР
Рис. 17 Большая течь с отказом САОЗ
|
Максимальная температура оболочек твэл.
1- с отказом активной части САОЗ
2- с подключением канала САОЗ
Рис. 18 Большая течь
|
Анализ аварии «Полное обесточивание станции»
Тяжелая авария в результате полного обесточивания станции характеризуется протеканием процессов осушения, разогрева и разрушения активной зоны при высоком давлении. Хронология событий аварии приведена в таблице 3. Длительность различных стадий аварии представлена в таблице 4.
После обесточивания станции, срабатывает аварийная защита, мощность реактора снижается до уровня остаточных тепловыделений и остаётся на уровне, определяемом остаточными тепловыделениям продуктов деления топлива до конца рассмотренного периода аварии.
Начинается выбег ГЦН, отвод остаточных тепловыделений осуществляется за счет работы БРУ-А. После снижения уровня воды в ПГ, уменьшается эффективность отвода тепла, давление в первом контуре начинает расти до открытия ИПУ КД. Потеря теплоносителя первого контура приводит к ухудшению охлаждения, нагреву и разрушению активной зоны.
Рассмотрены два сценария данной аварии: с учетом действий оператора и без учета действий оператора. Предполагается, что оператор открывает ИПУ КД на 5000 с с целью снижения давления в первом контуре, чтобы обеспечить срабатывание гидроемкостей САОЗ, и подачу воды в напорную и сборную камеры реактора. Работа гидроемкостей САОЗ ограничивает величину опорожнения реактора и к моменту отключения гидроемкостей обеспечивает существенное заполнение активной зоны и реактора охлаждающей водой.
Суммарный расход от САОЗ, приведенный на рисунке 30, отличается, что влияет на поведение параметров РУ, т.к. процесс развития аварии во многом определяется количеством теплоносителя в первом конуре. В расчетах по коду СОКРАТ/В1 гидроемкости опорожняются довольно быстро, тогда как по коду RELAP/SCDAPSIM подача борного раствора из гидроемкостей происходит периодически, что приводит к увеличению длительности стадии разогрева активной зоны и сдвигает время начала тяжелого повреждения активной зоны примерно на 6000 с.
После прекращения подачи воды из гидроемкостей происходит снижение объема и массы воды в первом контуре и осушение активной зоны. Осушение активной зоны приводит к её разогреву, плавлению и разрушению, окислению оболочек твэлов, выходу радиоактивных продуктов деления из-под оболочки.
В отличие от предыдущей аварии, в данном режиме процесс разрушения активной зоны протекает медленней, обрушение топлива происходит периодически, порциями (рис. 19).
Таблица 3
Хронология аварии «Полное обесточивание станции»
Событие
|
момент времени, с
|
СОКРАТ/В1
|
RELAP/SCDAPSIM
|
вар1
|
вар2*
|
вар1
|
вар2*
|
Исходное событие
|
0,0
|
0,0
|
0,0
|
0,0
|
Срабатывание АЗ
|
1,54
|
1,54
|
2,1
|
2,1
|
Опустошение ПГ (HПГ<0,2м)
|
4120
|
4120
|
6700
|
6700
|
Открытие ИПУ КД
|
4560
|
4560
|
3387
|
3387
|
Начало разогрева активной зоны (>600°C) (средний канал)
|
7705
|
10842
|
10000
|
12000
|
Температура теплоносителя на выходе из активной зоны превышает 600°C
|
8767
|
12562
|
-
|
11500
|
Начало плавления стержней СУЗ
|
9805
|
11545
|
11821
|
16300
|
Начало пароциркониевой реакции
|
8130
|
11355
|
11200
|
11600
|
Начало плавления оболочек твэл
|
9825
|
17580
|
11970
|
17600
|
Формирование бассейна расплава
|
-
|
-
|
12130
|
18400
|
Капли расплава попадают в НКР
|
9225
|
15550
|
12405
|
18600
|
Выход расплава за пределы активной зоны (массовое обрушение)
|
9835
|
18340
|
12600
|
25200
|
Формирование бассейна расплава в НКР
|
10542
|
18685
|
-
|
-
|
Разрушение шахты реактора
|
12436
|
19417
|
-
|
-
|
Разрушение корпуса реактора
|
15577
|
23956
|
18500
|
-
|
* - вар2 - управление аварией - открытие оператором контрольного ИПУ КД (RELAP/SCDAPSIM) и контрольного и рабочего ИПУ КД (СОКРАТ/В1) на 5000 с.
Таблица 4
Длительность стадий разрушения
Событие
|
Время, с
|
Полное обесточивание станции
|
Большая течь
|
СОКРАТ/В1
|
RELAP/ SCDAPSIM
|
СОКРАТ/В1
|
RELAP/ SCDAPSIM
|
вар1
|
вар2
|
вар1
|
вар2
|
Осушение ПГ, с
|
4100
|
4100
|
6700
|
6700
|
-
|
-
|
Осушение реактора, с
|
8400
|
11000
|
11000
|
22000
|
600
|
600
|
Плавление и стекание активной зоны, с
|
400
|
1000
|
600
|
7500
|
500
|
1000
|
Появление дебриса в НКР- разрушение корпуса реактора, с
|
6000
|
8500
|
6000
|
-
|
2500
|
4500
|
Таблица 4 показывает, что длительность стадии осушения, оцененная по кодам СОКРАТ/В1и RELAP/SCDAPSIM, в аварии большая течь совпадают, тогда как длительность фаз разрушения активной зоны и формирования расплава существенно различаются. В коде RELAP/SCDAPSIM расплав формируется в горячей части активной зоны (начиная со стержней СУЗ) и постепенно распространяется, как на участки с более низким энерговыделением в этом же горячем канале, так и на соседние каналы с твэлами, имеющими меньшее энерговыделение. Соответственно, и разрушение активной зоны происходит постепенно (рис. 14, 17). В коде СОКРАТ/В1не использовалось разбиение зоны на твэлы разной мощности и процесс обрушения активной зоны протекает существенно быстрее, что существенно снижает время расчета.
Более детально проработанная модель ПГ в коде RELAP/SCDAPSIM приводит к существенно более позднему осушению ПГ (примерно на 2500 с, табл. 3). Примерно эта разница во времени между моментами наступления различных процессов и сохраняется на протяжении всей аварии обесточивания без действий оператора (вар1).
В коде СОКРАТ/В1идет более быстрая потеря воды второго контура, что связано с моделированием парогенератора. Расход естественной циркуляции по петле показывает хорошее согласование по перепадам давлений на различных участках первого контура по двум кодам.
Поведение параметров РУ в данной аварии приведено на рисунках 19-30. Как видно из рисунка 30 имеются отличия в поведении расхода от гидроемкостей САОЗ, рассчитанном по кодам СОКРАТ/В1и RELAP/SCDAPSIM. Рисунок 20 позволяет заключить, что полная потеря теплоносителя описывается по этим кодам качественно согласованно. Мощность остаточного тепловыделения по коду RELAP/SCDAPSIM рассчитывается, а по коду СОКРАТ/В1задается табличной зависимостью от времени, поэтому несколько выше, чем расчетные данные. Модуль нейтронной кинетики будет включен в следующую версию кода СОКРАТ/В2.
Расход из гидроемкостей САОЗ влияет на залив активной зоны, температурный режим твэлов и, следовательно, на максимальную температуру оболочек твэл. По данным анализа с использованием RELAP/SCDAPSIM к 30000 с корпус реактора не разрушился, тогда как по оценке с применением кода СОКРАТ/В1 разрушение корпуса наступило на 24000 с. Анализ сценария данной аварии без учета управления показал, что характерные времена протекания процессов, оцененные по вышеназванным кодам схожи, за исключением скорости осушения парогенераторов. Анализ сценария с учетом управления аварией выявил существенные различия в оценках времен наступления характерных событий в ходе аварии. С учетом опыта анализов тяжелых аварий с использованием этих кодов, можно предположить, что оценка времени разрушения корпуса реактора по коду СОКРАТ/В1 более адекватна.
Из анализа аварии по коду RELAP/SCDAPSIM следует, что для продления стадии разогрева активной зоны достаточно открыть лишь контрольное ИПУ КД, тогда как анализ по коду СОКРАТ/В1показал, что контрольного ИПУ недостаточно, необходимо открывать также рабочие ИПУ. Эти различия указывают на важность детального моделирования первого и второго контура, а также системы КД, для получения реалистической оценки длительности стадий тяжелой аварии.
Вместе с тем, хотя методика анализа поздних стадий тяжелых аварий по коду RELAP/SCDAPSIM требует доработки, разработанные пакеты исходных данных для кодов СОКРАТ/В1и RELAP/SCDAPSIM будут использованы в дальнейшем для создания на их основе пакетов адаптированных для выполнения анализов тяжелых аварий в обоснование безопасности проектов современных РУ ВВЭР.
Относительная мощность энерговыделений
Рис. 19 Полное обесточивание станции
|
Давление на выходе из активной зоны
Рис. 20 Полное обесточивание станции
|
Давление в ПГ2
Рис. 21 Полное обесточивание станции
|
Весовой уровень в ПГ2
Рис. 22 Полное обесточивание станции
|
Масса теплоносителя в первом контуре
Рис. 23 Полное обесточивание станции
|
Масса воды в ПГ2
Рис. 24 Полное обесточивание станции
|
Максимальная температура оболочек твэл
Рис. 25 Полное обесточивание станции
|
Расход через петлю 1
Рис. 26 Полное обесточивание станции
|
Расход через ИПУ КД
Рис. 27 Полное обесточивание станции
|
Уровень теплоносителя в активной зоне
Рис. 28 Полное обесточивание станции
|
Уровень теплоносителя в КД
Рис. 29 Полное обесточивание станции
|
Расход борного раствора из САОЗ
Рис. 30 Полное обесточивание станции
|
|
|
|