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

Ключевые слова:
вероятностный анализ безопасности, вероятностные модели отказов, вероятностные модели восстановлений
Текст
Текст произведения (PDF): Читать Скачать

Введение

В соответствии с [1] понятие (или термин) «Безопасность АС» определено как «свойство АС при нормальной эксплуатации и нарушениях нормальной эксплуатации, включая аварии, ограничивать радиационное воздействие на персонал, население и окружающую среду установленными пределами».

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

По отношению к АС причинение вреда связано с радиационным воздействием. Поэтому приведенные в [1] определения понятия безопасности можно считать эквивалентными.

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

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

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

 

1. Независимые отказы оборудования

 

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

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

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

 

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

 

Параметры формул

Описание

Q

Вероятность отказа на требование

λ

Интенсивность отказов

F

Частота

W

Параметр потока отказов

µ

Параметр потока (частота) восстановлений

TR

Среднее время восстановления

TI

Интервал между испытаниями элемента

TF

Время первой проверки

TM

Время работы

 

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

 

1.1. Постоянно контролируемый, восстанавливаемый элемент (тип 1)

 

Модель предполагает экспоненциальное распределение, как для интенсивности отказов, так и для восстановления, т.е. интенсивность отказов и время восстановления являются константами. Неготовность  элемента такого типа рассчитывается по формуле

.                                                (1)

Требуемые параметры: λ, µ(r,TR).

Необязательные параметры: q.

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

Первое слагаемое является необязательным и в большинстве случаев не используется. Если q > 0, то неготовность имеет постоянную составляющую, и в нулевой момент времени имеет значение q, а затем понижается асимптотически до 0 со скоростью µ.

Среднее значение неготовности (при t → ∞) элемента типа 1 вычисляется по формуле

.                                                                        (2)

Частота W(t) для элементов типа 1 вычисляется по формуле

.                                                               (3)

При низкой неготовности частота приблизительно равна интенсивности отказов.

 

1.2. Периодически проверяемый элемент (тип 2)

 

Требуемые параметры: λ, TI(r,TI). .

Необязательные параметры: q, TR, TF.

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

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

,                                                    (4)

Функция неготовности в этом случае отображается пилообразной кривой.

Если задано время первой проверки TF>0, то неготовность будет рассчитываться так же по формуле (4), но моменты проведения проверок будут иметь другие значения:

Среднее значение неготовности элемента типа 2 получается интегрированием Q(t) по отрезку времени равному тестовому интервалу

.                                        (5)

Такая модель справедлива, если время восстановления незначительно (или TR=0). Если же это не так, то следует учитывать время восстановления (необязательный параметр TR).

В программе приняты следующие допущения:

  • восстановление предпринимается сразу после проведения теста,
  • восстановление происходит за фиксированное время TR.

Неготовность элемента типа 2 с учетом восстановления вычисляется по формулам

 

 для ,                                                    ….(6)

 для , n = 1, 2, 3,

 

      для ,

                               для .

Средняя неготовность с учетом восстановления вычисляется по формуле

.                                                           (7)

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

 

                                        для ,

                                     для ,

           для ,           (8)

                                                          для ,

                                  (9)

Частота в этом случае вычисляется по обычной формуле

                                                                                        (10)

 

1.3. Постоянная во времени неготовность, отказ на требование (тип 3)

 

Это самая простая и очень часто используемая модель, использующая единственный параметр q - вероятность отказа на требование.

В этом случае используются формулы

                                                                                                        (11)

 

1.4. Элемент с фиксированным временем работы (тип 4)

 

Неготовность в этом случае, так же как и для типа 3 является постоянной величиной. Однако ее значение не задается в исходных данных, а рассчитывается по заданным значениям интенсивности отказов и времени работы. Кроме того, неготовность может иметь и составляющую q - вероятность отказов типа на требование.

Требуемые параметры: λ, TM.

Необязательный параметр: q.

Используются следующие формулы

                                                                                       (12)

1.5. Постоянная частота (тип 5)

 

Эта модель используется, когда событие хорошо описывается пуассоновским

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

                                                                                                        (13)

Модель используется только для исходных событий, которые включаются в разрабатываемые деревья событий (рис. 1), а также могут использоваться в деревьях отказов (рис. 2).

1.6. Невосстанавливаемый элемент (тип 6)

 

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

Требуемый параметр: l (r).

Необязательный параметр: q.

                                                                                         (14)

 

2. Отказы общего вида

 

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

При расчете используется три параметрические модели учёта отказов по общим причинам: бета фактор, греческие буквы, альфа факторы.

Работа с группами ООВ может быть легко «прослежена» разработчиком модели. После определения группы определяются все возможные события комбинаций отказов по общим причинам и разрабатывает для каждого независимого отказа дополнительное дерево отказов по общим причинам. В этом дереве отказов логическим оператором OR (“или”) первичное событие независимого отказа объединяется со всеми возможными событиями комбинаций ООВ. И дополнительные события ООВ, и деревья отказов ООВ программа автоматически вводит в базу данных проекта, к которой имеется доступ для разработчика. Таким образом, разработчик модели может просмотреть логику деревьев отказов ООВ и численных значений параметров надежности отдельных дополнительных событий ООВ.

Три модели и соответствующие параметры, используемые для моделирования отказов по общим причинам приведены в таблице 2.

 

Таблица 2. Три модели и соответствующие параметры, используемые для моделирования отказов по общим причинам

 

Модель отказов по общим причинам

Параметры

Модель бета-фактора

Бета, â

Модель греческих букв

Бета, â
Гамма, â
Дельта, ä

Модель альфа-фактора

Альфа 2, á2
Альфа 3, á3
Альфа 4, á4

 

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

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

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

N - Количество первичных событий, включенных в одну группу отказов по общей

причине.

Таблица 3. Формулы, по которым определяются показатели неготовности для событий отказов по общим причинам

 

Модель бета-фактора

Модель греческих букв

Модель альфа-фактора

Обобщенные формулы

, k =1

 , 1<k<N

 , k = N

ρ1 = 1, ρ2 = β, ρ3 = γ,…, ρN+1 = 0,

,

 

Формулы для четырехэлементных групп отказов по общим причинам

 

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

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

 

Благодарности

 

Работа выполнена и опубликована при поддержке РФФИ, гранты РФФИ 17-07-01475, 19-07-00455 и 20-07-00577.

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

1. Федеральные нормы и правила в области использования атомной энергии «Общие положения обеспечения безопасности атомных станций» (НП-001-15): [Нормы и правила НП-001-15: утверждены приказом Федеральной службы по экологическому, технологическому и атомному надзору от 17 декабря 2015 г. № 522]. - М., 2015. - 74 с.

2. РД 03-418-01. Методические указания по проведению анализа риска опасных производственных объектов. Утверждены постановлением Госгортехнадзора России от 10.07.01 № 30. Введены в действие с 01.09.01.

3. NUREG/CR-2300, «PRA Procedures Guide», US NRC, January 1983.

4. NUREG/CR-2815, «Probabilistic Safety Analysis Procedures Guide», US NRC, August 1985.

5. NUREG/CR-4550. Analysis of Core Damage Frequency from Internal Events: Methodology Guidelines. Volume 1. US NRC, September 1987.

6. M.A. Berberova, A.V. Dmitriev, A.V. Golubkov, A.I. Elizarov, «Calculation of the probabilistic safety analysis and reliability by the fault trees and event trees methods», Graphi Con 2019. The 29th International Conference on Computer Graphics and Vision. Conference Proceedings (2019), в печати.