Опыт верификации прикладного программного обеспечения системы внутриреакторного контроля атомной электростанции «Куданкулам» в части реализации защиты по запасу до кризиса теплообмена

Размещено на http://www.allbest.ru/

ОПЫТ ВЕРИФИКАЦИИ ПРИКЛАДНОГО ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ СИСТЕМЫ ВНУТРИРЕАКТОРНОГО КОНТРОЛЯ АЭС «КУДАНКУЛАМ» В ЧАСТИ РЕАЛИЗАЦИИ ЗАЩИТЫ ПО ЗАПАСУ ДО КРИЗИСА ТЕПЛООБМЕНА

Автор доклада:

С.А. Еремеев

Руководитель темы:

М.А. Подшибякин

Верификация прикладного программного обеспечения (ППО) программного технического комплекса защиты (ПТК-З) системы внутриреакторного контроля (СВРК) АЭС «Куданкулам» проводилась в соответствии с требованиями Правил ядерной безопасности реакторных установок атомных станций НП-082-07 (п. 2.4.15) и требованиями технического задания на систему контроля управления и диагностики (СКУД).

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

Описание стенда для проведения верификации

В состав стенда для проведения верификации ППО ПТК-З входили два ноутбука моделей Sony Vaio VGN-FS315MR и Samsung P28. Характеристики ноутбуков представлены в таблице 1.

Таблица 1 Технические характеристики ноутбуков

Элементы	Samsung P28	Sony Vaio VGN-FS315MR
Центральный процессор	Intel Pentium M 1.5 ГГц	Intel Pentium M 1.73 ГГц
Оперативное запоминающее устройство	512 Мб	512 Мб
Жесткий диск	40 Гб	80 Гб
Оптический накопитель	DVD/CD-RW	DVD±RW
Видеокарта	128 Мб NVIDIA GeForce 6400	64 Мб ATI Radeon 9700

На ноутбуке модели Samsung P28 установлено ППО ПТК-З, которое выполняет следующие функции в части реализации защиты по запасу до кризиса теплообмена:

- обеспечивается прием массива данных об измерениях;

- прием измеренных соседними каналами сигналов;

- первичную обработку результатов измерения входных сигналов (отбраковку и сглаживание);

- расчет массовых расходов теплоносителя в петлях первого контура;

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

- расчет минимального запаса до кризиса теплообмена и формирование признака сигнала защиты по данному параметру;

- формирование сигнала АЗ при нарушении локальными параметрами заданных пределов;

- прием корректирующей и передачу диагностической информации и измеренных данных в сервисный компьютер по сети Ethernet.

На ноутбуке модели Sony Vaio VGN-FS315MR установлено сервисное ПО, которое выполняет следующие функции:

- формирование массивов данных об измерениях;

- предоставление информации о состоянии аппаратуры;

- отображение в реальном времени на экране монитора текущего состояния всех диагностических блоков;

- архивация данных на диске.

Принципиальная схема стенда для проведения верификации ППО ПТК-З приведена на рисунке 1.



Рис. 1 Принципиальная схема стенда для проведения верификации прикладного программного обеспечения программного технического комплекса защиты

Программа эксплуатации предназначена для обеспечения выполнения функций станции контроля нижнего уровня по сбору, обработке и представлению эксплуатирующему персоналу информации о состоянии аппаратуры ПТК-З.

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

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

Data Creator в свою очередь использует данные, созданные программой LENCO, которая перерабатывает результаты анализов безопасности в виде файла «lenta». LENCO преобразовывает записанные в файл данные об изменении линейного энерговыделения в токовые сигналы детекторов прямого заряда (ДПЗ) и выделяет следующие параметры:

- давление на выходе из реактора (1 шт.);

- температура теплоносителя в холодной нитке 1 контура i-ой петли (4 шт.);

- перепад давления на главном циркуляционном насосном агрегате (ГЦНА) i-ой петли (4 шт.);

- состояние ГЦНА i-ой петли (4 шт.)

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

Формирование сигналов защиты в ППО ПТК-З

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

Формирование сигналов защиты по локальным параметрам осуществляется двумя комплектами ПТК-З, каждый из которых состоит из трех каналов. Сигналы защиты, сформированные в каналах одного комплекта (три сигнала), поступают в оборудование первого комплекта аппаратуры логической обработки сигналов системы управления и защиты (АЛОС I), где обрабатываются по логике «два из трех». Аналогично, сигналы защиты из каналов другого комплекта ПТК-З поступают во второй комплект аппаратуры логической обработки сигналов системы управления и защиты (АЛОС II). Сигналы защиты, сформированные по логике «два из трех», в дальнейшем по логике «один из двух» поступают на управление приводами органов регулирования системы управления и защиты (СУЗ). Структурная схема формирования сигнала защиты приведена на рисунке 2.

Рис. 2 Структурная схема формирования сигнала защиты

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

Каналы ПТК-З переходят в основной режим работы после включения электропитания и работают асинхронно независимо друг от друга. Прием сигналов от собственных датчиков, обмен информацией между каналами по локальной сети нижнего уровня (сигналами ДПЗ, заводимых в разные каналы), расчет локальных параметров и выдача в АЛОС сигналов защиты при достижении параметрами допустимых пределов, в диапазоне мощности от 30 до 110 % от номинальной, осуществляется циклически.

Каждый канал ПТК_З имеет дублированную структуру, в которой блоки управления (БУП) и блоки ввода-вывода (БВВ) объединены дублированной магистралью. Ввод/вывод функции защиты осуществляется программным путем (автоматически) на уровне мощности реактора 30 % от номинальной, рассчитанной в ПТК-З. Сигналы защиты, сформированные в каналах ПТК-З, поступают в комплекты оборудования АЛОС, где обрабатываются по логике «два из трех» и далее поступают в комплекс электрооборудования (КЭ) СУЗ, где обрабатываются по логике «один из двух».

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

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

В случае отказа питания стойки ПТК_З, отказа питания плюс 24 В или отсутствия блока БВЦ-256Р (блок вывода релейных сигналов) на выходе ПТК-З формируется сигнал защиты, поступающий в соответствующий канал АЛОС.

Описание алгоритма расчетов переменных состояния РУ в ПТК-З

Перед расчетом физических параметров для каждого измеренного сигнала, кроме токов ДПЗ, выполняются процедуры отбраковки и проверки скорости изменения измеренного сигнала. Для нормированных сигналов после проведения процедуры отбраковки и проверки скорости дополнительно выполняется процедура сглаживания. Для токов ДПЗ после проведения процедуры отбраковки выполняется расчет нормализованных токов ДПЗ и затем для нормализованных токов ДПЗ, рассчитанных по фильтру Калмана-Цимбалова, выполняется процедура проверки скорости.

Массивы нормализованных токов ДПЗ поступают на обработку в блоки управления своей и других стоек ПТК-З. В указанных блоках проводятся расчеты линейной мощности максимально нагруженных тепловыделяющих элементов (ТВЭЛ) в 163 ТВС в семи слоях по высоте и сравнение с уставками для обеспечения защиты по линейной мощности максимально нагруженных ТВЭЛ в семи слоях по высоте активной зоны.

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

,

где KKij(0) - отношение линейного энерговыделения максимально нагруженного ТВЭЛ к среднему энерговыделению ТВЭЛ на i участке j ТВС в момент калибровки;

KVij(0) - относительная линейная мощность i участка j ТВС в момент калибровки; dfKVim (0) - относительная линейная мощность одного из ближайших к j ТВС i участка ТВС с m СВРД в момент калибровки (в месте реального размещения i ДПЗ, с учетом смещения сборок СВРД по высоте от проектного положения);

лim(0) - чувствительность i ДПЗ в m СВРД, приведенная к единице длины ТВС в момент калибровки;

NTVEL - число ТВЭЛ в ТВС;

(t) - текущее значение нормализованного тока i ДПЗ в m СВРД;

i = 17;

m = 154;

j = 1163.

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

Расчет статистически усредненной величины линейного энерговыделения максимально нагруженных ТВЭЛ в ТВС по показаниям ближайших ДПЗ проводится по формуле:

,

где Aijm - линейное энерговыделение максимально нагруженного ТВЭЛ на (i, j) участке активной зоны, рассчитанная по показаниям (i, m) ДПЗ;

qjm - статистический вес показаний m СВРД в линейной мощности j ТВС;

зim - достоверность величины (0, либо 1) определяется достоверностью сигнала (i, m) ДПЗ;

i = 17;

m = 154;

j = 1163.

Входными данными для расчета являются:

- значения нормализованных токов ДПЗ с признаками достоверности;

- номера ближайших ТВС с СВРД;

- статистические веса линейных энерговыделений максимально нагруженных ТВЭЛ, определенных в ближайших ТВС с СВРД.

Выходными данными являются:

- линейные энерговыделения максимально нагруженных ТВЭЛ на i участке (слое) j ТВС (i=17, j=1163).

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

Расчет массовых расходов теплоносителя в петлях первого контура производится по следующим формулам:

- при условии включенного ГЦН:

,

где Gi - массовый расход теплоносителя в i петле, кт/ч

Ai, Bi, Ci - коэффициенты, взятые из паспортных характеристик ГЦНА i петли;

i - плотность теплоносителя холодной нитки i петли, кг/м3;

fi - частота питания ГЦНА в i петле, Гц;

fi0 - тарировочная частота питания ГЦНА, Гц;

Pi - перепад давления на ГЦНА i петли, Па;

g - ускорение свободного падения, м/с2.

При расчетах признак достоверности устанавливается, если достоверен перепад давления на ГЦН, а также достоверна плотность теплоносителя. Кроме того, если какая-либо частота питания ГЦН имеет признак недостоверности, то ее значение заменяется средней частотой питания, рассчитанной по оставшимся достоверным частотам питания ГЦН. Если же не достоверна ни одна частота питания ГЦН, то все они заменяются тарировочным значением частоты питания ГЦН;

- при отключении ГЦН расчет массового расхода теплоносителя в петле первого контура производится по следующей формуле:

,

где i коэффициент гидравлического сопротивления обратного тока ГЦН в i петле;

FГЦН - площадь проходного сечения трубопровода, м2;

- при переходном режиме работы петли после отключения ГЦН в течение времени Твыбег, определяемого «выбегом» ГЦН, массовый расход теплоносителя во всех петлях первого контура определяется значением массового расхода теплоносителя до отключения ГЦН. По окончании выбега ГЦН расчет ведется по формуле обратного тока (4).

Значение времени Твыбег хранится в файле и может неоперативно корректироваться.

Проектное значение Твыбег = 30 с.

Расчет суммарного значения массового расхода производится суммированием всех значений массовых расходов теплоносителя в петлях по следующей формуле:

,

где G1K - суммарный массовый расход теплоносителя в петлях первого контура, кт/ч;

Gk - массовый расход теплоносителя в k-ой петле.

Выполняется расчет средней массовой скорости воды в ячейках ТВС (w) по формуле:

,

где Kпр - коэффициент, учитывающий протечки теплоносителя;

FТВС - площадь проходного сечения ТВС за вычетом центрального канала и каналов для ПЭЛ;

Kн - коэффициент, учитывающий неравномерность массовой скорости по сечению ТВС. Коэффициент Kн хранится в файле и может неоперативно корректироваться. Проектное значение Kн = 0,98.

Для каждой ТВС выполняется определение температуры теплоносителя на ее входе по формуле:

,

где j - номер ТВС;

- температура теплоносителя в холодной нитке k-ой петли;

Fk - признак включения ГЦНА в k-ой петле;

Gk - расход теплоносителя в k-ой петле;

m - код, равный F18 + F24 + F32 + F4;

Pkjm - коэффициент влияния температуры k-ой петли на входную температуру j-ой ТВС при m-ом состоянии значений признаков Fk..

Для каждого выбранного участка (i,j) вычисляется энтальпия Iij по формуле:

,

где Kt - коэффициент запаса по подогреву;

= I(Pвх,) - энтальпии на входе ТВС;

Hаз - высота активной зоны;

Pвх - давление на входе в активную зону;

Fяч - площадь проходного сечения рассматриваемой ячейки.

Для каждого выбранного участка (i,j) вычисляется относительная энтальпия Xij по формуле:

,

где I'(Pвх) - энтальпия воды на линии насыщения;

r(Pвх) - удельная теплота парообразования.

Для каждого выбранного участка (i,j) вычисляется величина критического теплового потока по формуле:

,

где fF - форм-фактор, вычисляемый по формуле:

,

если fF > 1,0, то fF = 1,0,

если fF < 0,85, то fF = 0,85,

где Pкр = 22,115 МПа - критическое давление воды.

Определяются запасы до кризиса теплообмена (DNBij) по формуле:

,

где Kf- погрешность формулы для критического теплового потока;

Kq- коэффициент запаса по мощности;

= 3,1416;

d = 0,0091м - диаметр твэла.

Из всех полученных значений DNBij выбирается минимальное DNBmin и сравнивается с установленными пределами (ПЗ, АЗ). При достижении или уменьшении ниже какого-либо предела выдается соответствующий сигнал защиты в систему АЗ-ПЗ.

Результаты верификации

Верификация ППО ПТК-З проводилась путем сравнения результатов расчетов, полученных по кодам ППО ПТК-З, с результатами анализов безопасности реакторной установки (РУ) АЭС «Куданкулам». При этом одним из условий проведения верификационных расчетов было соответствие исходных данных для ППО ПТК-З исходным данным, принятым в анализах безопасности. В частности: равенство профиля энерговыделения по высоте активной зоны, равенство температуры теплоносителя в холодных нитках, равенство давлений на выходе из реактора, равенство напоров ГЦН. В качестве приемочного критерия было принято время достижения уставки аварийной защиты параметром «минимальный запас до кризиса теплообмена» (МЗКТ) в ППО ПТК-З, которое должно было быть не позднее момента времени достижения уставки этим параметром в расчетных анализах безопасности.

Верификация проводилась в два этапа. Первый этап верификации был проведен по принципу «черного ящика». В результате было получено, что для ряда проектных режимов приемочный критерий не был выполнен. К таким режимам относятся: «нарушение или отказ регулятора давления пара, приводящее к уменьшению расхода пара»; «потеря внешней электрической нагрузки»; «потеря вакуума в конденсаторе и другие случаи, приводящие к останову турбины»; «нарушения в системе продувки-подпитки, приводящие к увеличению запаса теплоносителя первого контура в результате впрыска воды с температурой 20 оС»; «потеря нормального расхода питательной воды (за исключением разрыва трубопроводов питательной воды)»; «непреднамеренное закрытие отсечного клапана на паропроводе с последующей непосадкой импульсного предохранительного устройства (ИПУ) аварийного ПГ, приводящей к аварийному снижению давления в системе паропроводов свежего пара»; «непреднамеренное закрытие отсечного клапана на паропроводе»; «разрыв трубопровода питательной воды»; «нарушение в системе питательной воды, приводящее к увеличению расхода питательной воды»; «отключение различного числа ГЦНА». Некоторые результаты приведены на рисунках 3-6.

Рис. 3 Значение минимального запаса до кризиса теплообмена



Рис. 4 Значение минимального запаса до кризиса теплообмена

Отсутствие выполнения приемочного критерия объяснялось следующими причинами:

- применением формулы ОКБ «ГИДРОПРЕСС» /1/ для расчета критического теплового потока вне области её применения;

- некорректным расчетом расхода теплоносителя в аварийной петле в нестационарных режимах с отключением различного числа ГЦНА (рисунки 7 и 8).

Рис. 7 Расходы на входе в реактор при отключении одного ГЦНА из трех работающих (ТРАП)



Рис. 8 Расходы на входе в реактор при отключении одного ГЦНА из трех работающих (ПТК-З)

В результате были разработаны следующие основные требования по корректировке ППО ПТК-З:

1) за пределами области применимости корреляции ОКБ «ГИДРОПРЕСС» для области недогретой воды и малых паросодержаний применить формулу Смолина /2/:

,

где i' - энтальпия воды на линии насыщения, Дж/кг;

i'' - энтальпия пара на линии насыщения, Дж/кг;

с' - плотность воды на линии насыщения, кг/м3;

с'' - плотность пара на линии насыщения, кг/м3;

у - коэффициент поверхностного натяжения воды при температуре на линии насыщения, Н/м;

м' - динамическая вязкость воды на линии насыщения, Па•с;

м" - динамическая вязкость пара на линии насыщения, Па•с;

g - ускорение свободного падения, м/с2;

x - относительная энтальпия;

сщ - массовая скорость воды в ячейках ТВС, кг/(с*м2).

При расчете свойств воды и пара на линии насыщения, используемых в методике Смолина, были применены рассчитанные методом наименьших квадратов следующие зависимости (полиномы четвертой степени):

Динамическая вязкость воды на линии насыщения (Па•с):

м' = 1,025•10-4 + 1,759•10-10•Р - 4,161•10-7•Р2 + 2,672•10-8•Р3-5,848•10-10•Р4,

где Р - давление, МПа.

Динамическая вязкость пара на линии насыщения (Па•с):

м" = 1,664•10-5 + 5,801•10-9•Р + 7,875•10-8•Р2 - 6,01•10-9•Р3 + 1,704•10-10•Р4.

Плотность воды на линии насыщения (кг/м3):

с' = 799,13857 + 2,81899•Р - 2,59605•Р2 + 0,15691•Р3 - 3,62191•10-3•Р4.

Плотность пара на линии насыщения (кг/м3):

с'' = 75,70334 - 19,69858•Р + 3,00225•Р2 - 0,15772•Р3 + 3,4228•10-3•Р4.

Коэффициент поверхностного натяжения воды при температуре на линии насыщения (Н/м):

у = 0,03958 - 4,37034•10-3•Р + 2,1691•10-4•Р2 - 6,6697•10-6 Р3 + 9,61386•10-8•Р4.

Энтальпия воды на линии насыщения (Дж/кг):

i'= 9,37167•105 + 4,44044•104•Р + 1345•Р2 - 148,29622•Р3 + 4,04533•Р4.

Энтальпия пара на линии насыщения (Дж/кг):

i'' = 2,66578•106 + 5,07655•104•Р - 7377,32•Р2 + 362,78584•Р3 - 7,7 0182•Р4.

Погрешность определения параметров не превышает 0,1 %;

2) учесть выбег ГЦНА при его отключении. Расход в петле с неработающим ГЦНА в процессе его выбега может быть определен из решения уравнения:

,

где Т - постоянная времени, с;

g - расход в петле с отключенным ГЦНА, отн. ед.;

gр - расход в петле с работающим ГЦНА (усредненное значение по петлям с работающими ГЦНА), отн. ед.;

Дh = hр - ho - разность перепадов давления на ГЦНА, отн. ед.;

hр - перепад давления на работающем ГЦН (усредненное значение по петлям с работающими ГЦНА, отнесенное к среднему по петлям перепаду давления на ГЦНА в номинальном режиме), отн. ед.;

ho - перепад давления на отключенном (неработающем) ГЦНА (отнесен к среднему по петлям перепаду давления на ГЦНА в номинальном режиме), отн. ед.;

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

бо - коэффициент, характеризующий гидравлическое сопротивление петли с отключенным ГЦНА, без учета потерь в реакторе, отн.ед.

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

Условием для расчета выбега является изменение сигнала о состоянии ГЦНА (включен /отключен).

Переход на расчет расхода по формуле для обратного тока должен происходить по окончании выбега ГЦНА. Для всех режимов с отключением ГЦНА рекомендуется принять время выбега равным 100 с.

Уравнение расчета выбега ГЦНА является уравнением Риккати и при постоянном в пределах шага счета значении правой части имеет аналитическое решение.

Для удовлетворительного совпадения результатов расчета выбега по данному уравнению с результатами, приведенными в анализах безопасности, достаточно иметь два набора коэффициентов: один для режимов с отключением одного из четырех ГЦНА и другой для режимов с отключением одного из трех, двух из четырех или трех из четырех ГЦНА.

На рисунках 9-11 приведены результаты расчета выбегов по предлагаемой методике и рассчитанной по коду «ТРАП» для различных вариантов отключения ГЦНА.



Рис. 9 Сравнение расходов на входе в реактор при отключении одного ГЦНА из четырех работающих

Рис. 10 Сравнение расходов на входе в реактор при отключении двух ГЦНА из четырех работающих



Рис. 11 Сравнение расходов на входе в реактор при отключении одного ГЦНА из трех работающих

Режимы с отключением одного из четырех ГЦНА

Т = 2,5с; бо = 0,25; бр = 0,25 - при g > 0.58;

Т = 4,5с; бо = 0,25; бр = 0,25 - при 0,58 ? g > -0.05;

Т = 8,0с; бо= 1,546; бр = 0,25 - при g ? -0.05.

Режимы с отключением одного ГЦНА из трех или двух ГЦНА из четырех работающих

Т = 3,0с; бо = 0,25; бр = 0,25 - при g > 0.58;

Т = 8,5с; бо = 0,25; бр = 0,25 - при 0,58 ? g > -0.05;

Т = 16,0с; бо= 3,633; бр = 0,25 - при g ? -0.05.

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

Значения коэффициентов, определяемые по результатам расчета стационарных и переходных режимов с отключением ГЦНА, должны быть уточнены при проведении испытаний в процессе ПНР с учетом реальных характеристик ГЦНА и первого контура.

Предлагаемая методика расчета выбега ГЦНА может быть использована и для непроектных выбегов, поскольку при его расчете используются реальные значения напоров ГЦНА.

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

Как следует из приведенных результатов, учет дополнительных факторов (выбег ГЦНА и методика Смолина) в алгоритмах расчета МЗКТ в ПТК-З, существенно повышает точность расчета без привлечения дополнительных входных сигналов.

После внесения соответствующих изменений в ППО ПТК-З был выполнен второй этап верификации. Повторно, аналогично первому этапу, были рассмотрены все проектные режимы, требующие срабатывания аварийной защиты по запасу до кризиса теплообмена. На рисунке 12, в качестве примера приведены результаты расчета для режима «потеря внешней электрической нагрузки», приведены результаты расчета МЗКТ по ППО ПТК-З до (верифицируемая версия) и после (с методикой Смолина) корректировки.

Рис. 12 Значение минимального запаса до кризиса теплообмена

Из проведенных по второму этапу верификации расчетов следует, что МЗКТ, рассчитываемый по ППО ПТК-З, всегда превышает аналогичный параметр, полученный в анализах безопасности. При этом отклонения текущих значений МЗКТ, рассчитываемых в ППО ПТК-З, от аналогичного параметра в анализах безопасности связаны с различиями в алгоритмах расчета, основные из которых следующие:

- различное число расчетных участков по высоте активной зоны;

- в ППО ПТК-З параметры и свойства воды (пара) на расчетных участках рассчитываются при давлении на выходе из реактора (в анализах безопасности используется локальное давление);

- в ППО ПТК-З массовая скорость в горячем канале принимается пропорциональной средней по активной зоне массовой скорости теплоносителя с коэффициентом пропорциональности 0,98;

- в ППО ПТК-З теплопередача от топлива к теплоносителю не рассчитывается, а принимается равной линейному энерговыделению.

Таким образом, по результатам двух этапов верификации ППО ПТК-З был сделан вывод о том, что приемочный критерий может быть выполнен для откорректированного ППО ПТК-З при условии внесения в ППО ПТК-З корректирующего множителя Kкор в формуле расчета запаса до кризиса теплообмена. При этом для выполнения приемочного критерия при сохранении проектного значения уставки аварийной защиты 1,35 значение корректирующего множителя составляет 0,849. Полученное значение может быть уточнено после разработки и использования нормативной методики верификации программных средств управляющих систем нормальной эксплуатации и управляющих систем безопасности.

Список литературы

1. Программа для ЭВМ. Расчет нестационарных режимов энергетических установок с ВВЭР, ДИНАМИКА-97, 8624607.00467, ОКБ «Гидропресс», 1998.

2. Смолин В.Н., Поляков В.К. Семинар СЭВ ТФ-78, т.2, Будапешт, 1978.

3. АЭС «Куданкулам». Система контроля, управления и диагностики СКУД РУ В-412. Система внутриреакторного контроля СВРК-01Р. Программа функционирования ПТК-З. Частное техническое задание ТЗ 09.1101-2002,ООО «СНИИП-АСКУР», 2009.

Размещено на Allbest.ru

...