Авторефераты по всем темам >>
Авторефераты по разным специальностям
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Согласно данным Международного агентства по атомной энергии (МАГАТЭ), по профилю выработанного календарного ресурса энергоблоков атомных электростанций (АЭС) в настоящее время в эксплуатации находятся 148 блоков, проработавших более 30 лет, и 210 блоков, возраст которых составляет от 20 до 30 лет.
Суммарно эти энергоблоки представляют более 80% от 440 реакторов, эксплуатирующихся в мире. С целью безопасного и надёжного функционирования энергоблоков России в отрасли разработаны программы исследований по переоценке безопасности, управлению старением, оптимизации технического обслуживания и ремонта, оценке остаточного ресурса и продлению сроков службы АЭС. В рамках данных направлений возрастает ценность работ по анализу надёжности объектов ядерной энергетики и, как следствие, возникает необходимость разработки методов анализа ресурсных характеристик оборудования АЭС с целью определения их технического состояния.
Разработанные в представленной диссертации неасимптотические математические методы могут быть использованы для оценки ресурсных характеристик элементов и систем ядерных энергетических установок (ЯЭУ). Исследования ресурсных характеристик именно в неасимптотической постановке позволяют более достоверно описать процесс функционирования оборудования энергоблоков АЭС, учесть влияние различных режимов эксплуатации на работоспособность систем, оптимизировать периоды между плановыми профилактическими работами (ППР).
Кроме этого, описанные в диссертации неасимптотические модели являются более гибкими, чем большинство ранее известных математических моделей, и ориентированы на проведение исследований ресурсных характеристик и других показателей надежности уникального высоконадежного оборудования в условиях неполноты и ограниченности исходных данных об отказах.
Таким образом, актуальность работы заключается в том, что на сегодняшний день существует большое количество оборудования АЭС, для которого требуется проводить расчеты ресурсных характеристик с высокой достоверностью и точностью. Применение неасимптотических моделей является необходимым условием повышения точности расчетов показателей надежности элементов и систем ЯЭУ для решения задач по оценке технического состояния действующих энергоблоков при наличии ограниченной, неполной исходной информации, а также позволяет учесть уникальность и особенности эксплуатации оборудования, специфичные для АЭС.
Объектом исследования представленной работы являются элементы и подсистемы энергоблоков АЭС, рассматриваемые как сложные технические системы с различными стратегиями функционирования, в ходе эксплуатации которых имеют место контроль исправности составляющих элементов, проведение профилактических работ, наличие запасных элементов.
3 Предметом исследования являются неасимптотические модели расчета ресурсных характеристик, учитывающие различные стратегии функционирования элементов и систем ЯЭУ, и методы анализа надежности и технического состояния оборудования АЭС.
Цель исследования состоит в разработке неасимптотических моделей расчета ресурсных характеристик элементов и систем ЯЭУ, используемых для анализа технического состояния и надежности оборудования АЭС.
Для достижения данной цели в диссертации решены следующие задачи:
1. Сформулирована проблема необходимости разработки неасимптотических методов расчета ресурсных характеристик элементов и систем ЯЭУ, предоставляющих возможность оценки технического состояния и научного обоснования продления их проектного срока службы.
2. Разработан метод оценки средней остаточной наработки до отказа невосстанавливаемых объектов, учитывающий особенности сбора статистических данных об отказах оборудования на АЭС.
3. Разработаны неасимптотические модели для оценивания математического ожидания прямого остаточного времени - ресурсной характеристики восстанавливаемых систем, учитывающие особенности эксплуатации ЯЭУ как сложной технической системы (наличие встроенного контроля, планово-профилактического обслуживания).
4. Разработано программное обеспечение для расчета ресурсных характеристик с учетом различных видов представления исходных данных об отказах элементов ЯЭУ. Проведена его апробация на тестовых и реальных данных.
5. Выполнены расчеты ресурсных характеристик для электрооборудования системы управления и защиты (СУЗ) реакторных установок РБМК-1000 Смоленской АЭС (САЭС) и ЭГП-6 Билибинской АЭС (БиАЭС) с использованием разработанных моделей.
Научная новизна представленной работы состоит в том, что впервые были разработаны:
1. Неасимптотические методы оценки ресурсных характеристик, в отличие от ранее существующих подходов, позволяющие получать оценки показателей долговечности для невосстанавливаемых элементов ЯЭУ с учетом особенностей представления исходной статистической информации об отказах оборудования российских АЭС.
2. Неасимптотические модели расчета среднего прямого остаточного времени для восстанавливаемых элементов, учитывающие наличие встроенных систем контроля работоспособности оборудования, проведение обслуживающих мероприятий, и позволяющие, в отличие от известных на сегодняшний день асимптотических моделей, оценить техническое состояние на всех этапах эксплуатации объекта исследования и использовать в качестве исходных данных как параметрические, так и непараметрические модели надежности.
3. Метод построения доверительных интервалов для неасимптотических оценок математического ожидания прямого остаточного времени для случая представления исходных данных об отказах оборудования, в виде частот отказов группированных по временным интервалам.
Практическая значимость работы состоит в том, что:
1. Разработанные теоретические модели доведены до инженерных методик с соответствующей программной реализацией.
2. Разработанные методы использовались при диагностике и анализе технического состояния и надежности элементов и подсистем ЯЭУ реакторов ЭГП-БиАЭС и РБМК-1000 САЭС с учетом особенностей их функционирования и разнородной статистической информации об отказах, а также при обосновании решения о продлении установленных сроков службы электрооборудования СУЗ БиАЭС до 31.12.2016 (Приказ №283 филиала ОАО Концерн Росэнергоатом Билибинская атомная станция от 24.12.2011).
3. Научные результаты включены в методику и типовую программу оценки технического состояния и остаточного ресурса комплекса электрооборудования системы управления и защиты ядерных установок типа ВВЭР (РД ЭО 0294Ц01, изменение №1), в программу-методику работ по обследованию технического состояния и управлению ресурсными характеристиками электрооборудования СУЗ БиАЭС.
Достоверность научных положений обеспечена применением широко известных методов теории вероятности, теории надежности и математической статистики; использованием известных процедур преобразования функционалов при выводе уравнений восстановления и численных методов при их решении; результатами сравнения частных случаев полученных моделей с разработками других авторов; опытом внедрения и практическим использованием полученных результатов.
ичное участие автора Автор участвовал в качестве исполнителя на всех этапах разработки неасимптотических методов расчета и анализа ресурсных характеристик оборудования ЯЭУ. Лично автором получены неасимптотические модели для оценки математического ожидания прямого остаточного времени и других ресурсных характеристик для различных стратегий обслуживания с учетом особенностей представления исходных данных, разработан метод построения доверительного интервала для оценки математического ожидания прямого остаточного времени, показана применимость полученных моделей для различных типов оборудования, а также для решения оптимизационных задач надежности, например, определения оптимального периода между проведением ППР.
Положения, выносимые на защиту:
1. Неасимптотические методы оценки ресурсных характеристик невосстанавливаемых элементов ЯЭУ, учитывающие особенности представления исходных данных об отказах оборудования на российских АЭС.
2. Неасимптотические методы расчета математического ожидания прямого остаточного времени для восстанавливаемых элементов и систем ЯЭУ с учетом различных стратегий их функционирования и особенностей представления исходных данных об отказах оборудования на российских АЭС.
3. Метод построения доверительных интервалов для оценок математического ожидания прямого остаточного времени, основанный на применении процедур бутстреп-метода.
4. Расчетный анализ ресурсных характеристик и других показателей надежности элементов и подсистем СУЗ реакторов ЭГП-6 БиАЭС и РБМК-1000 САЭС с учетом ограниченной исходной информации и различных стратегий функционирования.
Апробация результатов работы Основные результаты диссертации опубликованы в российских научных журналах, докладывались и обсуждались на международных конференциях и получили одобрение российских и зарубежных специалистов по теории надежности.
Среди конференций, на которых представлялись результаты работы:
1. IX Международная научно-практическая конференция Системный анализ в проектировании и управлении, Санкт-Петербург, 2005;
2. Международная конференция Безопасность АЭС и подготовка кадров, Обнинск, 2005, 2007, 2009;
3. III международная конференция Математические идеи П.Л. Чебышева и их приложение к современным проблемам естествознания, Обнинск, 2006;
4. International Conference Degradation, Damage, Fatigue, and Accelerated Life Models in Reliability Testing (ALTТ2006), Angers, France, 2006;
5. V International Conference - Mathematical Methods in Reliability (MMRТ2007), Glasgow, UK, 2007;
6. VI International Conference - Mathematical Methods in Reliability. Theory.
Methods. Applications (MMRТ2009), Mosсow, 2009;
7. International Conference Accelerated Life Testing, Reliability-based Analysis and Design (ALTТ2010), Clermont-Ferrand, France, 2010;
8. International Workshop Applied Methods of Statistical Analysis. Simulation and Statistical Inference (AMSAТ2011), Novosibirsk, Russia, 2011.
Публикации Основные результаты диссертации опубликованы в 16 печатных работах, в том числе в 3 печатных работах в рецензируемых журналах из списка ВАК.
Структура и объем диссертации Диссертация состоит из введения, четырех глав и заключения со сквозной нумерацией внутри каждой главы. Общий объем диссертации составляет 211 страниц машинописного текста. Основной текст работы изложен на 144 страницах, включая 40 рисунков и 15 таблиц. Библиографический список литературы содержит наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введении отражена актуальность разработки неасимптотических методов оценки ресурсных характеристик элементов и систем АЭС на основании неполной, ограниченной информации об отказах, обозначены цели и задачи исследования.
Сформулированы научная новизна и практическая значимость работы, представлены положения, выносимые автором на защиту.
Первая глава посвящена анализу существующих литературных источников и анализу текущего состояния исследований по теме диссертации.
В первом параграфе проводится обзор литературы по проблеме оценивания ресурсных характеристик различного рода оборудования. Рассмотрены работы по параметрическим и непараметрическим методам оценки технического состояния и остаточного ресурса систем и составляющих компонентов. Показано, что решение этой задачи в большинстве случаев сводится к анализу определяющего параметра, в качестве которого выбирается один из показателей безопасности или комплексных показателей надежности. Кроме того, данные методы часто невозможно применить на практике из-за отсутствия необходимого количества исходных данных.
Во втором параграфе проведен обзор наиболее распространенных методов оценки технического состояния и остаточного ресурса оборудования ЯЭУ. В большинстве случаев для анализа используются методы оценивания вероятностных характеристик и проверяются различного рода гипотезы, связанные с предположением о закономерностях в распределении случайных наработок.
Согласно ГОСТ 27.002Ц89, основными показателями, отражающими изменение надежности элементов, являются вероятность безотказной работы (ВБР) P(t) (рис. 1 (а)), интенсивность отказов (t) (рис. 1 (б)) для невосстанавливаемых элементов и параметр потока отказов (t) (ППО) для восстанавливаемых элементов.
а) б) Рис. 1. Анализ вероятностных характеристик надежности с целью оценки остаточного ресурса В качестве одного из показателей надежности конкретного вида оборудования производителем часто указывается предельно допустимое значение ВБР Pпр. Данная характеристика может применяться при обосновании значения ресурса из условий риска (см. рис. 1 (а)). Момент исчерпания ресурса объектом определяется временем Tпр, при котором оценка P(t) станет меньше Pпр.
Также проектантом и эксплуатирующей организацией могут устанавливаются два уровня интенсивности отказов элементов ЯЭУ: контрольный КР, соответствующий наработке TНСт, выработав которую, объект начинает интенсивно стареть, и предельный уровень ПР, определяющий предельную наработку TНР, по достижению которой эксплуатация объекта нецелесообразна.
Выбор ресурсного определяющего параметра зависит от структуры и режима эксплуатации оборудования. При альтернирующем процессе для оценки технического состояния часто используются ведущая функция потока (ВФП) отказов, которая определяется как математическое ожидание числа отказов на заданном интервале времени [0; t] и нестационарный коэффициент готовности.
Методам оценки вероятностных характеристик надежности в условиях неопределенности исходных данных посвящен третий параграф первой главы. В большинстве случаев в процессе сбора эксплуатационной информации исследователь не имеет данных о значениях наработок отказавших объектов, доступны лишь сведения о количестве однотипных элементов одновременно функционирующих в составе системы и количестве их отказов, произошедших на некотором интервале времени. Говоря терминами теории надежности, статистические данные об отказах представлены в виде количества отказов i, совокупности m однотипных элементов, распределенных по s интервалам работоспособности [li; ri], где i=1..s. При этом неизвестно, какой именно объект из группы однотипных элементов оборудования отказал, когда отказавшее устройство было введено в эксплуатацию и т.д.
В параграфе описана процедура непараметрического оценивания плотности и функции распределения наработки до отказа по цензурированным данным об отказах элементов с помощью ядерной оценки ППО s i (t) = m(r li ) t - li G t - ri + (t), (1) G h - h i i=где h - параметр локальности (мера, зависящая от среднего квадратичного отклонения случайной величины наработки до отказа), x 1 u2 du G(x) = exp-, 2 - (t) - систематическая ошибка оценки ППО, которая определяется выражением 1 n m - t 2t n m + t , (t) = erfc + exp erfc 2 2 2n m 2n m где - математическое ожидание наработки до отказа, - настроечный параметр.
Как известно из теории восстановления, ППО связан с плотностью распределения интегральным уравнением Вольтерра второго рода t f (t) = (t) - f (u)(t - u)du. (2) Данный метод оценки плотности распределения наработки до отказа используется в представленной работе при расчетах ресурсных характеристик.
Таким образом, в первой главе сделан вывод о том, что при решении практически всех задач надежности используются классические теоретические модели оценки определяющего параметра, в качестве которого обычно используются показатели безопасности или комплексные показатели надежности. При этом в основе применяемых математических моделей лежит гипотеза о виде закона распределения наработки до отказа, что является существенным ограничением. Неполнота и ограниченность исходных данных, которые встречаются на практике в большинстве случаев, часто не позволяют получить необходимые результаты.
Вторая глава посвящена разработке методов оценки ресурсных характеристик для невосстанавливаемых элементов и систем с мгновенным восстановлением.
В первом параграфе выполнен обзор математических методов расчета ресурсных характеристик невосстанавливаемых элементов наиболее часто применяемых на практике в последние годы. Данные методы фактически представляют собой подходы классической теории надежности, адаптированные к условиям задачи практической оценки технического состояния реальных объектов.
Во втором параграфе представлена методика оценки остаточной наработки до отказа невосстанавливаемых элементов при однородном потоке отказов и ограниченной исходной информации об отказах.
Как известно, остаточная наработка до отказа является показателем долговечности невосстанавливаемых элементов и может быть вычислена как T (t) = P()d. (3) P(t) t Однако применить данное выражение на практике в случаях, когда исходные данные об отказах представлены цензурированными выборками, достаточно сложная задача. Даже если использовать ядерную оценку (1), то возможно получить оценку ВБР только на интервале наблюдений [0, t], в то время как для вычислений по формуле (3) необходимо определить ВБР на всей области определения [0, ).
В данном параграфе предложен подход, позволяющий для расчета остаточной наработки до отказа перейти от оценки (3) к оценке (4) t T (t) = Tср - ()d, (4) P P (t) где P(t) - ВБР элемента, а Tср - средняя наработка на отказ.
Оценки (1)Ц(2), (4) дают возможность проводить расчет функции средней остаточной наработки до отказа в случае, когда исследователю доступны лишь данные о количестве отказов совокупности однотипных элементов, произошедших на фиксированном интервале времени, что характерно для всех российских АЭС.
В третьем параграфе приведен метод оценки математического ожидания прямого остаточного времени MVt(t) - средней остаточной наработки работающего элемента, достигшего возраста t, для случая, когда восстановление происходит мгновенно. Данная стратегия характерна для элементов, которые в случае отказа заменяются на аналогичные из состава запасных изделий и приборов (ЗИП). При этом время замены элемента мало по сравнению с наработкой до отказа.
Согласно определению математического ожидания случайной величины, Vt (t) = i i i+1 i (X - t) I{X t < X } = (t), i=0 i=где i (t) = (Xi - t) I{Xi t < Xi+1}, Xi - момент i-го отказа элемента.
Тогда математическое ожидание i(t) равно t i (t) = + - t){u t < u + } fXi (u) f()ddu = f(x + t - u) fXi (u)dudx = (t; x)dx.
i (u x x 0 0 0 0 Преобразование Лапласа внутреннего интеграла i(t; x) равно t - pt i ( p; x) = f (x + t - u) f (u)dudt = fi ( p) g( p; x), Xi e 0 где g( p; x) - образ функции g(t; x)=f(t+x).
Применив преобразование Лапласа к i(t), можно получить i ( p) = fi ( p) p; x)dx.
xg( Тогда изображение по Лапласу функции MVt(t)будет определяться выражением Vt ( p) = fi ( p) p; x)dx = xg( xg( p; x)dx.
1- f ( p) i=0 Перейдя к оригиналам, можно получить интегральное уравнение t Vt (t) = (t + x)dx + f (t - х) Vt (х)du. (5) xf 0 Решение данного уравнения позволяет определить математическое ожидание прямого остаточного времени - среднего времени работы, оставшегося до очередного отказа, начиная с момента времени t, в который объект был работоспособен.
В четвертом параграфе решается проблема определения точности оценок MVt(t) восстанавливаемого оборудования для случая, когда исходная информация представлена в цензурированном виде. Автором предлагается метод построения доверительного интервала с заданной доверительной вероятностью для непараметрической оценки (5), основанный на бутстреп-процедуре. Основная идея разработанного метода состоит в отображении равномерно распределенной случайной величины, смоделированной на оси вероятности, на ось времени, разбитую на интервалы, с известным количеством попаданий наблюдаемой случайной величины.
Рассмотрим последовательность действий, выполняемых при определении точ- ности оценок с помощью бутстреп-метода.
Этап 1. Подготовка исходных данных для бутстреп-метода. Пусть имеется n от- казов, распределенных по s временным интервалам. Перенумеруем события, определяющие отказы объекта, в порядке возрастания времени их реализации, то есть каждому событию поставим в соответствие индекс i. Определим временной интервал, на который попало данное событие. Если на один интервал времени приходится несколько отказов, порядок нумерации событий внутри интервала произвольный.
Этап 2. Ось вероятностей делим на n равных непересекающихся интервалов [0, y1), [y1, y2),..., [yn-1, yn]. Моделируем случайную равномерно распределенную величину Uc[0,1] на оси вероятностей. Далее определяем, какой из зафиксированных отказов соответствует смоделированной переменной. Если смоделированная случайная величина принадлежит интервалу [yi-1, yi), это означает, что реализовалось событие с индексом i. Повторяем данную процедуру моделирования n раз, получаем выборку номеров событий отказов. Далее наблюдаемое событие должно быть поставлено в соответствие временному интервалу наблюдения согласно описанной на первом этапе процедуре. Используя полученную с помощью бутстреп процедуры выборку, производим вычисление оценки среднего прямого остаточного времени MV t1.
Этап 3. Повторяем многократно этап 2. Получаем множество оценок MVt j, где j=1,...,r, r - количество бутстреп повторений третьего этапа.
Этап 4. Непараметрические оценки границ доверительного интервала определяем следующим образом. Вначале задаем уровень значимости , соответствующий доверительной вероятности 1Ц2. Далее определяем границы интервалов, удовлетворяющие следующим соотношениям для заданных d(Vt i Vt н) = ; (6) r d(Vt i Vt в); (7) 1- = r где r - количество бутстреп-повторений, d(Vt i Vt н ) и d(Vt i Vt в ) определяют номер рассчитанной функции среднего прямого остаточного времени, для которого впервые выполнится условие (6), и номер рассчитанной функции среднего прямого остаточного времени, для которого впервые выполнится условие (7), соответственно. В этом случае, значения Vt н и Vt в, определенные выражениями (6) и (7), будут характеризовать приблизительные границы доверительного интервала, соответствующие доверительной вероятности 1Ц2.
Третья глава посвящена разработке неасимптотических моделей оценки ресурсных характеристик, учитывающих различные стратегии функционирования оборудования. Представлен метод вывода уравнений для математического ожидания прямого остаточного времени, описывающих данные модели.
В первом параграфе проведен анализ методов и моделей, используемых с последние годы на практике для оценки параметров надежности элементов и систем, функционирование которых описывается альтернирующими процессами восстановления. В результате отмечено, что для оценки технического состояния оборудования в большинстве случаев исследование строится на оценивании нестационарного коэффициента готовности (неготовности) при параметрически заданных законах распределения наработки до отказа и времени восстановления.
Во втором параграфе рассматривается модель, описывающая функционирование оборудования со встроенной аппаратурой постоянного контроля работоспособности, которая позволяет мгновенно определить событие отказа.
Пусть объект функционирует до момента отказа i, после чего проводится аварийное восстановление, которое длится промежуток времени i (см. рис. 2). После полного восстановления в момент времени si объект продолжает свою работу до очередного момента отказа. Такой цикл повторяется до выбранного момента времени t. Наработка до отказа на i-м шаге равна i.
1 1 2 л+ Ц л+ Ц л+ t 0 1 s1 t Vt 2 sРис. 2. Стратегия функционирования объекта с мгновенным обнаружением отказа Неасимптотическая оценка среднего прямого остаточного времени для данной стратегии функционирования может быть найдена с помощью математической модели, учитывающей взаимосвязь описанных случайных величин через определение математического ожидания остаточной наработки случайной величины Vt.
Применение преобразований Лапласа позволяет получить следующее интегральное уравнение Вольтерра второго рода, описывающее поведение системы с мгновенным обнаружением отказа:
t t-х Vt (t) = (t + x)dx + (х) f (t - х - ) f()ddх. (8) t xf V 0 0 Здесь f(t) - плотность распределения наработки до отказа и f(t) - плотность распределения времени восстановления, соответственно.
В третьем параграфе исследовались системы, для которых характерно наличие элементов контроля неисправностей без обнаружения места отказа (см. рис. 3). В случае отказа некоторого элемента данное событие незамедлительно идентифицируется, но требуется время для обнаружения отказавшего элемента.
В момент начала работы t0=0 система находится в работоспособном состоянии, и контроль осуществляется в системе постоянно. В случае отказа объекта в момент времени i, факт отказа становится известен (например, срабатывает аварийная индикация), однако какой именно элемент системы отказал неопределено, и оперативному персоналу требуется время Si на поиск места отказа. После обнаружения места отказа происходит восстановление системы в течение случайного времени ri. Наработка до отказа объекта при этом равна i.
1 S1 r1 2 S2 r2 л+ Ц Ц л+ Ц Ц л+ 0 1 s1 t Vt 2 s2 t Рис. 3. Стратегия функционирования объекта, учитывающая встроенный контроль с последующей идентификацией места отказа Воспользовавшись определением математического ожидания случайной величины, и выполнив преобразования как в предыдущем рассмотренном случае можно получить интегральное уравнение t t-x t-x-y Vt (t) = (t + x)dx + (x) f (x) f (y) f (t - x - y - z)dzdydx, (9) t s r xf V 0 0 0 где f(t) - плотность распределения наработки до отказа; fr (t) - плотность распределения времени восстановления; fs (t) - плотность распределения времени обнаружения места отказа. Решая уравнение (9), можно оценить математическое ожидание прямого остаточного времени MVt для данной стратегии.
В следующем параграфе рассматривается стратегия функционирования системы с учетом проведения ППР (см. рис.4).
T (1 < T) T (2 > T) T (3 < T) 1 2 t 0 1 f1 s1 2 r2 s2 t Vt 3 f3 s Простой системы Рис.4. Стратегия эксплуатации технической системы с возможностью проведения ремонтных работ только при остановленном оборудовании Отличие от предыдущих описанных случаев состоит в том, что оперативный персонал может обнаружить место отказа и приступить к ремонту через строго фиксированные промежутки времени длительностью T в период ППР, который длиться случайное время r i,если в период функционирования отказа не произошло. Если в течение периода до начала очередных профилактических работ происходит отказ, то система простаивает до начала ППР, во время которого проводится аварийное восстановление системы в течении случайного времени . После аваf i рийного или профилактического восстановления системы si происходит перепланировка момента профилактики.
В работе показано, что для данной стратегии MVt(t) определяется из решения интегрального уравнения t t t V (t) = (t + x)dx + (x) f(t - x)dx. (10) xf V 0 Здесь f(x) = f f (x -T )F (T) + fr (x -T )(1- F (T)), f (x) = (T - x) f (x), где f(t) и F(t)Ц плотность и функция распределения наработки до отказа; fr (t) и f f (t) - плотности распределения времени профилактического и аварийного восстановления, соответственно; H(x) - функция Хевисайда.
Проиллюстрируем поведение математического ожидания прямого остаточного времени для стратегии, учитывающей плановую профилактику (10), на тестовом примере для различных периодов ППР: T=1 год, T=2 года. Для распределения наработки до отказа, времен аварийного и профилактического восстановления было выбрано усеченное нормальное распределение с параметрами: m=1, =0,3;
m f =0,5, f =0,1; mr =0,1, r =0,03.
а) б) Рис.5. Среднее прямое остаточное время MVt(t) с учетом периода профилактического обслуживания: а) Т=1 год, б) Т=2 года Поведение MVt(t) имеет затухающий колебательный характер и переходит в асимптотический режим с увеличением времени эксплуатации системы (рис. 5).
Кроме зависимости от распределений наработки до отказа и времени восстановления, очевидна зависимость среднего прямого остаточного времени от периода проведения профилактических работ. Если проводить анализ функций MVt(t) для различных значений периода T, то следует обратить внимание на моменты времени, соответствующие локальным минимумам математического ожидания прямого остаточного времени. Интервал времени между локальными минимумами функции MVt(t), которые соответствуют максимальному значению вероятности отказа, больше при T=2 года (см. рис. 5 б)), чем при T=1 год (см. рис. 5 а).
При малых значениях T функция MVt(t) быстро убывает вследствие нахождения системы в ремонте по причине профилактики, при больших значениях система после отказа простаивает до момента начала ППР, что отрицательно сказывается на MVt(t). Однако можно подобрать такое оптимальное значение периода профилактики, что среднее значение интегральной характеристики будет максимально.
Таким образом, в третьей главе был представлен метод построения неасимптотических моделей оценки математического ожидания прямого остаточного времени систем с различными стратегиями функционирования. Описанный метод можно использовать для оценки технического состояния сложных технических систем, к которым можно отнести не только оборудование АЭС, но и сложные комплексы авиатехники, космических аппаратов, программные комплексы, компьютерные системы, сети, радиотехнику и электронику, узлы робототехники.
В четвертой главе приведены результаты анализа ресурсных характеристик и других показателей надежности электрооборудования СУЗ реакторов ЭГП-6 БиАЭС и РБМК-1000 САЭС в условиях ограниченности исходных данных.
В первом параграфе четвёртой главы представлено описание состава и принципа функционирования СУЗ реактора РБМК-1000 первой очереди САЭС, приведены функциональные и структурные схемы основных каналов СУЗ и построенные на их основе логические функции работоспособности (ЛФР).
Во втором параграфе представлен анализ результатов оценки ресурсных характеристик элементов и подсистем СУЗ реактора РБМК-1000. В качестве исходной информации использовались данные об отказах электрооборудования СУЗ энергоблока №1 САЭС, которые представляют собой группированные выборки малого объема, элементами которых являются число отказов элементов.
Для элементов СУЗ реактора РБМК-1000 энергоблока №1 САЭС были выполнены расчеты следующих вероятностных характеристик надежности: ППО (гистограммная и ядерная оценки), плотность и функция распределения наработки до отказа, ВБР, а также остаточный ресурс элементов СУЗ. Расчеты проводились почти по всем элементам СУЗ за исключением тех, для которых отсутствовала информация об их количестве в составе системы. Например, для блоков синхронного перемещения (БСП) на рис. 6 (а) представлен график оценки ВБР P(t) (сплошная линия), а также верхняя и нижняя границы доверительного интервала (пунктирные линии). На рис. 6 (б) представлены расчеты математического ожидания остаточной наработки до отказа T(t) и ее тренд (в годах). По оси абсцисс ведется отсчет времени в годах с 1990 по 2006 года.
Согласно полученным оценкам, ВБР элемента БСП, проработавшего 2 года без отказов (от начала установки этого элемента в систему), приближенно равна 0,8; лет без отказов - приближенно равна 0,72. Статистически данные оценки можно понимать следующим образом: из некоторого (большого) количества (допустим 100) одновременно начавших эксплуатироваться элементов БСП примерно 28 элементов через 10 лет откажут, 72 продолжат функционировать.
ВБР (БСП) Остаточный ресурс (БСП) 0,0,0,0,0,0,0,0,0,0 2 4 6 8 10 12 14 t, год 0 5 10 T(t) Tтр(t)=15-0,8t t, год a) б) Рис. 6. Оценка показателей надежности БСП Средняя остаточная наработка до отказа элемента, с увеличением времени эксплуатации снижается, но не обязательно линейно. Проанализировав линейный тренд для T(t), можно говорить о том, что вновь введенный в эксплуатацию БСП проработает до отказа в среднем 15 лет. В случае если элемент проработал 5 лет, то остаточная наработка будет равна 15Ц0,8*5=11 лет, если он проработал 10 лет, то T(10)=15Ц0,8*10=7 лет и т.д.
Также выполнены расчеты структурной надежности для четырех основных схем аварийной защиты (АЗ): аварийная защита по скорости в пусковом диапазоне (АЗСМ) и основном диапазоне (АЗСР), аварийная защита по превышению мощности в пусковом диапазоне (АЗММ) и в основном диапазоне (АЗМР), для которых были составлены схемы надежности, построены ЛФР и рассчитаны вероятностные характеристики надежности.
Полученные зависимости ВБР (рис. 7 (а)) и среднего прямого остаточного времени (рис. 7 (б)) для АЗСР, АЗСМ, АЗМР и АЗММ. Здесь тонкая сплошная линия соответствует графикам показателей надежности для канала АЗСР, тонкая пунктирная линия - каналу АЗММ, толстая сплошная линия - каналу АЗСМ, а толстая пунктирная линия - каналу АЗМР. При расчетах считалось, что время восстановления работоспособности элемента подчиняется нормальному усеченному (слева) закону с параметрами m f =8 ч, f =2 ч.
а) б) Рис. 7. Оценка показателей надежности подсистем СУЗ реактора РБМК-1000 на примере энергоблока №1 САЭС T(t), год P(t) Результаты расчетов показывают, что оборудование СУЗ реактора РБМК-10энергоблока №1 САЭС обладает достаточным ресурсом для продолжения безопасной эксплуатации системы. Из анализа оценок MVt(t) (см. рис. 7 (б)) видно, что наименьшим ресурсом обладает оборудование системы АЗСР, однако, значение математического ожидания прямого остаточного времени на 17 год эксплуатации составило, примерно, 10 лет.
Третий параграф данной главы посвящен анализу надежности электрооборудо- вания СУЗ реакторов ЭГП-6 БиАЭС, выполненного на основании расчетов ресурсных характеристик и других показателей надежности. Параграф содержит основные сведения об электрооборудовании СУЗ, необходимые для расчетов, включая ЛФР каналов СУЗ и статистику отказов. Оценка показателей надежности была выполнена для всех элементов электрооборудования СУЗ реактора ЭГП-6.
В табл. 1 приведены результаты расчета средней наработки элементов СУЗ.
Таблица 1.
Наработка на отказ элементов СУЗ реакторов ЭГП-6 БиАЭС Элемент Наработка Элемент Наработка Элемент Наработка на отказ, на отказ, на отказ, год год год УЗС-13 35,4 БПЧ 37,2 БРТЭ АЗ БПР 38,3 БУС 38,1 КПП 49,БП-38 35,3 БДС 37,5 АВ 31,БПДА 12,1 БРВ 33,1 БР РР 35,БПЛ 12,1 БР АЗ 42,4 ПР РС АЗ 17,БУ 12,1 БПМ 35,6 КНК53 44,СУ68 54,7 БРПЭ 41,1 КНК56 40,ПР РР 42 ПР АЗ 36,7 ПР АР 32,Для оценки структурной надежности подсистем АЗ реакторов ЭГП-6 БиАЭС были выбраны четыре основные подсистемы: канал Выработка сигнала АЗ в случае превышения нейтронной мощности на 20% (СПМ); формирование сигнала АЗ Снижение периода реактивности менее 20 сек (СПР); исполнительной части АЗ (СИАЗ); автоматический регулятор (АР). Был выполнен анализ структуры данных подсистем, построены ЛФР, выполнен анализ чувствительности и значимости элементов. Согласно ГОСТ 26843-86, произведена оценка вероятности несрабатывания на требование. Полученные характеристики структурной надежности удовлетворяют требованиям ГОСТ 26843-86 и равны: СПМ - 2,5*10-9; СПР - 4,6*10-9;
АР - 3,9*10-5; СИАЗ - 3,1*10-11.
Анализ результатов оценки ресурсных характеристик элементной базы и структурной надежности подсистем СУЗ ЭГП-6 позволил сделать заключение о том, что при сложившимся на БиАЭС режиме управления ресурсом техническое состояние элементов и систем СУЗ соответствует требованиям нормативов, что подтверждает обоснованность продления установленного срока службы электрооборудования СУЗ БиАЭС до 31.12.2016 (Решение ОАО Концерн Росэнергоатом БИЛАЭС14Р-701К/04-07-2011 от 23.12.2011).
Следует отметить, что результаты расчетов являются оценочными и получены на основании статистической информации малого объема. При поступлении новых данных оценки будут изменяться и уточняться. Поэтому работы по анализу характеристик надежности и остаточного ресурса необходимо проводить периодически, например, раз в год, особенно после проведения любых работ по модернизации, реконструкции и усовершенствованию оборудования СУЗ.
В заключении отмечено, что в представленной работе был разработан системный подход к анализу ресурсных характеристик оборудования АЭС, в основе которого лежат следующие результаты проведенных исследований:
1. Сформулирована проблема необходимости разработки методов оценки ресурсных характеристик, предоставляющих возможность адекватного и достоверного определения технического состояния элементов и систем ЯЭУ. В настоящее время оценка технического состояния оборудования производится по определяющему показателю. В большинстве случаев анализируются показатели безопасности или комплексные показатели надежности, в силу того, что предлагаемые методы оценки ресурсных характеристик как показателей долговечности часто невозможно применить на практике из-за сложности теоретических математических моделей или отсутствия необходимого количества исходной информации.
2. Разработаны методы анализа технического состояния невосстанавливаемых и восстанавливаемых элементов ЯЭУ, основанные на построении неасимптотических моделей расчета ресурсных характеристик с учетом особенностей представления исходных данных об отказах оборудования на российских АЭС. Для невосстанавливаемых элементов разработан метод расчета неасимптотической оценки математического ожидания остаточной наработки до отказа. Разработаны неасимптотические модели для оценки математического ожидания прямого остаточного времени невосстанавливаемого и восстанавливаемого оборудования.
3. Разработан метод построения доверительных интервалов для оценок математического ожидания прямого остаточного времени с использованием модифицированной бутстреп процедуры. Дано математическое обоснование применения представленного подхода для анализа точности оценок исследуемого показателя.
Приведены примеры доверительного оценивания математического ожидания прямого остаточного времени на тестовых данных.
4. Продемонстрирована возможность использования разработанных неасимптотических моделей для расчета ресурсных характеристик на примере электрооборудования системы управления защитой реакторов РБМК-1000 САЭС и ЭГП-БиАЭС. Результаты исследований включены в методики оценки технического состояния и остаточного ресурса электрооборудования СУЗ реакторов ВВЭР (РД ЭО 0294Ц01 (изменение №1), а также использовались при оценке текущего состояния электрооборудования СУЗ реакторов РБМК-1000 САЭС и ЭГП-6 БиАЭС.
Основные публикации по теме диссертации В рецензируемых журналах из списка ВАК:
1. Соколов С.В. Оценка остаточного ресурса невосстанавливаемых элементов электрооборудования СУЗ реактора РБМК-1000 первого блока Смоленской АЭС / С.В. Соколов, А.В. Антонов, В.А. Чепурко // Известия ВУЗов. Ядерная энергетика.
Ч2007. Ч№ 3. ЧВыпуск 1. ЧС.38-43.
2. Соколов С.В. Оценка остаточного ресурса подсистем СУЗ реактора РБМК1000 первого блока Смоленской АЭС / С.В. Соколов // Известия ВУЗов. Ядерная энергетика. Ч 2009. Ч№ 3. ЧС.37-43.
3. Антонов А.В. К вопросу оценки остаточной наработки восстанавливаемого оборудования ядерных энергетических установок / А.В. Антонов, С.В. Соколов, В.А. Чепурко // Надежность. Ч2011. Ч№ 4. ЧС.2-13.
В других изданиях:
4. Соколов С.В. Разработка методов поддержки принятия решений по определению остаточного ресурса оборудования АЭС / С.В. Соколов, А.В. Антонов // Системный анализ в проектировании и управлении: труды 9-й международной научно-практической конференции, Санкт-Петербург, 2005. ЧС.264-266.
5. Соколов С.В. Об одном методе оценки прямого остаточного времени восстанавливаемых систем / С.В. Соколов, А.В. Антонов // Математические методы в теории надежности. Теория. Методы. Приложения.: Расширенные тезисы докладов 6-й международной конференции MMR2009, Москва, 2009. ЧС.534-537.
6. Sokolov S.V. The model of residual lifetime estimation for NPPs equipmentТs operation strategy with heeding preventive maintenance / S.V. Sokolov, A.V. Antonov, V.A. Chepurko // Mathematical Methods in Reliability: Proceedings of International Conference MMR2007, Glasgow, UK, 2007. ЧP.216-218.
7. Antonov A. Assessment of residual lifetime of NPP equipment based on operational information specific type / A. Antonov, S. Sokolov // Accelerated Life Testing:
Proc. of the Third International Conference ALTТ2010, Clermont-Ferrand, France, 2010.
ЧP. 85-90.
8. Antonov A. Mathematical model of the residual lifetime of NPP equipment calculation based on operational information specific type / A. Antonov, S. Sokolov, V.
Chepurko // Applied Methods of Statistical Analysis. Simulation and Statistical Inference: Proceedings of the International Workshop AMSAТ2011, Novosibirsk, Russia, 2011. ЧP.54-62.
Компьютерная верстка С.В. Соколов ЛР № 020713 от 27.04.19Подписано к печати 30.04.2012 Формат бум. 60x84/Печать ризограф. Бумага МВ Печ. л. Заказ № 3040 Тираж 100 экз. Цена договорная Отдел множительной техники ИАТЭ 249035, г.Обнинск, Студгородок
Авторефераты по всем темам >>
Авторефераты по разным специальностям