Материал: Белозеров В.И., Жук М.М., Кузина Ю.А., Терновых М.Ю. Физика и эксплуатационные режимы реактора ВВЭР-1000

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

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

Программа ПРОРОК-А предназначена для решения в интерактивном и автоматическом режимах задач формирования топливных загрузок реактора ВВЭР. Программа ПРОРОК-А является составной частью комплекса программ КАСКАД, опирается на его базу данных и для расчетов нейтронно-физических характеристик в процессе оптимизации использует модули интегрированных в комплекс программ БИПР-7А и ПЕРМАК-А.

На основании накопленного опыта проектирования процесс решения задачи оптимизации топливного цикла представляется в виде следующей последовательности шагов:

-выбор стационарной загрузки и ее оптимизация;

-выбор и оптимизация переходных топливных загрузок. Очевидно, что процесс решения задачи является в общем случае

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

276

предыдущих шагах. В частном случае оптимизации одной или двух топливных загрузок шаг 1 может отсутствовать. Однако, всегда существует режим стационарных перегрузок (возможно, рассчитанный заранее и называемый проектным), на который пользователь ориентируется при решении задачи на шаге 2 с точки зрения номенклатуры загружаемых ТВС, их количества и принципов расстановки ТВС в активной зоне.

Выбор стационарной загрузки и ее оптимизация. Выбор топ-

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

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

В программе ПРОРОК-А предусмотрена возможность интерактивного выбора режима стационарных перегрузок с периодом повторения топливных загрузок равным единице.

Выбор и оптимизация переходных топливных загрузок. По-

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

277

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

Задача оптимизации очередной топливной загрузки реактора представляется в виде последовательности следующих подзадач:

-перегрузка топлива, заключающаяся в выборе количества и типов загружаемых свежих ТВС, определении совокупности выгружаемых из активной зоны ТВС и замене их свежими ТВС;

-оптимизация расстановки кассет c точки зрения минимума некоторой целевой функции при выполнении ряда ограничений;

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

ианализ полученных результатов.

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

При оптимизации топливной загрузки в программе ПРОРОК-А используются интерактивные методы решения в сочетании со специально разработанными оптимизационными алгоритмами.

Программа ПИР-A не входит в комплекс программ «Каскад». Она предназначена для автоматизации сопоставления результатов эксплуатационных измерений нейтронно-физичес-ких характеристик энергоблоков ВВЭР-1000 и решает следующие задачи:

-имитационный расчет выгорания топливной загрузки по реальному графику нагрузки;

-восстановление поля энерговыделения в реакторе по экспериментальным замерам системы внутриреакторного контроля (по показаниям термопар и ДПЗ);

-сопоставление измеренного поля энерговыделения с расчетным.

Врезультате расчета по имитационному графику нагрузки определяются:

-энерговыработка ТВС на конец кампании; выгорание ДПЗ (протекший заряд);

-расчетное распределение энерговыделения в активной зоне для сопоставления с экспериментально измеренными величинами и для восстановления поля энерговыделения;

278

-критическая для состояний на мощности концентрация борной кислоты.

Расчеты энерговыделения для сопоставления и восстановления поля энерговыделения производятся со значениями интегральных характеристик состояния реактора на момент измерения.

Технология проведения расчетов следующая:

-расчет с восстановлением поля энерговыделения:

-расчет выгорания в прогнозном режиме;

-расчет выгорания в имитационном режиме;

-расчет по ПИР с восстановлением поля энерговыделения;

-расчет с простым сопоставлением расчета с измерениями:

-расчет выгорания в имитационном режиме;

-расчет по ПИР с сопоставлением поля энерговыделения.

8.18. Пэльный эффект

УТВС имеют направляющие каналы из циркониевого сплава марки Э110 наружным и внутренним диаметрами 13,0±0,05 и 11,0+01 мм соответственно, а также имеют наконечник, снабжённый четырьмя проливными отверстиями диаметром 2,0+04 мм, входя-

щими в ступенчатое центральное отверстие, одна ступень которого имеет диаметр 2,5+04 мм, а другая – диаметр 4,5+01 мм.

Центральный канал – перфорированный. Чувствительный элемент ТП размещается в полости цилиндрической части головки кассеты (входит вниз от края головки на 25 мм).

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

Потоки теплоносителя, выходящие из трубок ПЭЛ и из отверстий Ø10 мм, имеют характер отдельных струй и до места установки чувствительного элемента ТП (на длине примерно 100 мм) практически не перемешиваются.

Вплане относительно трубок ПЭЛ и отверстий Ø10 мм в различных кассетах термопара может занимать одно из шести положений (на окружности 0,87 мм через 60°), то есть, может находить-

279

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

На рис. 8.20 приведено разбиение активной зоны реактора ВВЭР-1000 по орбитам симметрии в зависимости от обогащения топлива и удалённости ТВС от центра активной зоны.

Рис. 8.20. Картограмма разбиения активной зоны ВВЭР-1000 по орбитам симметрии

280