WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Александрович расчетный анализ методов измерения коэффициентов реактивности рбмк.

РОССИЙСКИЙ НАУЧНЫЙ ЦЕНТР «КУРЧАТОВСКИЙ ИНСТИТУТ»

На правах рукописи.

Балыгин Алексей Александрович

Расчетный анализ методов измерения коэффициентов реактивности РБМК.

Специальность 05.14.03 Ядерные энергетические установки включая проектирование, эксплуатацию и вывод из эксплуатации.

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук.

Москва.

2008.


Работа выполнена в ФГУ Российский научный центр «Курчатовский Институт»

Научный руководитель: доктор технических наук

Краюшкин Александр Викторович

ИЯР ФГУ РНЦ «Курчатовский институт»

Официальные оппоненты: доктор технических наук

Цибульский Виктор Филлипович

ИЯР ФГУ РНЦ «Курчатовский институт»

кандидат физико-математических наук

Попыкин Александр Иванович

НТЦ ЯРБ

Ведущая организация «Всероссийский Научно- Исследовательский Институт по эксплуатации атомных электростанций» (ОАО «ВНИИАЭС»)

Защита состоится « ___» _____________ 2009 г. в ___ ч. ___ мин.

на заседании диссертационного совета Д 320.009.06 в Российском научном центре «Курчатовский институт» по адресу 123182, г. Москва, пл. Курчатова, д1.

С диссертацией можно ознакомиться в библиотеке РНЦ «Курчатовский институт»

Автореферат разослан «____» _____________ 2009. г.

Ученый секретарь диссертационного совета,

доктор технических наук, профессор В.Г. Мадеев.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ.

Актуальность работы.

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

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

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

Цели и задачи работы.

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

2. Обеспечение корректного использования результатов измерений при валидации расчетных кодов.

3. Разработка новых методов измерений коэффициентов реактивности и обработки результатов измерений.

Результаты работы, выносимые на защиту.

  1. Методика и результаты стационарных и динамических расчетов для анализа причин различия величин коэффициентов реактивности, получаемых в эксперименте и стационарном расчете.
  2. Комплексный метод определения быстрого мощностного и графитового температурного коэффициентов реактивности.
  3. Пассивный метод определения графитового коэффициента реактивности и эффекта ксенонового отравления.
  4. Новая методика обработки результатов измерения парового коэффициента реактивности.


Методы исследования.

Стационарные и динамические расчеты проводились с использованием различных версий программы STEPAN.

Научная новизна.

Проанализированы методические расхождения в определении величин коэффициентов реактивности, полученных экспериментальным путем и из расчета по стационарной версии программы STEPAN:

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

- оценено влияние «медленных» обратных связей по температуре графита и концентрации ксенона-135 на результаты измерений парового () и быстрого мощностного (w) коэффициентов реактивности.

- дано расчетное обоснование «обратного хода» реактивности в процессе эксперимента по измерению.

Разработаны предложения по:

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

- пассивному методу определения графитового температурного коэффициента реактивности (с);

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

Достоверность.

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

Практическая значимость работы.

  • Проведенный анализ методик измерения коэффициентов реактивности обеспечивает корректное использование результатов измерений при валидации программ нейтронно-теплогидравлического расчета, способствуя повышению качества расчетного прогнозирования при обосновании мероприятий по повышению безопасности РБМК.
  • Объяснение «обратного хода» реактивности при измерениях дает лучшее понимание физики реактора.
  • Разработанные новый способ обработки результатов измерений и пассивный метод измерения графитового коэффициента реактивности позволяют сократить время, затрачиваемое на измерения на АЭС. Решением Совета по ядерной безопасности ОАО «Концерн Энергоатом» данные разработки рекомендованы для включения в «Комплексную методику определения физических и динамических характеристик реакторов РБМК-1000».

Апробация работы.

Материалы, представленные в диссертации были доложены на следующих научных конференциях:

1. XII семинар по проблемам физики реакторов. «Волга - 2002» Москва, МИФИ, 2-6 сентября 2002 г.

  1. IX Международная Конференция «Безопасность и подготовка кадров». Обнинск, 24-28 октября 2005 г.
  2. Пятая Международная научно-техническая конференция «Безопасность, экономика и эффективность атомной энергетики МНТК 2006» Москва, Росэнергоатом, 19-21 апреля 2006.
  3. Шестая Международная научно-техническая конференция «Безопасность, экономика и эффективность атомной энергетики МНТК 2008» Москва, Росэнергоатом, 21-23 мая 2008.
  4. XV семинар по проблемам физики реакторов «Волга - 2008», Москва, МИФИ, 2 – 6 сентября 2008 г.

Личный вклад автора.

    • Проведены динамические расчеты, моделирующие эксперименты по измерению коэффициентов реактивности РБМК.
    • Выполнена оценка влияния различных факторов, не учитываемых в стационарных расчетах, на величины коэффициентов реактивности.
    • Получено объяснение «обратного хода» реактивности при измерениях парового коэффициента реактивности.
    • Разработаны новые методики экспериментального определения быстрого мощностного, парового и графитового температурного коэффициентов реактивности.

Публикации.

Список основных публикаций приведен в конце автореферата.

Структура и объем диссертации.

Материал диссертационной работы изложен на 164 страницах, содержит список литературы из 63 наименований, 41 таблицу и 76 рисунков.

СОДЕРЖАНИЕ РАБОТЫ.

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

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

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

В диссертации приводятся результаты расчетного моделирования экспериментов по программе STEPAN/KOBRA. На основании полученных результатов анализируются возможные причины расхождения экспериментальных и расчетных результатов.

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

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

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

В стационарном расчете получаются ”чистые” значения величин и w, не учитывающие влияния изменения таких параметров как температура графита и концентрация ксенона. А в эксперименте изменение температуры графита и концентрации ксенона могут оказывать влияние на измеряемые значения и w.

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

Таким образом, в эксперименте непросто правильно определить «теоретические» величины коэффициентов реактивности и w. Эти величины «в чистом виде» можно получить только в результате стационарного расчета. При этом критерием правильности стационарного расчета может служить совпадение измеренных и полученных в результате расчетного моделирования измерений величин и w, при условии, что стационарный расчет и динамическое моделирование эксперимента проводились в рамках одной модели.

Численное моделирование экспериментов по измерению и w.

Приводится сравнение величин и w, полученных экспериментально, в результате расчетного моделирования экспериментов и в результате стационарных расчетов. Полученные результаты представлены в таблицах 1 и 2.

Таблица 1.
Величины w, полученные из эксперимента, стационарного расчета и динамического расчета.

Блок, дата w, 10-4/МВт
wст wврд wбик wэкс
ЛАЭС-1, 30.01.01 -3.4(62) -2.6(23) -2.3(9) -2.1
ЛАЭС-2, 22.08.03 -4.4(52) -3.1(7) -2.8(-3) -2.9
ЛАЭС-4, 18.11.02 -3.9(70) -2.88(25) -2.2(-4) -2.3
КАЭС-1, 24.01.04 -3.3(14) -2.4(-17) -2.1(-28) -2.9
КАЭС-4, 07.04.03 -3.7(37) -2.6(-3) -2.05(-24) -2.7
САЭС-1, 20.01.04 -3.1(26) -2.02(-17) - -2.48


Таблица 2.
Величины , полученные из эксперимента, стационарного расчета и динамического расчета.

Блок, дата , 10-2
ст дин экс
ЛАЭС-1, 30.01.01 0.56 (-31) 0.79 (-4) 0.82
ЛАЭС-2, 22.08.03 0.01 (-97) 0.37 (12) 0.33
ЛАЭС-4, 18.11.02 0.27 (-64) 0.47 (-37) 0.75
КАЭС-1, 21.01.04 0.74 (-27) 1.55 (51) 1.02
КАЭС-4, 07.04.03 0.70 (-1) 0.88 (24) 0.71
САЭС-1, 20.01.04 1.03 (30) 1.34 (70) 0.79

В первых столбцах таблиц 1 и 2 представлены номера блоков и даты состояний реактора, для которых проводились эксперименты и расчеты.

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

1 – расчет w, полученный с учетом величины введенной реактивности, «измеренной» внутризонными датчиками (врд);

2 - расчет w, полученный с учетом величины введенной реактивности, «измеренной» боковыми ионизационными датчиками (бик);

В скобках приводится отличие в процентах абсолютной расчетной величины w от экспериментальной величины wэкс.

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

Влияние пространственной неоднородности нейтронно-физических свойств на измеряемые величины и w.

Величины w, полученные из статического расчета, систематически больше по абсолютной величине, чем величины w, полученные экспериментально и в результате динамического моделирования эксперимента. Это расхождение связано с тем, что в статическом расчете w мощность меняется равномерно (увеличивается на 4% в каждой расчетной ноде). В эксперименте и в динамическом расчете мощность меняется в результате частичного погружения в активную зону 4 стержней. Это приводит к перераспределению энерговыделения по объему активной зоны.

Для оценки влияния перераспределения энерговыделения по объему активной зоны проводились сравнительные расчеты, моделирующие измерение w на 2 блоке ЛАЭС по состоянию на 22 августа 2003 года. При этом рассматривалось 2 варианта:

  1. Реактивность вносилась перемещением 4 стержней, расположенных в центральной части активной зоны. Начальная глубина погружения стержней около 300 см.
  2. Реактивность вносилась одновременным перемещением всех 136 стержней РР на глубину 4 см со своего начального положения. В этом варианте обеспечивается более равномерное внесение возмущения по объему активной зоны.

Были получены следующие результаты:

первый вариант - w=-3.14·10-4 /МВт;

второй вариант - w=-3.92·10-4 /МВт.

Величина w, полученная при одновременном погружении всех стержней РР (2 вариант) на 25% выше по абсолютной величине, чем величина w, полученная для варианта 1. Таким образом, при более равномерном вводе реактивности по объему активной зоны измеренная величина w была бы больше по абсолютной величине. Величина w, полученная во втором варианте, лучше согласуется с результатами статических расчетов, где реактивность вводится равномерно по объему активной зоны.

Использование датчиков БИК, расположенных за пределами реактора для измерения введенной реактивности при взвешивании стержней приводит к занижению значения введенной реактивности и абсолютной величины w по сравнению с реактивностью, полученной с учетом изменения нейтронного потока во всей активной зоне. Величина занижения для разных вариантов может меняться от 12 до 30% (см. табл. 1).

Анализируется различие расчетных величин и w, полученных на основании моделирования экспериментов с использованием для оценки изменения суммарной мощности реактора боковых датчиков, расположенных за пределами активной зоны (датчики БИК) и внутризонных датчиков (датчики ВРД). Для моделирования экспериментов было выбрано состояние 4 блока Курской АЭС на 7 апреля 2003 года. На рисунке 1 представлены рассчитанные временные зависимости изменения мощности реактора при моделировании измерений w.

Рис 1. Расчетное изменение мощности при моделировании эксперимента.

Сплошной линией обозначена зависимость интегральной мощности реактора, полученная на основании расчета. Интегральная мощность практически совпадает с мощностью, измеренной одной или несколькими группами внутризонных датчиков ВРД. Символами обозначается зависимость мощности реактора, полученная на основании сигналов БИК. Видно, что мощность реактора, восстановленная по показаниям БИК, меняется медленнее, чем реальная мощность реактора. Это вызвано тем, что внереакторные датчики БИК расположены далеко от погружаемых стержней.

На рис. 2 представлена зависимость изменения мощности реактора при моделировании измерений.

Рис. 2. Изменение мощности при моделировании эксперимента по измерению .

Представлены две зависимости мощности от времени: интегральная мощность реактора, полученная прямым расчетом и мощность реактора, полученная с учетом изменения суммарного сигнала 4 датчиков БИК. Видно, что мощность реактора, восстановленная по суммарному сигналу четырех БИК, меняется больше чем мощность реактора, полученная по показаниям внутризонных датчиков.

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

На основании приведенных расчетов показано, что использование внереакторных датчиков БИК для определения уровня мощности в процессе экспериментов по измерению величин и w приводит к тому, что абсолютная величина w занижается примерно на 10%, а величина завышается на 5080%.

Исследуется влияние изменения формы аксиального распределения энерговыделения на величины и w. По стационарной версии программы STEPAN были проведены расчеты для двух состояний 1 и 3 блоков САЭС. При этом вместо реальных показаний 12 аксиальных датчиков, по которым в программе STEPAN восстанавливается высотное поле энерговыделения задавались искусственные для формирования нового поля энерговыделения. Таким образом, было сформировано по 4 искусственных высотных поля энерговыделения для двух рассматриваемых состояний (варианты 14).

Вариант 1 имитирует одногорбое поле с максимумом в центре.

Вариант 2 имитирует одногорбое поле с максимумом, сдвинутым к верху активной зоны.

Вариант 3 имитирует одногорбое поле с максимумом, сдвинутым к низу активной зоны.

Вариант 4 имитирует двугорбое поле энерговыделения.

В таблице 3 представлены значения и w для четырех вариантов для состояний Смоленской АЭС 1 и 2 блоков.

Таблица 3.

Величины и w для 4 искусственных аксиальных распределений энерговыделения.

Состояние Вариант
1 2 3 4
САЭС-1 06.03.06 , 0.83 0.56 0.87 0.75
w, /МВт -3.310-4 -3.710-4 -3.310-4 -3.410-4
САЭС-3 13.09.06 , 0.33 0.12 0.34 0.27
w, /МВт -3.710-4 -4.110-4 -3.710-4 -3.810-4

Согласно представленным результатам расчетов, изменение аксиального профиля энерговыделения может оказывать существенное влияние на величину парового коэффициента реактивности. Измеряемая величина может меняться в пределах 0.20.3 для одного и того же блока без существенного изменения загрузки, уровня мощности и других параметров реактора. Показано, что смещение аксиального поля энерговыделения вверх приводит к уменьшению величин и w. Известно, что при изменении температуры теплоносителя на входе в активную зону паросодержание () меняется по высоте неравномерно и имеет выраженный максимум на расстоянии около 2.5 м от низа активной зоны (на границе экономайзерного и испарительного участков). Смещение поля энерговыделения вверх приводит к тому, что максимум энерговыделения сдвигается в область меньшего, что и приводит к уменьшению и w. Приведенные результаты подтверждают известную корреляцию между величинами и w : w·104=( - 3.0)±0.3 /МВт.

Исследуется влияние температуры теплоносителя на входе в активную зону на величину парового коэффициента реактивности. Согласно технологическому регламенту по эксплуатации РБМК-1000 температура теплоносителя на входе в активную зону на номинальном уровне мощности может меняться в пределах 265270 0С. При этом уменьшение температуры воды на 1 градус приводит к увеличению длины экономайзерного участка примерно на 10 см. Таким образом, изменение температуры на входе в активную зону меняет форму зависимости от высоты участка активной зоны. Это может приводить к некоторому изменению величины парового коэффициента реактивности.

Для оценки влияния температуры теплоносителя на входе в активную зону на величину были проведены расчеты по стационарной версии программы STEPAN состояния 3 блока САЭС на 13.09.06 с температурами теплоносителя на входе в активную зону 265 и 270 0С. В таблице 4 представлены результаты расчетов парового коэффициента реактивности для температуры теплоносителя на входе 265 и 270 0С.

Таблица 4.

Величины парового коэффициента реактивности для 4 искусственных аксиальных распределений энерговыделения. Температура воды на входе в активную зону Твх=265 и 270 0С.

Параметр Вариант
1 2 3 4
265, 0.38 0.11 0.43 0.32
270, 0.33 0.12 0.34 0.27
, 0.05 -0.01 0.09 0.05

В строках 2 и 3 в величинах приведены значения парового коэффициента реактивности для температуры теплоносителя на входе в активную зону 265 и 270 0С (265 и 270). В четвертой строке приведена разность этих величин.

Согласно представленным в этом разделе результатам расчетов, при уменьшении температуры теплоносителя на входе в активную зону с 270 до 265 0С, величина увеличивается в пределах 0.1. При уменьшении температуры теплоносителя увеличивается длина экономайзерного участка и максимум смещается выше, в область более высокого энерговыделения. Это и приводит к росту. При этом на величину изменения оказывает влияние аксиальное распределение энерговыделения в активной зоне.

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

В настоящее время в РБМК загружается топливо, содержащее выгорающий поглотитель (Er2O3). Нейтронно-физические свойства такого топлива сильно зависят от глубины выгорания. Согласно результатам расчетов «индивидуальный» такой канала со свежим топливом отрицательный, но по мере выгорания растет и становится положительным.

Обычно каналы максимальной мощности - это каналы с относительно свежим «эрбиевым» топливом, имеющие большой отрицательный эффект обезвоживания. В каналах максимальной мощности устанавливаются большие расходы теплоносителя (6-8 кг/сек). Это приводит к тому, что скорость движения теплоносителя в этих каналах максимальная. Таким образом, каналы с большой мощностью быстрее заполняются «холодным» теплоносителем, имеющим бльшую плотность (при увеличении расхода питательной воды). Это приводит к начальному росту реактивности и мощности реактора.

На рисунке 3 представлена временная зависимость реактивности от времени при измерении на 4 блоке ЛАЭС 18 ноября 2002 года.

Вариант 1 – это стандартный расчет.

Вариант 2 – это расчет того же состояния, но с измененным распределением расхода теплоносителя –6 кг/сек на входе в каждый канал.

Рис. 3. Зависимость реактивности от времени.

Во втором варианте скорость движения теплоносителя во всех каналах одинаковая. Следовательно, заполнение всех каналов «холодной водой» происходит одновременно. Поэтому во втором варианте роста реактивности в начальный момент времени нет.

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

Исследуется влияние изменения температуры графита и концентрации ксенона на результаты измерения w. Предполагается, что изменение этих параметров не оказывает существенного влияния на измеряемую величину w, так как сами измерения проводятся в течение достаточно короткого времени – 20-30 секунд. Для определения роли обратной связи, связанной с изменением этих параметров были проведены расчеты для трех вариантов:

  1. При моделировании эксперимента учитывалось изменение концентрации ксенона и температуры графита.
  2. При моделировании эксперимента учитывалось изменение температуры графита. Концентрация ксенона в топливе предполагалась постоянной и равной начальной величине.
  3. При моделировании эксперимента предполагалось, что температура графита и концентрация ксенона в топливе постоянны и равны начальным величинам.

На рисунке 4 представлены расчетные временные зависимости изменения мощности реактора, полученные при моделировании эксперимента для одного из рассмотренных состояний – 1 блока КАЭС 24 января 2004 года.


Рис. 4 Изменение мощности при моделировании измерения w на 1 блоке КАЭС 24.01.2004.

Видно, что в течение первых 20-30 секунд влияние изменения температуры графита и концентрации ксенона действительно пренебрежимо мало. Как показывают расчеты, величины w, полученные с учетом и без учета влияния изменения температуры графита и концентрации ксенона (варианты 1 и 3) различаются на 1.5-2.0%, что намного меньше точности измерения w. Однако видно также, что в третьем варианте мощность реактора через 30 секунд после начала эксперимента вновь начинает снижаться и стабилизируется на более низком уровне только через 150-200 секунд.

Согласно расчетам реальная величина изменения мощности на 20-25 МВт больше, чем величина, взятая в области локального максимума (как это принято при обработке измерений). Однако к этому моменту времени влияние обратной связи по температуре графита и концентрации ксенона уже настолько значительно, что стабилизация мощности в первом и втором вариантах не достигается. Согласно полученным результатам, использование экспериментально измеренной величины изменения мощности в области локального максимума приводит к завышению абсолютной величины w примерно на 15%.

Исследуется влияние изменения температуры графита и концентрации ксенона на измеряемую величину. Измерение проводится на интервале времени около 120 секунд. Этого времени может оказаться достаточно для того, чтобы начала проявляться обратная связь по температуре графита и концентрации ксенона. Измерение температуры графита и концентрации ксенона в процессе эксперимента может быть связано также с нестационарной мощностью реактора перед началом измерений. Согласно «Комплексной методике» измерение парового коэффициента реактивности должно проводиться на реакторе, в котором мощность поддерживалась на постоянном уровне в течение достаточно долгого времени (не менее суток до начала измерения). Это необходимо для того, чтобы исключить влияние «фонового» изменения мощности на результаты измерения. Однако это требование не всегда соблюдается. Для оценки влияния обратных связей по температуре графита и концентрации ксенона по программе STEPAN/KOBRA были проведены расчеты, моделирующие первые четыре опыта эксперимента по измерению методом малых возмущений. Для расчетов было взято состояние 4 блока ЛАЭС на 18 ноября 2002 года. Рассматривались три варианта:

  1. При моделировании эксперимента предполагалось, что температура графита и концентрация ксенона в топливе постоянны и равны начальным величинам. Перед измерением мощность реактора не меняется, температура графита постоянная, концентрация ксенона равновесная. (Вариант 1).
  2. При моделировании эксперимента учитывалось изменение концентрации ксенона и температуры графита. Перед измерением мощность реактора не меняется, температура графита постоянная, концентрация ксенона равновесная. (Вариант 2).
  3. При моделировании эксперимента учитывалось изменение концентрации ксенона и температуры графита. За 200 секунд перед измерением в результате перемещения 4 стержней АР с начальной глубины погружения 350 см до 354 см была внесена отрицательная реактивность, которая вызвала переходный процесс со снижением мощности реактора (около 2.0 МВт/мин) и изменением температуры графита и концентрации ксенона. (Вариант 3).

Результаты расчетов для трех вариантов представлены в таблице 5.

Таблица 5.
Результаты моделирования эксперимента по измерению..

№ опыта Вариант 1 Вариант 2 Вариант 3
1 0.38 0.42 0.42
2 0.38 0.35 0.33
3 0.36 0.41 0.45
4 0.38 0.30 0.26

Видно, что во 2 и 3 вариантах наблюдается систематическое занижение величин в четных опытах по отношению к нечетным. Это связано с тем, что на двухминутном интервале уже начинает проявляться обратная связь по температуре графита и концентрации ксенона. Такая же картина наблюдается в реальных экспериментах. Видно также, что фоновое снижение мощности (вариант 2) приводит к некоторому увеличению «раскачки». Результаты расчета варианта 1 с постоянными величинами температуры графита и концентрации ксенона показали, что при «отключенной» обратной связи по температуре графита и концентрации ксенона величины, полученные в четных и нечетных опытах практически совпадают.

Исследуется влияние изменения давления в барабан-сепараторе на результаты измерения w. В стационарных расчетах w давление в барабан-сепараторе предполагается постоянным. Согласно экспериментальным результатам при измерении w давление в барабан сепараторе снижается на 0.2 – 0.4 атмосферы. Сравниваются результаты расчетов, моделирующих измерение w в предположении, что Рбс не меняется в процессе эксперимента и расчетов, учитывающих изменение Рбс в процессе эксперимента. На рисунках 5 и 6 приведены результаты моделирования экспериментов по измерению w на 1 блоке САЭС 06.03.2006. Согласно результатам эксперимента при моделировании предполагалось, что давление в барабан сепараторе линейно снижается на 0.2 атмосферы (с 70 до 69.8) в течение 40 секунд после начала погружения стержней.

Рис. 5. Изменение реактивности при моделировании измерения w.

Рис. 6. Изменение мощности при моделировании измереняи w.

Видно, что учет изменения давления в барабан-сепараторе приводит к заметному изменению временных зависимостей реактивности и мощности реактора. Более быстрый рост реактивности в случае учета изменения давления приводит к более выраженному локальному максиму мощности в момент времени 3540 секунд после начала движения стержней. Учет изменения давления в процессе измерения w приводит к тому, что величина w оказывается больше по абсолютной величине на 20 30% по сравнению с расчетом варианта с постоянным давлением в барабан сепараторе. Таким образом, при моделировании эксперимента по определению значения w необходимо учитывать изменение давления в барабан-сепараторе.

Численное моделирование экспериментов по измерению с

Далее моделируются эксперименты по измерению графитового температурного коэффициента реактивности с. проведенные на блоках Ленинградской, Смоленской и Курской АЭС. Моделирование эксперимента проводилось по специальной версии программы STEPAN, предназначенной для расчетов медленных переходных процессов (выгорание и перегрузка топлива, ксеноновые переходные процессы). Как и в эксперименте, в расчете моделировалось изменение температуры графита. Внесенная при этом реактивность компенсировалась перемещением регулирующих стержней. Результаты расчетов представлены в таблице 7.

Таблица 7.

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

Параметр, размерность Станция, номер блока и дата
ЛАЭС-1, 18.12.03 КАЭС-1, 30.03.01 КАЭС-2, 30.06.04 САЭС-2, 01.11.01 САЭС-2, 11.01.03 САЭС-3, 19.09.01
c, 1/0С (мод) 2.610-5 2.910-5 3.510-5 5.710-5 3.410-5 4.210-5
c, 1/0С (ст.) 4.510-5 4.610-5 3.310-5 4.310-5 4.410-5 4.310-5
c, 1/0С (экс.) 2.710-5 3.7410-5 4.710-5 4.810-5 3.3410-5 4.410-5

В таблице 7 приведены значения с, полученные при моделировании эксперимента с (мод.), полученные из стационарного расчета по стационарной версии программы STEPAN с (ст.) и полученные из эксперимента с (экс.). Величины с, полученные экспериментальным путем и в результате численного моделирования эксперимента, имеют довольно большой разброс, что, по видимому, объясняется трудностью точного измерения параметров, необходимых для определения с. Величины с, полученные в результате стационарного расчета имеют гораздо меньший разброс и лежат в пределах изменения экспериментальных значений с. Это говорит об отсутствии систематической ошибки в стационарном расчете.

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

Комбинированный метод определения с и w.

Представлено предложение по методу измерения быстрого мощностного и графитового коэффициентов реактивности, основанному, как и штатный метод измерения w, на вводе отрицательной реактивности в результате частичного погружения четверки регулирующих стержней. В штатном методе определения быстрого мощностного коэффициента тепловая мощность реактора после достижения локального максимума (в период времени 20-25 секунд после начала погружения стержней) начинает снова монотонно снижаться. Это обусловлено тем, что стабилизация мощности на новом уровне при отключенных медленных обратных связях наступает лишь через 100-200 секунд (рисунок 4). Можно записать следующую систему уравнений:

i=wWi+cTi+XXei ; i=1,3

В левой части (i) стоит реактивность, внесенная в результате изменения мощности, температуры графита и концентрации ксенона (1 2 и 3 члены в правой части уравнения). Три уравнения составляются для трех разных моментов времени. К сожалению, решить эту систему уравнений достаточно сложно. Дело в том, что вклады в реактивность в результате изменения температуры графита и концентрации ксенона (величины cTc и XXe соответственно) имеют похожие временные зависимости.

Можно переписать систему уравнений следующим образом:

i=wWi+c*Ti ; i=1,2

где c* - эффективный коэффициент реактивности по температуре графита, учитывающий также вклад в реактивность от изменения концентрации ксенона. Это можно сделать потому, что если изменение мощности по форме близко к ступеньке, то в течение некоторого времени после изменения мощности временные зависимости изменения температуры графита и концентрации ксенона близки к линейным. По сути, коэффициент реактивности c* можно рассматривать как характеристику медленной положительной обратной связи по мощности реактора. Для проверки метода было проведено численное моделирование эксперимента по измерению w на 3 блоке ЛАЭС 15 ноября 2005 года. В результате были получены следующие значения:

w=-2.410-4 -2.610-4 /МВт и c*=8.910-3 9.910-3 /0С. Полученная величина быстрого мощностного коэффициента реактивности хорошо согласуется с экспериментальной величиной w=-2.310-4 /МВт. Согласно стационарному расчету нейтронно-физических характеристик реактора для данного состояния, эффект разотравления реактора на мощности (-X) равен 5.5. Подставляя эту величину в уравнение с*Tc= сTc-5.5Xe, можно получить величину графитового температурного коэффициента реактивности с. Таким образом, получаются значения с =4.110-3 4.310-3 /0C. Согласно стационарному расчету, проведенному по стационарной версии программы STEPAN, величина с=3.710-3/0C.

Пассивный метод определения с и эффекта ксенонового отравления.

Рассмотрим теперь возможный альтернативный метод определения c без проведения специального эксперимента. Он может быть использован на любом энергетическом уровне мощности реактора, в том числе, на номинальном. Определение c проводится из анализа процесса перевода реактора с одного уровня мощности на другой. При этом, как и в традиционном методе, измеряется изменение ОЗР и определяется величина реактивности, внесенной в результате изменения уровня мощности, температуры графита и концентрации ксенона. Можно записать следующее уравнение:

Где - суммарное изменение реактивности;

wW – изменение реактивности за счет изменения мощности;

сTc – изменение реактивности за счет изменения температуры графита;

XXe - изменение реактивности за счет изменения концентрации ксенона.

Суммарное изменение реактивности () может быть выражено через изменение ОЗР:

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

К достоинствам этого метода можно отнести:

- этот метод пассивный (не требует проведения специальных экспериментов);

- существенно возрастают возможности регулярного измерения с;

- этот метод позволяет изучить зависимость с от уровня мощности (как показывают расчеты, с имеет существенную зависимость от уровня мощности реактора);

- одновременно экспериментально определяется величина отравления Xe135, ранее экспериментально не определявшаяся.

Усовершенствованная методика определения .

Далее обсуждается новая методика обработки экспериментальных данных при измерении парового коэффициента реактивности. При определении меняется температура теплоносителя на входе в реактор. Считается, что внесенная при этом реактивность полностью компенсируется в течение 2 минутного интервала в результате перемещения стержней и изменения мощности реактора и мощность реактора стабилизируется на новом уровне. Однако, как показывает практика, полной компенсации за этот интервал времени не происходит. Это может приводить к значительной «раскачке» результатов при задании возмущения разного знака. Предлагается вносимую реактивность измерять непосредственно реактиметром, а не оценивать по перемещению стержней и изменению мощности реактора. При этом время одного измерения снижается с 2 минут до 30 секунд. Показывается, что таким образом можно полностью избавиться от зависимости измеряемой величины

 
Похожие работы:

«УДК 662.997:537.22. НУРМАТОВ ШАВКАТ РАСУЛМАТОВИЧ РАЗРАБОТКА однозеркальных солнечных высокотемпературных технологических установок и технологии СИНТЕЗА КАРБИДОВ 05.14.08 - Энергоустановки на основе возобновляемых видов энергии АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата...»

«ВАСИЛЬЕВ Владимир Владимирович разработка автоматики комплексного аварийного управления нагрузкой Специальность 05.14.02 – Электрические станции и электроэнергетические системы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2011 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет Научный руководитель:...»

«Латыпов Руслан Назымович ТЕПЛОПРОВОДНОСТЬ БЕРИЛЛИЯ ПОСЛЕ ОБЛУЧЕНИЯ ДО ВЫСОКОЙ ПОВРЕЖДАЮЩЕЙ ДОЗЫ Специальность: 05.14.03 - Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации. Автореферат диссертации на соискание ученой степени кандидата технических наук Москва, 2012 Работа выполнена в отделении реакторного материаловедения ОАО Государственный научный центр Научно-исследовательский институт атомных реакторов (ОАО ГНЦ НИИАР) Научный...»

«БОРУШ Олеся Владимировна Эффективность применения парогазовых установок в условиях топливно-энергетического баланса реги о на Специальность 05.14.14 – тепловые электрические станции, их энергетические системы и агрегаты АВТОРЕФЕРАТ Диссертации на соискание учёной степени кандидата технических наук Новосибирск – 2008 Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет Научный...»

«АЮЕВ БОРИС ИЛЬИЧ МЕТОДЫ И МОДЕЛИ ЭФФЕКТИВНОГО УПРАВЛЕНИЯ РЕЖИМАМИ ЕДИНОЙ ЭЛЕКТРОЭНЕРГЕТИЧЕСКОЙ СИСТЕМЫ РОССИИ Специальность 05.14.02 – Электростанции и электроэнергетические системы Автореферат диссертации на соискание ученой степени доктора технических наук Новосибирск – 2008 Работа выполнена в ГОУ ВПО Уральский государственный технический университет – УПИ имени первого Президента России Б.Н. Ельцина Научный консультант: доктор технических наук, профессор Бартоломей Петр...»

«РЫЖКИНА Александра Юрьевна АНАЛИЗ И СОВЕРШЕНСТВОВАНИЕ ХРОМАТОГРАФИЧЕСКИХ МЕТОДОВ ДИАГНОСТИКИ ВЫСОКОВОЛЬТНОГО МАСЛОНАПОЛНЕННОГО ЭЛЕКТРООБОРУДОВАНИЯ Специальность 05.14.12 – Техника высоких напряжений АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2012 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет Научный...»

«Луконин Алексей Валерьевич Р елейная защита закрытых электроустановок напряжением 0,4-10 кВ с распознаванием повреждений, сопровождаемых электрической дугой Специальность 05.14.02 – Электростанции и электроэнергетические системы Автореферат диссертации на соискание ученой степени кандидата технических наук Новочеркасск 2009 г. Работа выполнена в ГОУ ВПО Южно-Российский государственный технический университет (Новочеркасский политехнический институт) на кафедре Электрические...»

«УДК 662.997:537.22. ТЎРАЕВА ЎЛМАСОЙ ФАРМОНОВНА ЭФФЕКТИВНОСТЬ СОЛНЕЧНЫХ теплоэнергетических УСТАНОВОК С СЕЛЕКТИВНЫМИ ПРИЕМНИКАМИ ИЗЛУЧЕНИЯ 05.14.08 - Энергоустановки на основе возобновляемых видов энергии 05.14.05 - Теоретические основы теплотехники АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук...»

«УДК 574:539.1.04+621.039.7 Семенов Сергей Геннадьевич РАЗРАБОТКА И ПРИМЕНЕНИЕ МЕТОДОВ ВЫВОДА ИЗ ЭКСПЛУАТАЦИИ И РЕАБИЛИТАЦИИ ЯДЕРНО- И РАДИАЦИОННО-ОПАСНЫХ ОБЪЕКТОВ НИЦ КУРЧАТОВСКИЙ ИНСТИТУТ С ИСПОЛЬЗОВАНИЕМ ДИСТАНЦИОННО УПРАВЛЯЕМЫХ РОБОТОТЕХНИЧЕСКИХ МЕХАНИЗМОВ Специальность: 05.14.03 – Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации АВТОРЕФЕРАТ диссертации на соискание научной степени кандидата технических наук Москва - 2013...»

«СТЕПЕННОВ Дмитрий Борисович СОЗДАНИЕ И ПРИМЕНЕНИЕ центра технической поддержки объектов использования атомной энергии в ниц курчатовский институт Специальность 05.14.03 – Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации Автореферат диссертации на соискание ученой степени кандидата технических наук Москва – 2012 Работа выполнена в Национальном исследовательском центре Курчатовский институт (НИЦ Курчатовский институт) (г. Москва)....»

«Алексеев Александр Вениаминович ИСПЫТАНИЯ В РЕАКТОРЕ МИР ТВЭЛОВ ВВЭР-1000 В РЕЖИМЕ АВАРИИ С ВВОДОМ ПОЛОЖИТЕЛЬНОЙ РЕАКТИВНОСТИ Специальность 05.14.03 – ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации Автореферат диссертации на соискание ученой степени кандидата технических наук Димитровград – 2011г. Работа выполнена в открытом акционерном обществе Государственный научный центр - Научно-исследовательский институт атомных реакторов,...»

«Кузин Юрий Сергеевич ПОВЫШЕНИЕ ЭКСПЛУАТАЦИОННОЙ НАДЕЖНОСТИ РЕГУЛИРУЮЩИХ КЛАПАНОВ ДИСКОВОГО ТИПА ДЛЯ ТРУБОПРОВОДОВ ТЭС И АЭС 05.14.14 Тепловые электрические станции, их энергетические системы и агрегаты Автореферат диссертации на соискание ученой степени кандидата технических наук Новочеркасск-2012 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Южно-Российский государственный технический университет...»

«БЕЛОГЛАЗОВ Алексей Владимирович Разработка адаптивных средств выявления неисправностей и стратегии обслуживания гидроагрегатов Специальность 05.14.02 – Электрические станции и электроэнергетические системы АВТОРЕФЕРАТ диссертация на соискание ученой степени кандидата технических наук Новосибирск – 2011 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет...»

«Айзатулин Амир Исмаилович СРЕДСТВА МОДЕЛИРОВАНИЯ И ПРОЕКТИРОВАНИЯ АЛГОРИТМОВ АСУ ТП ЭНЕРГОБЛОКА АЭС И СИСТЕМА ВИЗУАЛИЗАЦИИ И УПРАВЛЕНИЯ ДЛЯ МОДЕЛИРУЮЩИХ ПРОГРАММНЫХ КОМПЛЕКСОВ Специальность 05.14.03. – ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва 2006 Работа выполнена в ОАО Всероссийский научно-исследовательский институт по эксплуатации...»

«АГЕЕВ Михаил Александрович ОПТИМИЗАЦИЯ СХЕМ И РАБОЧИХ ПАРАМЕТРОВ СИСТЕМ ГЛУБОКОЙ ОЧИСТКИ ДЫМОВЫХ ГАЗОВ ПРИ МОДЕРНИЗАЦИИ ТВЕРДОТОПЛИВНЫХ ПРОМЫШЛЕННЫХ ЭНЕРГОУСТАНОВОК Специальность 05.14.04 – Промышленная теплоэнергетика Автореферат диссертации на соискание ученой степени кандидата технических наук Саратов – 2010 Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Саратовский государственный технический университет Научный...»

«СКЛАДЧИКОВ Александр Александрович ОЦЕНКА НАДЕЖНОСТИ И УПРАВЛЕНИЕ РИСКАМИ ТЕХНОЛОГИЧЕСКИХ НАРУШЕНИЙ НА ВОЗДУШНЫХ ЛИНИЯХ ЭЛЕКТРОПЕРЕДАЧИ Специальность 05.14.02 – Электрические станции и электроэнергетические системы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Чебоксары – 2012 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Самарском государственном техническом...»

«Тутундаева Дарья Викторовна МОНИТОРИНГ ДОПУСТИМОСТИ ПОСЛЕАВАРИЙНЫХ РЕЖИМОВ ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ СИСТЕМ Специальность 05.14.02 – Электрические станции и электроэнергетические системы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2011 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет

«Михеев Павел Александрович ПРОДОЛЬНОЕ И ПОПЕРЕЧНОЕ ТОКООГРАНИЧЕНИЕ В ЭЛЕКТРИЧЕСКИХ СИСТЕМАХ С ПОМОЩЬЮ СВЕРХПРОВОДНИКОВЫХ УСТРОЙСТВ Специальность 05.14.02 – Электростанции и электроэнергетические системы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск - 2008 Работа...»

«Малков Андрей Павлович Обеспечение ядерной безопасности водоохлаждаемых исследовательских реакторов Специальность: 05.14.03. Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации Автореферат диссертации на соискание ученой степени доктора технических наук Димитровград 2012 Работа выполнена в Государственном научном центре Научно-исследовательском институте атомных реакторов Научный консультант: доктор технических наук, профессор...»

«МЯТЕЖ аЛЕКСАНДР ВЛАДИМИРОВИЧ РЕГУЛИРОВАНИЕ НАПРЯЖЕНИЯ В СИСТЕМАХ ЭЛЕКТРОСНАБЖЕНИЯ С ИСПОЛЬЗОВАНИЕМ НЕЧЕТКОЙ ЛОГИКИ Специальность 05.14.02 – Электростанции и электроэнергетические системы АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2009 Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Новосибирский государственный технический университет Научный руководитель: доктор технических...»








 
2014 www.avtoreferat.seluk.ru - «Бесплатная электронная библиотека - Авторефераты диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.