Методические и вычислительные аспекты значимости риска компонентов сложных систем. Ч. I Текст научной статьи по специальности «Компьютерные и информационные науки»

УДК 631.3

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

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

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

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

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

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

Abstract. Methodical and computing aspects of technogenic hazard significance are considered. Algorithmic complex for the rate computation of different significance types is described. It is targeted for the usage in technogenic hazard automated systems.

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

1. Введение

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

Показатели значимости являются основными критериями оценки техногенной опасности при:

- анализе проектов потенциально опасных систем сложной структуры;

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

- анализе риска реальных ПОО и выявлении слабого звена в их структуре;

- установлении режимов и очередности технических осмотров элементов ПОО;

- планировании и организации мероприятий предупреждения отказов элементов

ПОО.

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

- выполнить систематизированное описание комплекса понятий значимости структурных компонент ПОО и причинных факторов в возникновении отказов сложных систем, которое включает, в том числе, классификацию типов значимости, перечень условий вычисления ПЗ, графическую интерпретацию ПЗ риска;

- описать комплекс алгоритмов вычисления ПЗ элементов в системе сложной структуры.

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

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

2. Основные понятия и показатели значимости

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

Термин «сечение» определяет комбинацию наименьшего числа базисных событий, достаточных для возникновения системного отказа.

Г рупповой отказ дизъюнктивной группы - отказ хотя бы одного из элементов группы.

Примечание. Для компактности изложения будут применяться сокращения. Например, полное наименование «значимости риска элемента системы в возникновении системного отказа» в сокращенном виде будет иметь вид: «значимость элемента в системе».

Понятие значимости включает в себя две составляющие:

- «субъект» - событие, которое влияет на возникновение других событий;

- «объект» - событие, на которое влияет субъект.

Одно и то же событие может быть как субъектом, так и объектом. Например, «групповой отказ» является объектом для «отказа элемента» и субъектом для «отказа системы». Каждой паре «субъект - объект» устанавливается определенный тип понятия значимость. На рис. 1 приведены основные типы значимости, используемые в вероятностном анализе безопасности (ВАБ).

Определение. Значимость (по Бирн-бауму) элемента Ек в системе «С» есть величина [2]:

(Бс.Д)

(Бс.Б)

вк =

ЛРК (ак )

(2.1)

Рис. 1 . Перечень типов понятия значимость

базисного события на элементе Е

йР(ак )

где Р(ак) - вероятность возникновения

Рск (ак ) - вероятность отказа системы, зависящая от Р(ак ); ак е {а. / = 1,п) п - количество элементов системы. Показано [2], что

вк = рк (1) - рк (0),

(2.2)

Р(ак ) = 0 соот-

где Рк (1), Рк (0) - вероятность отказа системы при условиях Р(ак) = 1 ветственно.

В терминах «частота отказа системы» показатель в к характеризует величину уве личения частоты отказа системы при переходе элемента Ек из состояния абсолютной на

дежности в состояние абсолютного отказа. Низкий показатель значимости элемента Ек

(по Бирнбауму) в системе свидетельствует о высоком уровне функционального резервирования этой системы.

Определение. Значимость (по Фусселу-Веселы) элемента EK в системе «С» есть величина [4]:

FVK = [PCK(N) - PCK(0)]/PCK(N), (2.3)

где PC (N) - вероятность отказа системы при номинальном значении вероятности отказа элемента EK .

Показатель FVK характеризует относительное увеличение частоты отказа системы при переходе элемента Ek из состояния абсолютной надежности в «номинальное» состояние (состояние, соответствующее регламентируемым нормам).

Кроме показателей BK, FVK, в ВАБ используются показатели «Коэффициенты» и «Интервалы» изменения риска [5]:

- коэффициент увеличения риска (Risk Increase Ratio), (RIR);

- коэффициент уменьшения риска (Risk Reduction Ratio), (RRR);

- интервал увеличения риска (Risk Increase Interval), (RII);

- интервал уменьшения риска (Risk Reduction Interval), (RRI).

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

RIR = PC (1)/ PC (N), (2.4)

RRR = PC (N)/ PC (0), (2.5)

RII = Pc (1) - Pc (N ^ (2.6)

RRI = PC (N) - PC (0). (2.7)

На рис. 2 представлена графическая интерпретация показателей значимости и их соотношения.

B = RII + RRI. (2.8)

Из подобия треугольников AEH и ADS, с учетом AS = 1, EH = RRI, AH = N, DS = B, следует:

RRI = B x N. (2.9)

Из подобия треугольников AEH и EDQ, с учетом (2.8), AH = N, EQ = LC = 1 — N, DQ = RII, следует:

RII = B(1 - N). (2.10)

Используя выражения (2.3), (2.7), (2.9), получим

FV = (B x N)/PC (N), (2.11)

B = FV x PC (N). (212)

N

Выражения (2.11) и (2.12) формализуют связь между показателями В и FV.

Примечание. На рис. 2 в обозначениях показателей значимостей К-го элемента с целью визуального упрощения не проставлен индекс «К», который следует подразумевать. Как следствие, то же самое сделано в выражениях (2.4) - (2.7). Р(Е) - вероятность отказа элемента Е ;

Рс [ Р( Е )] - вероятность

отказа системы, соответствующая значению вероятности отказа элемента Е ;

N - номинальная вероятность отказа элемента Е .

3. Вычисление значимости элементов системы в возникновении отказа системы

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

а - вычисление статистической значимости;

Ь - вычисление структурной значимости;

с - вычисление номинальной значимости;

ё - вычисление ситуационной значимости.

Далее рассмотрен процесс вычисления значимости для каждого приведенного случая.

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

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

вк = пк /NC-nK /^, (3.1)

Рис. 3. Условия вычисления значимости риска

Рис. 2. Графическая интерпретация показател ей значимости риска

где B’K - статистическая значимость элемента EK (по Бирнбауму);

NC - число всех отказов системы;

nK - число отказов системы при условии отказа элемента Ek ;

nK - число отказов системы при условии безотказности элемента EK .

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

В ВАБ модели дерева отказов (ДО) и дерева событий (ДС) являются базовым средством моделирования возникновения и развития техногенной опасности. Применение моделей ДО и ДС позволяет получить формализованное описание аварии в виде дизъюнктивнонормальной формы (ДНФ), переменными которой являются базисные события [5-7].

C = G({a,}(i = 1П)) = US, (32)

reR

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

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

{a } (i = 1, n) - базисные события на элементах системы E.;

Sr - сечения в ДНФ;

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

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

Pc = Q[ P(aJ, P(a2),..., P(an)], (3.3)

где P(a.) - вероятности отказов элементов Et (i = 1, n);

Pc = Q({P(a,-)} (i = Vn) = 2P(Sr) = X ПP(ar, )• (3 4)

reR reR ir e Ir

где P(a. ) - вероятности базисных событий на элементах E. из сечений Sr;

P(Sr) - вероятности сечений Sr;

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

Описание Процедуры оценки значимости элемента в системе (“Significance Assessment Procedure ”) (SAP)

Процедура SAP содержит следующие шаги.

A. Определение оценок вероятностей базисных событий P(ai) (i = 1, n) на элементах системы.

B. Вычисление значений вероятностной функции PC (3.4) в случаях, когда:

* элемент EK находится в состоянии отказа (P(aK ) = 1);

* элемент EK абсолютно безопасен (P(aK ) = 0);

* вероятность отказа элемента EK имеет номинальное значение (P(aK ) = N) . Обозначим эти значения функции PC соответственно P^K (1); PK (0); PCC (N).

PcK(1) = Q[P(aJ,...,P(aK) = 1,...,P (an)]. (3.5)

РСК(0) = 0Р(А),..,Р(ак) = 0,...,Р (ая)]. (3.6)

РК(М) = б№),..,Р(ак) = М,...,Р (ап)]. (3.7)

С. Вычисление значимостей элемента Ек в системе по Бирнбауму и по Фусселу-

Веселы подстановкой значений Р^К (1), Р^К (0), Р^К (N) в выражение (2.2) и (2.3).

Примечание. Формула (3.4) дает удовлетворительный результат в случае, когда

"Е—Р(а.) £ 10-3, т.е. когда 3^—Р(а.) > 10-3 [8].

В противном случае необходимо использовать точные методы, в частности, метод ортогонализации, который позволяет представить ДНФ G({ai}) в виде совокупности ортогональных (несовместимых) конъюнкций отказов элементов

G«а})=и К,

!еЬ

где К{ - ортогональные конъюнкции базисных событий .

В результате вероятность отказа системы может быть вычислена точно:

Рс = 0({Р(а,)}) = X Р( К),

!еЬ

где Ь - количество ортогональных конъюнкций [3].

При реализации шага «А» процедуры БАР возможны следующие случаи:

- значения Р(а.) (' = 1, п) априори неизвестны;

- P(ai) (' = 1, п) имеют номинальные значения (Ni) ;

- P(ai) определяются на основании описания конкретной ситуации, сложившейся на объекте.

Значимости риска, вычисленные в этих условиях, назовем, соответственно, «структурной», «номинальной» и «ситуационной» значимостями.

3.2.1. Вычисление структурной значимости элемента Ек

Определение. Структурная значимость элемента в системе есть значимость, определяемая исключительно на основании формализованного описания структуры системы (логиковероятностной модели взаимосвязей элементов системы).

Шаг «А» процедуры «БАР» заключается в присвоении вероятностям отказов всех элементов системы значений, равных 0,5.

Шаг «В» процедуры “БАР” реализуется при условии

" — Р'(а ) = 0,5. (3.8)

«Е{(1,п )\К} У '/ ’ У '

Примечание. На шаге «В» процедуры “БАР” предварительно выполняется ортогона-лизация ДНФ отказа системы, и только после этого производятся расчеты.

3.2.2. Вычисление номинальной значимости элемента Ек в системе

На шаге «А» процедуры “БАР” значения Р(а.) (' = 1,п)\ К устанавливаются равными своим соответствующим номинальным значениям (N.), которые получены из стандартов и других документов, нормирующих безопасность.

Шаг «В» процедуры “БАР” реализуется следующим образом:

- значения РСК (1) и РСК (0) ((3.5), (3.6)) вычисляются при условии "■ (7л—Р'(а ) = N.;

'е{(1,п)\К} у ^ '’

- значения РСК(N) (3.7) вычисляются при условии

" — Р (а ) = N.. (3.9)

'Е{(1,п) \ I/ I V /

3.2.3. Вычисление значимости элемента на основании описания конкретной ситуации, сложившейся на объекте

Описание ситуации представляет собой совокупность сведений о состоянии каждого элемента системы, включенного в выражение (3.2). Оно формируется по результатам мониторинга состояния объекта и включает в себя для каждого элемента:

- описание реальных событий, происшедших на элементе системы и влияющих на вероятность отказа системы;

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

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

- текущие значения причинных факторов риска, влияющих на возникновение базисных событий.

Значимость элемента, вычисленного на основе подобных данных, будем называть ситуационной значимостью риска элемента. Рассмотрим реализацию шага «А» процедуры “БАР” для каждого указанного случая.

3.2.3.1. Определение вероятности отказа каждого элемента Ei на основании конкретных событий, происшедших на этом элементе

К моменту фиксации состояния элемента возможны следующие события:

а) отказ элемента;

б) ремонтно-профилактические работы (РПР) на элементе не завершены (неработоспособное состояние);

в) РПР завершены или выполнена замена элемента на новый.

В случаях а и б вероятность базисного события устанавливается равной единице

(Р(а) = 1).

В случае в) вероятность отказа устанавливается равной номинальному значению

(Р(а,) = Nl).

3.2.3.2. Определение вероятности отказа элемента на основании определяющего параметра

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

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

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

3.2.3.3. Оценка вероятности отказа элемента на основании наработки элемента системы

Для многих типов элементов сведения о наработке являются наиболее информативными данными, на основании которых оценивается риск отказа. При этом используются методы теории надежности. По данной проблематике имеются достаточно обстоятельные публикации, в частности, [9-11].

В кратчайшем виде оценка вероятности отказа элемента на заданном временном интервале [т0,т0 + Ат] сводится к реализации следующих этапов.

Этап 1. Построение модели отказа элемента Е.. Данный этап включает в себя выбор

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

позволяет вычислить вероятность отказа элемента Е. до момента, обусловленного наработкой /.

Данный этап выполняется в период подготовки базы данных автоматизированной системы оценки и анализа безопасности.

Этап 2. Определение наработки элемента Е1, соответствующей левой и правой границам заданного прогнозного интервала [т0,т0 + Ат].

Моменту т0 соответствует реальная наработка 10, которая определяется:

- по результатам мониторинга за фактической эксплуатацией анализируемого элемента;

- либо с использованием выражения

где К3 - коэффициент усредненной загрузки элемента Ei.

Моменту т0 + Ат соответствует наработка /0 + Аt, где А/ - прогнозируемая наработка за период Ат. А/ устанавливается на основании экспертных оценок либо с использованием коэффициента усредненной загрузки элемента Е..

Этап 3. Вычисление оценки вероятности отказа элемента Е{ в заданном прогнозном интервале [т0,т + Ат].

Оценка вероятности отказа выполняется с помощью выражения [12]:

(3.10)

(3.11)

А/ = Ат0 х К3.

(3.12)

Р(а,) = [^ (/0 + А/) - ^ (0]/[1 - ^ (/0)],

(3.13)

где /0, А/ - наработки, определенные на Этапе 2;

¥' (ї) - аналитическое представление модели отказов, определенное в виде выражения 3.10 на Этапе 1.

Э.2.3.4. Оценка вероятности отказа элемента системы на основании текущих значений причинных факторов риска

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

Для решения проблемы вычисления оценок вероятностей возникновения базисных событий, возникающих под влиянием совокупности независимых факторов риска, разработан «Метод экспертных оценочных шкал» (МЭОШ), детально изложенный в [13, 14]. Здесь приведем его кратчайшее описание.

В МЕОШ выделены два элемента: формирование базы знаний (БЗ), позволяющей формализовать описание ситуаций, возможных на объекте, и вычисление оценок вероятностей базисных событий. БЗ включает в себя описания факторов и их влияний. Каждый фактор представлен множеством его упорядоченных возможных значений. Влияние фактора представлено функцией влияния (ФВ). ФВ фактора X. на базисное событие а. есть

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

а.») в базе знаний формируется Ь. кортежей следующего вида:

Л хм, Л(x))), (314)

где і - индекс базисного события а. (і = 1, п);

. - индекс фактора X. (. = 1,М);

х. - одно из возможных значений фактора X (I = 1, Ь );

Л' (х1. ) - значение функции влияния фактора X. на вероятность возникновения события а. при условии

(X = х )П (V,,. X, = х?), Р.15)

где X? - значение нормы фактора X .

Подобная база знаний позволяет вычислять оценки вероятности базисного события а., возникшего при совокупном влиянии независимых факторов {X.} (. = 1,М) в ситуации, когда

X = Х1С; X2 = х2с;...; Xм = хМ. (3.16)

Оценка вероятности базисного события а. вычисляется с помощью выражения

М

Р(а,) = 1 -П[1 - Л (хс,)]. (317)

.=1

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

Таким образом, шаг «А» процедуры “БАР” при вычислении ситуационной значимости сводится к последовательному применению действий, изложенных в п.п. 3.2.3.1— 3.2.34.

Шаг «В» процедуры “SAP” выполняется при условии

P ) = Ps (a'),

где Ps (a.) - оценки вероятностей возникновения a., вычисленные из условий контролируемой ситуации (п.п. 3.2.3.1-3.2.3.4).

Примечание. Шаг «С» процедуры “SAP” остается инвариантным во всех описанных случаях 3.2.1-3.2.3.

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

В работе рассмотрены методические и вычислительные аспекты ключевого инструмента анализа риска значимости.

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Барлоу Р. Статистическая теория надежности и испытания на безотказность / Р. Барлоу, Ф. Про-шан. - М.: Наука, 1984. - 328 с.

2. Birnbaum Z.W. On the importance of different components in a multi-component system. Multivariate Analyses - 2 / Birnbaum Z.W. - N.Y.: Academic Press, 1969. - Р. 581 - 592.

3. Рябинин И.А. Надежность и безопасность структурно-сложных систем / Рябинин И. А. - СПб.: Издательство Санкт-Петербургского университета, 2007. - 278 с.

4. Fussel J.B. How to Hand-Calculate system reliability characteristics / J.B. Fussel // IEEE Transactions on Reliability. - 1973. - N 3. - R-24.

5. Можаев А.С. Автоматизированное структурно-логическое моделирование систем / А.С. Можаев.

- СПб.: ВМА им. Н.Г. Кузнецова, 2006. - 577 с.

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

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

- 2009. - № 1. - С. 165 - 172.

8. Risk Spectrum PSA Professional 1.10 Theory Manual/RELCON AB, 1998. - 57 p. - Режим доступу: http: www. ri skspectrum .com.

9. Стрельников В.П. Оценка и прогнозирование надежности электронных элементов и систем //

В.П. Стрельников, А.В. Федухин. - К.: Логос, 2002. - 486 с.

10. ГОСТ 27.005-97. Надежность в технике. Основные положения. Модели отказов. - Введ. 05.12.97, № 732. - Киев: Госком Украины по стандартизации, 1997. - 45 с.

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

В.П. Стрельников // Математичні машини і системи. - 2007. - № 1. - С. 137 - 143.

12. Серебровский А.Н. Методы оценки вероятностей отказов в процессах прогнозирования техногенных чрезвычайных происшествий / А.Н. Серебровский // Математичні машини і системи. -2007. - № 2. - С. 111 - 116.

13. Серебровский А.Н. Об одном методе вероятностного анализа безопасности потенциально опасных объектов / А.Н. Серебровский // Математичні машини і системи. - 2002. - № 1. - С. 41- 48.

14. Serebrovsky A.N. Models and Algorithms of Probabilistic safety assessment of potentially hazardous objects/ A.N. Serebrovsky // Материалы “6-th International Conference on Imformation System Technology and Application”, (Kharkiv, may 23-25, 2007). - Kharkiv, 2007. - P. 127 - 134.

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