Методические и вычислительные аспекты значимости риска компонентов сложных систем. Ч. ІІ

МОДЕЛЮВАННЯ І УПРАВЛІННЯ

УДК 631.3

А.Н. СЕРЕБРОВСКИЙ

МЕТОДИЧЕСКИЕ И ВЫЧИСЛИТЕЛЬНЫЕ АСПЕКТЫ ЗНАЧИМОСТИ РИСКА КОМПОНЕНТОВ СЛОЖНЫХ СИСТЕМ. Ч. П

Анотація. Наведено комплекс алгоритмів обчислення показників значущості групи елементів системи, причинних чинників ризику. Комплекс призначений для використання в автоматизованих системах аналізу техногенного ризику.

Ключові слова: значущість елементів системи, чинники ризику, структурна значущість, Бірнба-ум, Фуссел-Весели.

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

Ключевые слова: значимость элементов системы, факторы риска, структурная значимость, Би-рнбаум, Фуссел-Веселы.

Abstract. Algorithmic complex for the rate computation of system group components significance and casual hazard factors are described. It is targeted for the usage in automated systems for the technological hazard analysis.

Keywords: system components significance, hazard factors, structural significance, Birnbaum, Fussell-Vesely.

1. Введение

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

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

«Вывод оборудования в плановый ремонт есть упреждение отказа оборудования, и, как следствие, упреждение отказа системы. Степень упреждения должна зависеть от значимости оборудования. Чем значительней оборудование влияет на безопасность эксплуатации, тем с большим запасом времени должен назначаться плановый ремонт этого оборудования. Если же оборудование имеет пренебрежимо низкое влияние, то ремонт может проводиться по факту отказа» [2].

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

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

© Серебровский А.Н., 2012

ISSN 1028-9763. Математичні машини і системи, 2012, № 3

Определение основных понятий и показателей значимости описано в п. 2 [1].

В данной работе предполагается, что перед вычислением показателей значимости предварительно была проведена формализация связей между элементами системы. Для этого применяются модели дерева отказов и дерева событий. В результате формализованное описание отказа системы (в частности, аварии) можно получить в виде дизъюнктивнонормальной формы (ДНФ), переменными которой являются отказы элементов системы (базисные события) [3-6].

где С - опасное событие, вызванное отказом системы;

G - логическая функция в виде ДНФ опасного события;

{а.} (/ = 1, п) - базисные события на элементах системы Е{;

£,г - сечения в ДНФ;

Я - множество индексов сечений в ДНФ.

ДНФ дает возможность представить вероятность опасного события «С» как аналитическую функцию, у которой аргументами являются вероятности базисных событий.

где P(ai) - вероятности базисных событий элементов Е1;

Р(£ ) - вероятности сечений £г ;

1г - множество индексов базисных событий из сечения £ .

2. Вычисление значимости риска элемента в возникновении отказа группы элементов

2.1. Группа представляет собой конъюнкцию событий из множества {а.} . = 1, п

Пример. Вычислить значимости по Бирнбауму (ВА ) и Фусселу-Веселы (РУА ) [1] элемента «А», входящего в состав сечения £ = АIВIЕ .

ВА = РА (1) - РА (0) = [1X Р(В) х Р(Б)] - [О х Р(В) х Р(Б)] = В( В) х Р( Б).

В общем виде значимость элемента а в сечении

вычисляется со

4^0.5

гласно выражениям:

Вк = П Р(% ),

4^{0з \к3}

так как

= [ Р/ (N) - Р5 (0)]/ Р5 (N ^

Р5А (0) = 0, то ру/ = 1.

Примечание. Из (2.3) следует, что для оценки значимости риска элемента в возникновении сечений не имеет смысла применять показатель FV.

2.2. Группа представляет собой дизъюнкцию независимых событий из множества

{ai} i = 1 n

Пример. Вычислить значимость элемента «А» дизъюнкции G = A U B U D .

FVG4 = [PGA

FVGA = {[NA + P

FVGA = Na /[Na + P

где NA - номинальная вероятность отказа элемента A.

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

G = {aq} q е Qg вычисляется согласно выражению

FVtK = NK /[NK + £P

q

qe{Qo \ к}

\tK

где N - номинальная вероятность aK .

Bg = PgA

Примечание. Из (2.8) следует, что для оценки значимости элемента в составе дизъюнкции не имеет смысла применять показатель «В».

3. Вычисление значимости группы элементов в системе

3.1. Группа представляет одно из сечений Sk формализованного описания отказа системы «С»

Согласно

FVSK = NSk /

nSk + £ P{st)

qeQR

где Он ={(1, Я) \к};

Я - количество сечений в ДНФ.

= П Р" (а,).

.,^°к

где - номинальная вероятность сечения £к;

Ок - множество индексов базисных событий из сечения £к ; Р" (а{ ) - номинальные вероятности событий а{ ;

р

а,

9 ' 9 6°9

О9 - множество индексов сечения !;

Р(а,,) - вероятности событий а. из сечения !ч.

Значимость по ВСк вычислять не имеет смысла (см. Примечание к разд. 2.2).

3.2. Группа, представляющая собой дизъюнкцию независимых событий в системе «С»

Допустим, что во множестве базисных событий {а.} (/ = 1, п) выделена группа О1 = и а],

дополнением к которой является = {а.} \ О1 . Известны номинальные вероятности событий из О1, которые обозначим "'7. Требуется определить значимость по ЯУ группы О1 в системе С = О1 и О2 .

РУС° = [РОСО1

РсО1

с У ' Р(а )="] 4 1 2) Р(а] )=" У 1У Р(а )=Ы] 4 2' 4 '

Рс°‘

Так как события из группы О1 независимы и составляют дизъюнкцию, то [5]

РР,,,=»](О1)=1 -ГК1 - " )• <35>

] }=!

где "} - номинальная вероятность а] (] е J).

Так как события из группы О2 независимы от событий группы О1 , то вероятность РО2 ) инвариантна вероятностям событий из группы О1, в частности

РР

Очевидно, что РР

Подставив (3.3), (3.4), (3.5), в (3.2) и учитывая (3.6) и (3.7), получим

ру° = [1 -П

jеJ jеJ ]

Данное выражение дает возможность вычислять групповую значимость группы О1 в двух случаях:

- когда для событий группы О2 заданы их номинальные вероятности Р(а1) = "1

I е Ь = {1, п}\ J;

- когда для событий группы О2 заданы их ситуационные вероятности Р(а1) = Р!,, где Р^ - ситуационное значение вероятности события а1 I е Ь .

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

4. Значимость причинных факторов риска

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

4.1. Определение значимости фактора в возникновении отдельного базисного события

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

- экспертное оценивание значимостей;

- оценка на основании статистики базисных событий;

- оценка на основании функций влияния факторов риска.

4.1.1. Экспертная значимость фактора устанавливается на основании экспертных оценок об относительной важности фактора риска для вероятности возникновения конкретного базисного события. При этом возможно как прямое установление оценок значимости, так и с применением метода анализа иерархий (составление матрицы парных сравнений, вычисление собственных значений матрицы парных сравнений [7]).

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

7;

где 7; (X) - статистическая значимость фактора " X";

р;

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

4.1.3. Значимость фактора риска, вычисленная на основе функций влияния фактора на базисное событие

Обозначим значимость фактора Xв возникновении базисного события а, через 7\. Подобная значимость, согласно определению функции влияния, вычисляется с помощью выражения

г;

где /а (хкр) , (х*) - значения функции влияния фактора X на возникновение события

"а", при условии, когда фактор X принимает критическое и нормальное значения соответственно.

Примечание. В случае допущения, что хкр и х* соответствуют правой и левой границам интервала возможных значений фактора X, т.е., когда /а (хкр) = 1, /^ (х*) = 0, выражение (4.2) представляет собой значимость фактора X в возникновении события "а" по Бирн-бауму. Таким образом, значимость фактора по Бирнбауму является частным случаем значимости фактора, выраженной через функции влияния.

4.2. Определение значимости заданного фактора Xк в возникновении отказа системы

Значимость фактора в отказе системы вычисляется при условиях:

- создана модель отказа системы в дизъюнктивно-нормальной форме, переменными которой являются базисные события а. (, = 1,п), а также функция вычисления вероятности отказа системы ((1.1), (1.3));

- создана база знаний, содержащая описание факторов риска X. (, = 1,М) и функций влияния факторов на базисные события

/ (х)(/=іі,),

где х, - одно из возможных значений фактора X.,

I. - количество возможных значений фактора X.;

- описана ситуация, в условиях которой находится объект, в том числе, определены ситуационные значения факторов риска (х,,..., хси).

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

I. Вычисление значений вероятностей каждого базисного события а. (, = 1,п) при двух условиях, когда фактор Xк принимает критическое и нормальное значения.

р'к’

,е((1,М )\к}

р~

,е{(1и )\к}

где Гк (х‘к), /к (х* ) - значения функции влияния фактора Xк на событие аі при критическом и нормальном значениях соответственно;

/, (х(С ) - значение функции влияния фактора X. на событие а., когда X, = хСС;

с

х, - ситуационное значение фактора X,.

II. Вычисление значений вероятностей отказов системы при критическом и нормальном значениях фактора Xк. Для этого используется вероятностная функция

Q[{P (а,) , = 1, п}] (1.3), в которую в качестве аргументов подставляются рк(а.) и р* (а,), вычисленные на шаге I.

рк = Q[{PKP

р,с = Q[{PN

III. Вычисление значимости фактора риска в возникновении отказа системы

7,

5. Заключение

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

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

- показано, что значимость фактора риска по Бирнбауму является частным случаем функций влияния, используемых в Методе экспертных оценочных шкал (п. 3.2.3.4 [1]).

Практическое значение данной работы и работы [1].

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

Предлагаемый комплекс алгоритмов вычисления показателей различных типов значимости может быть использован разработчиками:

- ситуационных центров по чрезвычайным ситуациям;

- автоматизированных систем оценки и анализа техногенного риска потенциально опасных объектов;

- гибридных экспертных систем анализа и предотвращения техногенного риска. СПИСОК ЛИТЕРАТУРЫ

1. Серебровский А.Н. Методические и вычислительные аспекты значимости риска компонентов в сложных системах. Ч. І / А.Н. Серебровский // Математичні машини і системи. - 2012 - № 2. -С.145 - 154.

2. Комаров Ю.А. Методология проведения ремонта по техническому состоянию оборудования АЭС / Ю. А. Комаров // Проблеми безпеки атомних електростанцій і Чорнобиля. - 2010. - Вип. 13.

- С. 48 - 59.

3. ГОСТ Р 51901.13-2005. Менеджмент риска. Анализ дерева неисправностей (IEC б 1025:1990, Fault Tree Analysis). - Введ. 31.05.2005. - М.: Стандартинформ, 2005. - 1б с. (Национальный стандарт Российской Федерации).

4. Вероятностный анализ безопасности атомных станций (ВАБ) [В.В. Бегун, О.В. Горбунов, И.Н. Каденко и др.]. - К.: НТУУ «КПИ», 2000. - 5б8 с.

5. Серебровский А.Н. Алгоритм формирования и минимизации логического представления дерева отказов / А.Н. Серебровский, Л.П. Ситниченко, В.Г. Пилипенко // Математичні машини і системи.

- 2009. - № 1. - С. 1б5 - 172.

6. Культура безопасности на ядерных объектах Украины / [В.В. Бегун, С.В. Бегун, В.В. Литвинов и др.]. - К.: НТУУ «КПИ», 2009. - 38б с.

7. Саати Т.Л. Принятие решений. Метод анализа иерархий / Саати Т.Л. - М.: Радио и связь, 1989. -31б с.

Стаття надійшла до редакції 21.11.2011