CC BY
38
ЭКСПЛУАТАЦИЯ ОБЪЕКТОВ АТОМНОЙ ОТРАСЛИ
УДК 621.039.51
ОЦЕНКА ТОЧНОСТИ МЕТОДА ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТОВ РЕАКТИВНОСТИ И ЭФФЕКТИВНОСТИ ТВЕРДЫХ ПОГЛОТИТЕЛЕЙ НА ОСНОВЕ РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧИ ДИНАМИКИ ТОЧЕЧНОГО РЕАКТОРА
© 2018 В.К. Семенов, М.А. Вольман
Ивановский государственный энергетический университет имени В.И. Ленина, Иваново, Россия
На основе сравнения результатов решения прямой и обратной задачи динамики точечного ядерного реактора произведена оценка точности методики определения коэффициентов реактивности и эффективности твердых поглотителей, апробированной на полномасштабном компьютерном имитаторе энергоблока с реактором ВВЭР-1000. Показано, что оптимизационный метод решения обратной задачи динамики реактора обладает хорошей точностью.
Ключевые слова: эффекты и коэффициенты реактивности, обращенные решения, динамика реактора, численные эксперименты, точечный реактор, метод оптимизации.
Поступила в редакцию 25.06.2018 г.
В процессе эксплуатации ядерного реактора производится экспериментальное уточнение его нейтронно-физических характеристик, в частности, коэффициентов реактивности и эффективности групп твердых поглотителей. Эксперименты проводятся на основе стандартной методики, суть которой состоит в следующем [1, 2]. При постоянном отравлении реактора ксеноном и самарием и неизменной концентрации борной кислоты изменение реактивности реактора при малых возмущениях определяется выражением:
(ОР)ДТВ + (|Р)А3 ДМ + Адр, + (|Р)ДИ = 0 . (1)
ОТ Ж Ор, ОИ
Ор Ор Ор „
где — = — н---суммарный температурный коэффициент реактивности по
ОТ ОТ в 0Ти
температуре топлива Тв и теплоносителя Ти;
(—Р)аз =—— — - мощностной коэффициент реактивности при постоянной ОМ 0Ти №
температуре теплоносителя в активной зоне; N - нейтронная мощность реактора;
к - эффективный коэффициент теплопередачи от топлива к теплоносителю; Б - эффективная поверхность теплообмена;
—— барометрический коэффициент реактивности; 0р,
р1 - давление теплоносителя в первом контуре;
- дифференциальная эффективность группы органов регулирования;
0И
И - положение регулирующей группы.
© Национальный исследовательский ядерный университет «МИФИ», 2018
Рассматриваются отдельно возмущение реактора изменением температуры теплоносителя на входе в реактор ДТв при постоянных р1 и Н, возмущение изменением давления в первом контуре при неизменном положении регулирующей группы и возмущение изменением положения регулирующей группы при постоянстве давления в первом контуре. В результате получается система трех уравнений с четырьмя неизвестными. Для замыкания системы уравнений выполняется отдельный
эксперимент по определению коэффициента . Осуществить указанные
ЭН
эксперименты при варьировании одних параметров и поддержании других неизменными весьма непросто.
Нами предлагается другая методика определения вышеназванных характеристик, основанная на решении обратной задачи динамики реактора методом оптимизации. [3, 4] Решение обратной задачи динамики реактора позволяет найти входные параметры модели по имеющимся выходным динамическим характеристикам, полученным экспериментально. Примером успешной реализации подобного подхода является известное обращенное решение уравнений кинетики реактора, позволяющее построить реактиметр - прибор, с помощью которого рассчитывается реактивность реактора. Эту математическую модель с определенными мерами и ограничениями удалось применить к реальному объекту [5-7]. Обращенное решение уравнений кинетики нашло применение и в расчётах эффективности аварийной защиты реакторов ВВЭР-1000 [8].
Оценку точности метода проведем на основе численного эксперимента, основанного на решении прямой задачи с известными коэффициентами реактивности и известным законом возмущения реактора по реактивности. На основе полученных решений методом оптимизации решим обратную задачу по нахождению значений введенной реактивности и коэффициентов реактивности. Сравнение полученных значений с входными параметрами этих величин, заданными в прямой задаче, позволит оценить точность предлагаемой методики.
В качестве объекта моделирования будем иметь в виду реактор ВВЭР-1000. Для реактора с сосредоточенными параметрами (точечный реактор) математическая модель может быть представлена следующей системой нелинейных жестких дифференциальных уравнений [5]:
аК_ 5р(Х) + а (Тц - Ти (0)) + а 2 (Тв - Тв (0)) - р АТ | 6
- _ !> + Ь (2)
ах х т р!
-г1 -—К-Х^, ] = 1, ..., 6, (3)
1 ГТ-1
тцСц-Т1 - N - кБ(Тц - Тв ), (4)
ах
1гр
т в С в - кБ(Тц - Тв ) - у вСв 0(ТвЫХ - ТвХ ). (5)
ах
где N - нейтронная мощность реактора;
N - мощность, вносимая]-ой группой запаздывающих нейтронов (] = 1, .,6); Р] и ^ - соответственно доля запаздывающих нейтронов и постоянная времени
распада ядер-предшественников ]-ой группы; Р - X Р] - суммарная доля запаздывающих нейтронов;
тт - приведенное время жизни одного поколения мгновенных нейтронов;
ти и си - соответственно масса и удельная теплоемкость топливной загрузки;
к, Б, Ти - соответственно коэффициент теплоотдачи, поверхность теплоотдачи и
температура топлива; Тв, Твых и Твх - соответственно средняя температура теплоносителя, температура на
выходе и входе в реактор; ув, Св, О - соответственно плотность, удельная теплоемкость теплоносителя
при постоянном давлении и расход теплоносителя; 5р(Х) - внешнее возмущение реактора по реактивности;
а1 и а2 - соответственно температурные коэффициенты реактивности по температуре топлива и теплоносителя.
Поскольку уравнения записаны для точечного реактора, то, тем самым, предполагается, что за время переходного процесса пространственное распределение нейтронного поля реактора не успевает заметно измениться. Во-вторых, также считается, что за время протекания вышеназванных процессов не изменяются концентрация борной кислоты и отравление реактора ксеноном и самарием. В последнем уравнении необходимо указать связь между указанными тремя температурами теплоносителя. Как показал проведенный нами анализ, с достаточной степенью точности можно воспользоваться соотношением [9]:
Т + Т
Тв = 'вых ^ ' вх . (6)
Дополняя уравнения (2-6) соответствующими начальными условиями:
в ¡N0
N = N0, N. = -¡-0-, Ти = Ти(0), Тв(0), Твх = Твх(0) и бр(1), (7)
можно найти зависимости всех динамических величин от времени при заданном законе изменения реактивности реактора 5р(1:).
В качестве примера на рисунках 1-5 приведены расчетные кривые при переводе реактора на более низкий уровень мощности за счет введения отрицательной реактивности брО) = - 1,21 • 10-3-(1 - ехр (-0,51) за время ^ = 20 с, имитирующей действие твердых поглотителей. Принятые для расчета значения коэффициентов реактивности имели следующие значения: ах = - 1,22 • 1СГ5 1/К и а2 = - 1,00 • Ю-4 1/К.
Рисунок 1 - Возмущение реактора введенной реактивностью [The perturbation of the reactor by the
introduced reactivity]
Рисунок 2 - Зависимость мощности реактора от времени на начальной стадии: 1 - нейтронная мощность; 2 - тепловая мощность [Dependence of reactor power on time at the initial stage:
1 - neutron power; 2 - thermal power]
Рисунок 3 - Зависимость температуры топлива от временинав начальной стадии [Fuel temperature versus
time in the initial stage]
Рисунок 4 - Зависимость температуры теплоносителя от времени на начальной стадии [Temperature dependence of the coolant temperature at the initial stage]
Рисунок 5 - Зависимость реактивности реактора от времени на начальной стадии процесса [Dependence of reactivity of the reactor on time at the initial stage of the process]
Рисунок 6 - Зависимость реактивности реактора от времени после окончания ввода возмущения [Dependence of reactivity of the reactor on time at the initial stage of the process]
Далее по полученным результатам решим обратную задачу. Вначале найдем максимальное значение введенной отрицательной реактивности. Так как реактор переводится из одного стационарного состояния в другое стационарное состояние, то введенная максимальная реактивность 5p(t1) компенсируется температурными эффектами реактивности:
5p(t1) = a1(Tu(0) - ВД)+ а2(Тв(0) - Тв^)). (8)
где tR - время окончания переходного процесса, равное в нашем примере 400 с.
Зависимость полной реактивности реактора от времени, от момента t1 до окончания переходного процесса tk, определяется выражением:
Ap(t) = a1ATu(t) + а2ЛТв(0, (9)
где ATu(t) = Tu(t) - Tu(U ЛТв(0 = Твф - Тв(^.
Решая задачу оптимизации, составим функционал:
1(аь а2) = | [Др(1) - а1ДТи(1) - а2 ДТВ (1))]2ё1, (10)
минимизируя который по параметрам а, получим систему двух алгебраических уравнений, определяющих искомые коэффициенты реактивности:
I к
| [Др(1) - а: ДТи (1) - а2 ДТв (1) ]ДТв (1)ё1 = 0, (11)
1 к
| [Др(1) - а1 ДТи (1) - а2 ДТв (1)]ДТи (1)ё1 = 0 . (12)
Вычислив соответствующие интегралы, и решив уравнения (11) и (12), получили значения а1 = - 1,22 • 10-5 1/К и а2 = - 1,00 • 10-4 1/К, что полностью совпадает со входными параметрами. Подставляя найденные значения коэффициентов реактивности и рассчитанные температуры в формулу (8), получим 5рмакс = 5р(1) = - 1,209 • 10 , что также с хорошей точностью согласуется со входной реактивностью - 1,21 • 10 . Численные эксперименты, выполненные для других значений коэффициентов реактивности и значений 5р(1) также показывают хорошее согласие результатов решения прямой и обратной задачи.
Различного рода тренажеры в настоящее время нашли широчайшее применение [10-14], причем они могут использоваться не только для обучения персонала существующим методикам физических экспериментов, но и для разработки и апробации новых методик измерений до выхода на реальные объекты. Апробация предлагаемой методики выполнена на полномасштабном аналитическом тренажере атомного блока с реактором ВВЭР-1000, являющегося программным продуктом ООО «Вестерн Сервисес» [15], реализованном на персональном компьютере в программной среде 3КеуМаБ1ег и являющегося компьютерной моделью функционирования энергоблока в различных режимах в реальном времени.
При симуляции экспериментов на тренажере при отключенном автоматическом регуляторе мощности возмущение реактора проводилось регулирующей группой на фиксируемую глубину со штатной скоростью 2 см/с. Предварительно за счет борной перекомпенсации регулирующая группа выводилась на верхний конечный выключатель. В дальнейшем концентрация борной кислоты оставалась неизменной, отравление реактора ксеноном за время проведения численных экспериментов не изменялось, и для «чистоты» эксперимента рассматривалась модель неотравленного реактора. В соотношения (8) и (9) были добавлены слагаемые, учитывающие положительный барометрический эффект реактивности. В ходе переходного процесса измерялись интересующие нас временные зависимости: положение регулирующей группы нейтронная мощность, реактивность реактора, температуры теплоносителя и давление теплоносителя в первом контуре. Так как зависимость температуры топлива от времени непосредственно на эксперименте не определяется, то её можно найти из уравнения теплового баланса (4), решение которого выражается в квадратурах:
1 1 1 — 1' — 1' --е Ти 1 --е Ти 1 --/-10Ч
Ти(1) = Ти(0)е ти + ^-1N(1 ')е^1'+ е-1Тв(1')е^1'. ( )
тиси 0 Ти 0
где ти = —= 3,6 + 4,2 с - характерное время теплопередачи от топлива к кБ
теплоносителю.
В остальном методика обработки результатов экспериментов оставалась в соответствии с вышеизложенным.
Таким образом, оптимизационный метод решения обратной задачи позволил за одну серию экспериментов определить коэффициенты реактивности и введенную реактивность, что позволяет для поэтапного возмущения реактора погружением твердых поглотителей построить их интегральную и дифференциальную характеристики.
СПИСОК ЛИТЕРАТУРЫ
1. РД ЭО 0151-2004. Методики расчета нейтронно-физических характеристик по данным физических экспериментов на энергоблоках атомных электростанций с реакторами ВВ ЭР-1000 [Текст]. - Москва : Росэнергоатом, 2005. - 101 с.
2. РД ЭО 0150-2004. Типовые программы и методики проведения физических экспериментов на энергоблоках атомных электростанций с реакторами ВВЭР-1000 [Текст]. - Москва : Росэнергоатом, 2005. - 273 с.
3. Вольман, М.А. Имитационное моделирование нейтронно-физических и теплогидравлических процессов в реакторах ВВЭР-1000: автореф. дис. канд. техн. наук [Текст] / М.А. Вольман. -Москва, 2017. - 20 с.
4. Семенов, В.К. Методика определения коэффициентов реактивности и эффективности групп твердых поглотителей на аналитическом тренажере энергоблока ВВЭР-1000 [Текст] /
B.К. Семенов, М.А. Вольман, А.А. Беляков // Вестник Ивановского государственного энергетического университета. - 2018. - № 1. - С. 19-24.
5. Казанский, Ю.А. Кинетика ядерных реакторов. Коэффициенты реактивности. Введение в динамику: учебное пособие [Текст] / Ю.А. Казанский, Я.В. Слекеничс - Москва : МИФИ, 2012. - 300 с.
6. Кошелев, А.С. О возможности создания специализированного ОРУК-реактиметра с токовым детектором нейтронов [Текст] / А.С. Кошелев, А.В. Арапов, М.А. Овчинников // Вопросы атомной науки и техники. Серия : Физика ядерных реакторов. - 2015. - № 4. - С. 39-48.
7. Колесов, В.Ф. Истоки неточностей в реактивности, определяемой с помощью обращенного решения уравнений кинетики [Текст] // Вопросы атомной науки и техники. Серия : Физика ядерных реакторов. - 2013. - № 3. - С. 30-45.
8. Зизин, М.Н. О трактовке обращенного уравнения кинетики и пространственно-временных расчётов эффективности аварийной защиты ВВЭР-1000 [Текст] / М.Н. Зизин, Л.Д. Иванов // Вопросы атомной науки и техники. Серия : Физика ядерных реакторов. - 2012. - № 2. -
C. 28-43.
9. Семенов, В.К. Обоснование математической модели теплообмена для реактора с сосредоточенными параметрами [Текст] / В.К. Семенов, М.А. Вольман // Глобальная ядерная безопасность. - 2015. - № 4 (17). - С. 35-42.
10.Дугинов, О.Б. SSL DYNCO LAB SYSTEM - Готовые решения в области обучающих тренажерных систем [Текст] / О.Б. Дугинов, А.В. Левченко, Д.С. Самохин // XIII Международная конференция «Безопасность АЭС и подготовка кадров». - Обнинск : ИАТЭ, 2013 - С. 195-198.
11. Чернов, Е.В. Применение компьютерных тренажеров в подготовке специалистов для ядерной энергетики [Текст] / Е.В. Чернов, С.Б. Выговский, Д.И. Макаун, В.Е. Ямный // Энергетическая Стратегия. - 2010. - № 5 (17). - С. 58-60.
12. Лощаков, И.И. Компьютерный тренажер «АЭС с РБМК-1000» [Текст] / И.И. Лощаков, К.М. Мазурик, Г.А. Ромахова, Г.П. Шаргин // Безопасность АЭС и подготовка кадров. Х Международная конференция : Тезисы докладов (Обнинск, 1-4 октября 2007 г.). Часть 1. -Обнинск: ИАТЭ, 2007. - С. 57-58.
13. Росляков, М.В. Опыт внедрения тренажера оперативного моделирования аварийных ситуаций «ТОМАС- 1А» в учебно-методический процесс кафедры «Атомная энергетика» УГТУ-УПИ [Текст] / М.В. Росляков, Г.П. Титов, С.Е. Щеклеин, А.И. Айзатулин, Е.Ф. Селезнев, И.П. Федоров // Перспективные энергетические технологии. Экология. Экономика, безопасность и подготовка кадров. Сбор. науч. трудов. - Екатеринбург, 2006. - С. 150.
14. Чернаков, В.А. Особенности современных моделирующих комплексов сложных технологических объектов (на примере анализатора режимов АЭС с ВВЭР) [Текст] / В.А. Чернаков, М.А. Осадчий // Приборы. - № 7. - 2002. - С. 12.
15. Western Services Corporation [Электронный ресурс]. - URL: http://www.ws-corp.com.
REFERENCES
[1] RD EO 0151-2004. Metodiki rascheta neytronno-fizicheskikh kharakteristik po dannym fizicheskikh eksperimentov na energoblokakh atomnykh elektrostantsiy s reaktorami VVER-1000 [Calculation Methods of Neutron-Physical Characteristics According to Physical Experiments on Nuclear Power Units with WWER-1000]. Moscow: Rosenergoatom. 2005. 101 p. (in Russian).
[2] RD EO 0150-2004. Tipovye programmy i metodiki provedeniya fizicheskikh eksperimentov na energoblokakh atomnykh elektrostantsiy s reaktorami VVER-1000 [Typical Programs and Methods of Physical Experiments on Nuclear Power Units with WWER-1000]. Moscow: Rosenergoatom. 2005. 273 p. (in Russian).
[3] Volman M.A. Imitacionnoe modelirovanie nejtronno-fizicheskih i teplogidravlicheskih processov v reaktorah VVER-1000 [Simulation Modeling of Neutron-Physical and Thermal-Hydraulic Processes in WWER-1000 Reactors. Thesis Abstract of Ph. D. in Engineering]. Moscow. 2017. 20 p. (in Russian).
[4] Semenov V.K., Volman M.A., Belyakov A.A. Metodika opredeleniya koeffitsientov reaktivnosti i effektivnosti grupp tverdykh poglotiteley na analiticheskom trenazhere energobloka VVER-1000 [Method for Determination Of Reactivity Coefficients and Efficiency of Solid Absorbers on Computer Simulator WWER-1000]. Vestnik Ivanovskogo gosudarstvennogo energeticheskogo universiteta [Bulletin of Ivanovo State Power University]. 2018. issue 1. Р. 19-24 (in Russian).
[5] Kazanskij Yu.A., Slekenichs Ya.V. Kinetika yadernykh reaktorov. Koeffitsienty reaktivnosti. Vvedenie v dinamiku: Uchebnoe posobie [The Coefficients of Reactivity. Introduction to Dynamics: Training Course]. Moscow: National Research Nuclear University MEPhI. 2012. 300 p. (in Russian).
[6] Koshelev A.S., Arapov A.V., Ovchinnikov M.A. O vozmozhnosti sozdaniya spetsializirovannogo oruk-reaktimetra s tokovym detektorom neytronov [Possibility of Creation of a Specialized ORUK-Reactimeter with a Current Detector of Neutrons]. Voprosy atomnoy nauki i tekhniki. Seriya: Fizika yadernykh reaktorov [Questions of Atomic Science and Technology. Series: Physics of Nuclear Reactors]. 2015. issue 4. Р. 39-48 (in Russian).
[7] Kolesov V.F. Istoki netochnostey v reaktivnosti, opredelyaemoy s pomoshch'yu obrashchennogo resheniya uravneniy kinetiki [Background of Uncertainties in Reactivity Determined with the Aid of Inverse Solution of Kinetics Equations]. Voprosy atomnoy nauki i tekhniki. Seriya: Fizika yadernykh reaktorov [Questions of Atomic Science and Technology. Series: Physics of Nuclear Reactors]. 2013. issue 3. Р. 30-45 (in Russian).
[8] Zizin M.N., Ivanov L.D. O traktovke obrashchennogo uravneniya kinetiki i prostranstvenno-vremennykh raschetov effektivnosti avariynoy zashchity VVER-1000 [About Interpretation of the Inverted Equation of Kinetics and Space-Time Calculations of the VVER-1000 Reactor Scram System Effectiveness]. Voprosy atomnoy nauki i tekhniki. Seriya: Fizika yadernykh reaktorov [Questions of Atomic Science and Technology. Series: Physics of Nuclear Reactors]. 2012. issue 2. Р. 28-43 (in Russian).
[9] Semenov V.K., Volman M.A. Obosnovanie matematicheskoy modeli teploobmena dlya reaktora s sosredotochennymi parametrami [Argumentation of Heat-Exchange Mathematical Model for Reactor with Lumped Parameters]. Global'naya yadernaya bezopasnost' [Global Nuclear Safety]. 2015. issue 4 (17). Р. 35-42 (in Russian).
[10] Duginov O.B., Levchenko A.V., Samokhin D.S. SSL DYNCO LAB SYSTEM - Gotovye resheniya v oblasti obuchayushchikh trenazhernykh system [SSL DYNCO LAB SYSTEM - Ready-Made Solutions in the Field of Training Simulator Systems]. XIII Mezhdunarodnaya konferentsiya «Bezopasnost' AES i podgotovka kadrov» [XIII International Conference "NPP Safety and Training"]. Obninsk: IATE. 2013. Р. 195-198 (in Russian).
[11] Chernov E.V., Vygovskiy S.B., Makaun D.I., Yamnyy V.E. Primenenie komp'yuternykh trenazherov v podgotovke spetsialistov dlya yadernoy energetiki [The Use of Computer Simulators in Training Specialists for Nuclear Power]. Energeticheskaya Strategiya (Energy Strategy). 2010. issue 5 (17). Р.58-60 (in Russian).
[12] Loshchakov I.I., Mazurik K.M., Romakhova G.A., Shargin G.P. Komp'yuternyy trenazher «AES s RBMK-1000» [Computer Simulator "NPP with RBMK-1000"]. X Mezhdunarodnaya konferentsiya «Bezopasnost' AES i podgotovka kadrov» [Safety of Nuclear Power Plants and Training of Personnel. X International Conference: Abstracts]. Obninsk: IATE. 2007. Р. 57-58 (in Russian).
[13] Roslyakov M.V., Titov G.P., Shcheklein S.E., Ayzatulin A.I., Seleznev E.F., Fedorov I.P. Opyt vnedreniya trenazhera operativnogo modelirovaniya avariynykh situatsiy «TOMAS-1A» v uchebno-metodicheskiy protsess kafedry «Atomnaya energetika» UGTU-UPI [Experience in Implementing the Simulator of Operational Modeling of Emergency Situations "ТОМАС-1А" in the Educational-
Methodical Process]. Perspektivnye energeticheskie tekhnologii. Ekologiya. Ekonomika, bezopasnost' i podgotovka kadrov ["Perspective Energy Technologies. Ecology. Economics, Security and Training ". Collection of scientific papers]. Ekaterinburg. 2006. P. 150 (in Russian).
[14] Chernakov V.A., Osadchiy M.A. Osobennosti sovremennykh modeliruyushchikh kompleksov slozhnykh tekhnologicheskikh ob"ektov (na primere analizatora rezhimov AES s VVER) [Features of Modern Modeling Complexes of Complex Technological Objects]. Pribory [Devices]. 2002. issue 7. P. 12 (in Russian).
[15] Western Services Corporation. Available at: http://www.ws-corp.com.
Estimation of the Accuracy of the Method of Reactivity Coefficient Identification and Solid Absorbers Efficiency Based on Reverse Problem Solution of Point Reactor
Dynamics
V.K. Semenov1, M.A. Volman2
Ivanovo State Power Engineering University, 34, Rabfakovskaya st., Ivanovo, Russia 153003
'ORCID iD: 0000-0001-9765-2096 WoS Researcher ID: T-1066-2017
e-mail: Semenov [email protected] 2ORCID iD: 0000-0001-6805-6287 WoS Researcher ID: M-1304-201 e-mail: [email protected]
Abstract - The estimation of the method of the reactivity coefficient determining and the absorbers efficiency is evaluated. The results of solving the direct and reverse task of the dynamics of a point reactor are compared. The technique is tested on a full-scale computer simulator of a power unit with a WWER-1000 reactor. It is shown that the optimization method for solving the inverse problem of reactor dynamics has good accuracy.
Keywords: effects and coefficients reactivity, reversed solutions, reactor dynamics, numerical experiments, point reactor, optimization method.
