CC BY
155
УДК 004.492.2
ОЦЕНКА РИСКОВ ФИЗИЧЕСКОЙ БЕЗОПАСНОСТИ ЛИНЕЙНОЙ ЧАСТИ МАГИСТРАЛЬНОГО НЕФТЕПРОВОДА
Кукало Иван Анатольевич,
аспирант кафедры радиоэлектроники и защиты информации ФГБОУ ВПО «Томский государственный университет систем управления и радиоэлектроники», Россия, 634050, г. Томск, пр. Ленина, 40. E-mail: i@kukalo.ru
Гривцов Сергей Николаевич,
ведущий специалист отдела охраны магистральных нефтепроводов Службы безопасности ОАО «Центрсибнефтепровод», Россия, 634050,г. Томск, ул. Набережная р. Ушайки, д. 24. E-mail: GrivtsovSN@tom.transneft.ru
Актуальность работы обусловлена необходимостью защиты объектов топливно-энергетического комплекса от совершения актов незаконного вмешательства.
Цель работы: обобщение и развитие методики оценки рисков в иерархических структурах критически важных объектов для пространственно-распределенных линейных объектов топливно-энергетического комплекса типа линейная часть магистрального нефтепровода.
Методы исследования: Математическое моделирование распределения вероятности рискового события, экспертная оценка методом Дельфи, методом ранжирования, метод непосредственной оценки, численное интегрирование двумерной функции риска во времени и пространстве, интерполяция функции траектории проведенного мероприятия.
Результаты: Предложена методика оценки риска физической безопасности для линейной части магистрального нефтепровода, реализующая концепцию ГОСТР ИСО 31000-2010. Проведена идентификация рисков физической безопасности линейной части магистрального нефтепровода, разработана математическая модель оценки рисков для пространственно-распределенных линейных объектов топливно-энергетического комплекса типа линейная часть магистрального нефтепровода для проведения анализа и сравнительной оценки риска. Модель позволяет рассчитать пространственно-временное распределение величины риска физической безопасности, которое основывается на предположении о двумерном дискретном случайном распределении вероятности рискового события. Результирующее значение величины риска физической безопасности определяется на основе данных о проведенных мероприятиях по контролю защищенности линейного объекта. Произведен расчет выходных параметров модели для участков магистрального нефтепровода «Александровское~Анжеро-Судженск» с учетом входных параметров, предоставленных ОАО «Центрсибнефтепровод».
Ключевые слова:
Оценка рисков, магистральный нефтепровод, физическая безопасность, матрица распределения риска, пространственно-временное распределение, накопленный риск, уровень риска.
Основные меры по обеспечению безопасности линейных объектов топливно-энергетического комплекса (ТЭК) в Российской Федерации устанавливаются Федеральным законом [1]. Известные подходы к управлению безопасностью систем масштаба ТЭК в значительной мере базируются на процедуре категорирования опасных объектов, в результате которой объектам одной категории предъявляются единые требования по обеспечению безопасности [2]. Альтернативой процедуре категорирования, как известно, является подход по управлению рисками систем физической безопасности (УР-СФБ). Создание методик, реализующих данный подход, сопряжено с необходимостью решения нескольких взаимосвязанных задач, включающих установление области применения, оценку риска и обработку риска [3]. Пример формальной реализации такого подхода для иерархических структур критически важных пространственно-локализованных объектов дан в работе [4]. В работах [5-7] предложены аналогичные концепции управления рисками различных систем -
информационных, кибернетических, транспортных, таможенных, банковских и др.
К сожалению, указанные методики управления рисками не могут охватить всего многообразия действующих систем физической безопасности (ФБ). Одним из факторов, ограничивающих метод оценки рисков [4], является допущение о точечном характере объекта УР-СФБ. Это допущение исключает из рассмотрения множество пространственно-распределенных систем. К таковым, в частности, относится линейная часть магистрального нефтепровода (ЛЧ МН), характеризуемая пространственно-протяженной и многосвязной инфраструктурой [8]. Особенности ЛЧ МН осложняют задачу обеспечения безопасности ТЭК, не позволяя использовать для этой цели стандартные методики УР-СФБ, ориентированные на точечные объекты.
Целью настоящей работы является обобщение и развитие методики оценки рисков в иерархических структурах критически важных объектов [4] для пространственно-распределенных линейных
объектов ТЭК типа ЛЧ МН, учитывая зарубежные [9-11] и отечественные [12, 13] разработки по управлению эксплуатационными рисками МН.
Процедура оценки риска, как известно, является ключевым этапом управления рисками и представляет собой процесс, объединяющий идентификацию, анализ и сравнительную оценку риска
[14]. Данная процедура основывается на моделях злоумышленника и угроз СФБ, которые для объектов ЛЧ МН были описаны в работе [15]. Таким образом, разрабатываемая методика является основой для создания специализированной модели управления рисками линейных объектов ТЭК, которая позволит реализовать заключительный этап процесса управления рисками - обработку риска ФБ [3].
На первом этапе оценки риска системы ЛЧ МН рассмотрим перечень источников риска и событий. Считаем, что множество источников рисковых событий для ЛЧ МН совпадает с множеством угроз
[15].
Возможный перечень источников рисковых событий для типового линейного объекта - ЛЧ МН, представлен в табл. 1.
Таблица 1. Перечень источников рисковых событий
п Источник рискового события Описание источника
п Группа преступных элементов Лицо или группа лиц, действующие согласованно с целью совершения акта незаконного вмешательства для получения выгоды
П Внутренние нарушите- ли Лицо или группа лиц, обладающих правом доступа на объект в силу выполнения служебных или иных обязанностей
п-3 Диверси-онно-тер-рористиче-ская группа Подразделение, используемое для совершения диверсионно-террористических актов уничтожения или временного выведения из строя важнейших объектов топливно-энергетического комплекса
Анализ известных угроз [15] и рисковых событий позволил выявить множество возможных рисковых событий для линейного объекта 5 как т$. Отметим, что для ЛЧ МН каждому рисковому событию соответствует одна угроза безопасности.
Рассмотрим возможный вариант задания элементов множества рисковых событий, представленный в табл. 2.
Перейдем далее к анализу риска, включающего оценку вероятности и последствий идентифицированных опасных событий с учетом наличия и эффективности применяемых способов управления [14]. Для учета особенностей пространственно-протяженных линейных объектов, руководствуясь общими принципами моделирования [16], опишем концептуальную модель оценки риска. С этой целью выделим основные понятия:
• акт незаконного вмешательства - противоправное действие (бездействие), угрожающее безопасному функционированию линейного объек-
та, повлекшее за собой повреждение или уничтожение имущества либо создавшее угрозу наступления таких последствий [1];
• риск ФБ - количественная оценка сочетания вероятности и последствий совершения акта незаконного вмешательства [3];
• линейный объект - линия электропередачи, линия связи (в том числе линейно-кабельные сооружения), трубопровод, автомобильная дорога, железнодорожная линия и другие подобные сооружения [17];
• контрольные мероприятия - мероприятия по контролю защищенности линейного объекта, направленные на выявление совершенных актов незаконного вмешательства.
Таблица 2. Перечень рисковых событий системы 5
Г Наименование рискового события Источники рисковых событий Описание угрозы
г Хищение нефти с использованием действующей врезки в ЛЧ МН либо регулярный отбор нефти через вдоль-трассовое оборудование п, п-2 Установка патрубка с задвижкой с последующим нарушением целостности трубопровода для несанкционированной кражи нефти, взлом вантузов, камер приема-пуска средств очистки и диагностики, запорной арматуры и др.технологических устройств с целью несанкционированной кражи нефти
Г-2 Кража вдоль- трассового оборудования п, п-2 Кража вдольтрассового оборудования и его компонентов (высоковольтные линии,компоненты из цветных металлов и др.)
Гз Установка взрывчатого устройства на ЛЧ МН Пз Реализация террористического акта с целью нанесения материального, экологического и другого ущерба объекту топливно-энергетического комплекса
Далее, следуя [16], в качестве входных параметров математической модели принимаем параметры линейного объекта й, интенсивность инцидентов нарушения ФБ линейного объекта, характеристики проводимых контрольных мероприятий.
Важной характеристикой рассматриваемой модели является ее динамический характер, учитывающий суточные и сезонные колебания риска вдоль различных участков ЛЧ МН.
Расчет этой динамики основывается на априорном показателе интенсивности инцидентов нарушения ФБ и степени опасности для того или иного участка линейного объекта. Модель должна также учитывать информацию о мероприятиях по контролю защищенности линейного объекта, снижающих величину риска ФБ в местах их проведения на величину их рископонижающего потенциала [4]. Выходным параметром модели управления рисками й является динамика пространственновременного распределения величины риска ФБ ЛЧ МН. Аналогичные подходы по использованию про-
странственно-временного распределения вероятности событий применяются в информационной безопасности [18], экологии [19], криминалистике [20] и др.
Каждое рисковое событие имеет три основные характеристики: ср - цену риска (оценку ущерба, который может быть нанесен событием риска), рр -вероятность рискового события и величину риска wp, рассчитываемую по формуле [4]
^ = с * • Р* ■
На основе входного параметра модели Я - периода повторяемости инцидентов нарушения физической безопасности линейного объекта - зададим границы области определения функции рр во времени teЛte(tнач,tкон). Оценку величины риска для каждого события определим как функцию
(1) = с,* • р,* (1)> (1)
где wf(t) - значение величины риска ФБ в момент времени Р, ср - цена риска, т. е. оценка ущерба, который может быть нанесен событием риска; р^) -значение вероятности реализации рискового события в момент времени t.
Определим линейный объект как отрезок числовой прямой 5(хнач,хкон), где хнач - начальная граничная точка отрезка, хкон - конечная граничная точка отрезка. Зададим границы области определения функции рр вдоль оси линейного объекта как хеЛхе(хнач,хкон). Следуя подходу авторов [12], добавим второй аргумент в формулу (1) и определим функцию Wf(x,t)
w.R (х, 1) = с.* • р.* (х, 1 )■
Рассмотрим рисковое событие как двумерную дискретную случайную величину, определенную на равных интервалах длиной ht и Нх. Тогда интегральное значение риска на интервале ЛtxЛx будет определяться уравнением
w .* (Лх, Л1) = с,*р.* (х, 1) =
п т __ __ п т __
= ЕЕ w,*(х, о)=с * ЕЕ р, *(х, ь)’
7=1 ,=1 7=1 ,=1
х — х 1 1
пттл 7„ кон нач 7„ кон нач „
где пх =----------------------------------------, ц =-шаги сетки; п -
т п
количество шагов сетки для Л^ т - количество
шагов сетки для Лх; з^м - множество индексов,
определенных на интервале Л^ £е7м - множество
индексов, определенных на интервале Лх;
- e((/-1)ht,^ht), ^((/-^хДх) - интервалы;-^-,-) -
значение величины-риска на интервале -,
Значения w1—x1,-j) могут быть записаны в виде матрицы ^»(х";,-)|| (табл. 3).
В качестве допущения модели, значение цены риска линейного объекта считаем постоянной величиной, не зависящей от времени и ее конкретного расположения на линейном объекте. Для линейных объектов транспортной инфраструктуры цену риска ср (величину ущерба) можно определить как обобщающий показатель некоторых частных зна-
чений ущерба физической безопасности ЛЧ МН для рискового события [2] по следующей формуле
С■ * С-К + С-К + С*& + ’
і /л /гр і бал і экол
где врл - финансовый ущерб, определяемый численностью погибших и пострадавших; с,«гр - финансовый ущерб от кражи транспортируемого сырья и уменьшения грузопотока; врбал - балансовая стоимость сооружения или стоимость восстановления; врэкол - стоимостное выражение ожидаемого экологического ущерба.
Таблица 3. Матрица распределения рисков \\м’і^(хі,і)\\
А/ Аг
гі г Хт
А у р (хі, Ь) у р (г, іі) уі* (Хт > І1)
Ь у р ( гі, і3) ™,-р (г, 17 ) ( Хт > {7 )
Іп у р (г1> Іп ) у Р ( г, Іп ) у ,р ( гт > Іп )
Оп-ределим матрицу распределения величины |р;»(х;,-)|| (табл. 4), на основе которой можно вычислить Мх;,-)||
^(х,, ^)
Таблица 4. Матрица распределения \ р‘(хі,і])\\
м Аг
хі г Хт
н Рі? (гі, ?і) Рі* (гі > ?1) Рі* (Хт > І1)
Ті Рік (гі, ) Р,к (гі, ) Рі* (Хт > <7 )
К Рі* (г1> <п ) Р,і (Хі> <п ) Рі* (Хт > <п )
Примем, что события Р;»(-) и р^х) независимы, тогда
р.* (х, г3) = р.* ^) • р., (х, ).
Указанное допущение свидетельствует о том, что рисковое соб-ытие может произойти в любой момент времени - на произвольном участке линейного объекта х и корреляция между функциями вероятности р^-) и рр(-) отсутствует. Соответственно, для вероятности справедливы формулы
__ т _ __ __ к __ _
р,*(17) = Ер*(х,ьх р*(х)=Ер*(х, 1).
,=1 7 =1
По условию нормировки
п т _ __ п __ т __
Е Е р ,*(х, )=Е р,* (17)=Е р* (х)=1 (2)
7=1 (=1 7=1 7 =1
= С,
Пусть в качестве допущения модели величина £ имеет равномерное дискретное распределение. Тогда значение вероятности составляющей р^х) системы случайных величин р^х-,-) будет определяться уравнением
1
Р* (17 ) = _ ■
' 7 п
Вероятность составляющей р^х-) системы случайных величин р^,-) определим экспертным методом [15]. Для всех участков л-е(х1,^,-т) определим среднее значение вероятности реализации рискового события р^Х.р)
р.* (хср) = — ■
' р т
Для того чтобы учесть факторы, влияющие на вероятность реализации угрозы на том или ином участке линейного объекта, введем нормирующие коэффициенты 5, определяющие степень опасности конкретного участка ЛЧ МН
p. (xi ) = £г p.,« (xcp).
(З)
В работе [15] приведены группы факторов, факторы и их веса для определения значения 5 типовой модели угроз ЛЧ МН. Факторы и веса для рассмотренных рисковых событий идентичны факторам модели угроз ЛЧ МН, поскольку каждому рисковому событию ставится в соответствие только одна угроза безопасности. Группы факторов и веса рисковых событий представлены в табл. 5.
Таблица 5. Группы факторов влияния
Обозначение и наименование группы факторов Доля группы, pa
ris Г2 Гз
Грі Предпосылки акта незаконного вмешательства 0,20 0,20 0,80
Гр2 Реализация акта незаконного вмешательства 0,40 0,40 0,10
Грз Извлечения выгоды в результате акта незаконного вмешательства 0,40 0,40 0,10
Для определения коэффициентов 5 в формуле
(3) для каждого участка л-е(х:1,^,-т) линейного объекта воспользуемся методом балльной оценки факторов [13, 21]. На первоначальном этапе определим группы факторов влияния Гра с весовыми коэффициентами р, определяющими вероятность их влияния на реализацию угрозы ФБ. В пределах каждой группы Гра определим факторы влияния в количестве Ва. Каждый фактор имеет буквенноцифровое обозначение ЕЛ, где а - номер группы, Ь - номер фактора в группе. Относительный вклад фактора ¥аЬ внутри своей группы учитывается с помощью весового коэффициента (доли) qaЬ.
Расчет коэффициента 5 производился нами с использованием балльной оценочной системы, при которой каждому значению фактора ¥аЬ ставится в соответствие определенное, назначаемое на основании расчета или экспертной оценки, количество баллов ВаЬ (по 10-балльной шкале), отражающее интенсивность его влияния. При рассмотрении конкретного !-го участка линейного объекта последовательно оценивается степень влияния каждого из факторов. Полученные для всех факторов влияния балльные оценки {ВаЬ,ае(1...А),Ье(1...В)} подставляются в формулу
^,
Вср
где В1 - балльная оценка £-го участка линейного объекта, определяемая как:
А В (а)
В =Е ЕРа • ЯаЬ • ВаЬ >
а=1 Ь=1
где Вср - средняя балльная оценка линейного объекта (среднее арифметическое), полученная на основе балльной оценки каждого участка линейного объекта £е(1...т), а т - общее количество рассматриваемых участков линейного объекта.
Факторы влияния и их весовые коэффициенты представлены в табл. 6.
Таблица 6. Факторы влияния
Обозначение и наименование фактора влияния Доля в группе qab
r{ r-S r3
Группа 1: Предпосылки незаконного вмешательства
Fii Криминогенная обстановка 0,20 0,30 0,70
Fi2 Населенность прилегающих территорий 0,60 0,40 0,20
F3 Уровень занятости населения прилегающих территорий 0,20 0,30 0,i0
Группа 2: Реализация незаконного вмешательства
F2i Защищенность объекта 0,i0 0,i0 0,60
F22 Наличие подъездных путей 0,60 0,60 0,i0
F23 Наличие специального оборудования, персонала и транспорта для реализации акта незаконного вмешательства 0,30 0,30 0,30
Группа 3: Извлечения выгоды в результате незаконного вмешательства
Fi Наличие в прилегающих территориях пунктов сбыта похищенной продукции 0,90 0,80 0
F32 Репутационные, экологические и политические последствия 0,i0 0,20 i
В качестве примера в табл. 7 приведены значения отдельных составляющих фактора F31 и соответствующие им балльные оценки B3f, где m - номер составляющей. Итоговая балльная оценка для данного фактора рассчитывается как сумма балльных оценок нижеприведенных пяти составляющих. Если сумма баллов превышает 10, то B31=10.
В работе [22] предложен способ хранения информации о GPS траекториях движущихся объектов в виде двумерного пространственно-временного распределения в системах обработки неструктурированных данных. Следуя данному подходу и переходя от абсолютных GPS координат к относительному позиционированию на ЛЧ МН, представим данные о мероприятиях по обеспечению безопасности объекта S в виде следующего системного множества
М = {т м (г), к. м},
где ішєїш - множество индексов мероприятий; кмєКш - множество типов мероприятий; т1 м(г) - относительная координата мероприятия на ЛЧ МН.
Таблица 7. Балльная оценка фактора В
т Наименование составляющей т Взі(т)
1 В районе пролегания участка нефтепровода находятся потребители, использующие нефть в качестве топлива (малые котельные, теплицы, хлебозаводы, асфальтобитумные заводы и др.) 5
2 В районе пролегания участка нефтепровода находятся нелегальные нефтеперерабатывающие заводы 10
3 В районе пролегания участка нефтепровода находятся частные малые нефтеперерабатывающие заводы 7
4 В районе пролегания участка нефтепровода находятся крупные нефтеперерабатывающие заводы 2
5 В районе пролегания участка нефтепровода осуществляют предпринимательскую деятельность юридические лица, занимающиеся торговлей нефтепродуктами 3
Функция т1 м: м^хэх,м является таблично
заданной (табл. 8).
Таблица 8. Значения функции т()
г гі" '-2 ... ... гіШ 2 ... ... ... л Ч
тДі) х{‘ х2м ... ... Xі ... ... ... х\х
Функция т1 м(£) представляет собой траекторию мероприятия во времени и пространстве линейного объекта, а множество ге(1.Д«) представляет собой индексы координат мероприятия. В качестве допущения определим функцию т1 м(£) как непрерывную.
Таким образом, каждому мероприятию бу-дут соответствовать дискретные значения р;»(л-,-) в области определения ЛххЛ£ (табл. 2-. Для вычисления промежуточных значений (х;,-) для каждого мероприятия т1 м(£) значение функции т1 м(£) аппроксимируем в виде непрерывной кусочно-линейной функции, построенной по значениям функции в двух соседних точках
т. м (г) =
•М ;М ;М
(г - г; )(х2 - х; )
, .-М
(А - г; )
;М ;М ;М
(г - г\ )(х;+; - х; )
^ Iм /М ч
(С; - )
•М -М -М
(г - К м -;)(хп м - хп м -;)
+ х; , г; < г < г2
+ х , г < г < г'
; "> ; ; -
(К, м - гп м -;)
•м -м -м
■ + хп _;, гп _; < г < гп
п.м ; п.м ; п.\
Тогда каждое контрольное мероприятие в каждой точке его проведения снижает вероятность рискового события до момента предыдущего мероприятия на величину
(х/, г,) = Рі* (х, г,) •(; - X
•м .М .М
х/ = х; , гі Є(г; , г; X
.пред
где г;Р
- момент времени проведения предыдуще-
)м".
го мероприятия; (1 , ) - интервал времени с
момента проведения предыдущего мероприятия; стк мб(0,1) - коэффициент эффективности сниже-н^ я рискового события для типа мероприятия к м.
Коэффициент ак м определим с применением методов экспертного оц; енивания путем непосредственного оценивания эффективности снижения рискового события для мероприятия типа к м. Чем выше коэффициент ак м, тем эффективнее мероприятие снижает величину риска на временном интервале до момента проведения предыдущего мероприятия, поскольку в результате контрольного мероприятия мы получаем информацию об отсутствии реализованных рисковых событий на обследуемом участке. Общий рископонижающий потенциал меро-
;м -I
приятия п\К можно вычислить по формуле
.М г—I .м — — . м . м . м
и' =Ер'* (х,tj),х = х ,Ь е(С„,,С ).
С учетом основного свойства закона распределения (2) значение вероятности рискового события перераспределится в «будущие» участки времени и на «соседние» участки линейного объекта
Р ',*' (хі, гі) =
и.* —
Iм
г г
; * * * І
, х = х; , гі
Iм Iм Iм
РК (х; , г; )
Р '* '* (х', ) =
■ -К—^ х+1
х е(х1 ...^:Р х^... хтх ^ = с.
Каждое мероприятие снижает отдельны-е элементы общей матрицы распределения |р;»(х;,-)|| на величину |р;*гМ(х;,-)|| и повышает другие элементы матрицы |р;»(х;,-)|| на величину |р"гМ(х;,-)||. Для за/м /м
данной координаты х2 , ^ мероприятия т1 м(£)
может быть матрица корректировки значения вероятности ||р|(х".,-)|| на всей области определения
функции р;»(х;,-) (табл. 9).
Для оценки рископонижающего потенциала множества мероприятий М={т1 м(0,к; м} необходимо осуществить последовательный расчет значений
|^|(;-,-5)|| для каждого х1 = т.и (1) = х2'м, упорядо-
„ „ чивая множество значений всех мероприятии по возрастанию, для того чтобы достоверно определить дату проведения предыдущего мероприятия.
Таблица 9. Матрица корректировки значения вероятности | \$|(х^-)\\
м Ат
.м х{ .м х 2 хГ хпм п м
Ч 0 0 0 0
Ч пред 3 0 0 Р * * ( х м, і;) 0
ім Ч Пі м р " ( х1 М 4і М ) Рі* (х1 , ‘п1 м) р" (Xім Чім ) Рі* ( 2 п м) ім ім Р ,* ( хі , ^ м) р" (хім 4ім ) ^і* ( Пім 5 Пім )
ім Пі м +1 0 0 р" (х1 м І1 м ) "і* (Хі пм +1) 0
Чп 0 0 0 Р-* (х і ,) 0 0
Если траектории нескольких мероприятий про/ #М , #М\
ходят через точку (х\ ,Ц ), то для мероприятия с на- -М -М
ибольшим ок (*2 , *2 ) > 0, для остальных меро-
приятий ? = ? . Таким образом, при расчете ри-
2пред 2
скообразующего потенциала учитывается наиболее эффективное мероприятие из всех проведенных.
Множество мероприятий М образуют совокупную матрицу корректирующих мероприятий |рм(х,-)||, на основе которой возможно определить результирующее значение риска матрицы распределения рисков ||^рл'з(.х^,—)|| как набора арифметических операций с матрицами:
УГ( х ’ )
где Аі'є(інач,іі) - интервал времени для расчета накопленного риска.
Введенный выше показатель Wf(xi,tj) основан на накопленной вероятности события риска, отражает финансовую составляющую риска и может применяться при обосновании затрат по проведению контрольных мероприятий, а также для оценки необходимого количества денежных средств для обеспечения безопасности линейных объектов [13].
Также введем понятие уровня риска Е^(х,,-), выраженного в долях для участка линейного объекта х на момент времени - как:
____^рез( х, А? ')
Е,(х, і }) = ■ '
= с1« 'IIР“(х,ь)| =
=с!» (|| д*(х, о ^ -1 \р™(х, Ь ф =
(________Iм СМ . м_________3
=с^ (X, *1 ^- рМ (X, *1 ^ |
Результирующая матрица ||^!Т(Х,-)|| определяет динамическое, пространственно-временное распределение риска ФБ. Матрица является источником формирования дополнительных интегральных количественных показателей, участвующих в процессе поддержки принятия решений по распределению имеющихся ресурсов для снижения риска.
В качестве примера такого показателя введем понятие накопленного риска WR(Xt,-), описывающего уровень риска в монетарном выражении для участка линейного объекта х1 на момент времени - как
w^s (х, ^) = ^р;( х, а*') =
= ХХ“( X, = сР Хр* (X , Ь'Х
Т <
І'=1
Т Р
. 7'=1
Ч х >00
7'=1
І '=1
Т (х> *у) Т р? (х> оо
І '=1 ' І'=1
Ер(х1,їі) позволяет оценить динамику риска произвольного участка линейного объекта с учетом проведенных мероприятий. Показатель можно выразить в процентном отношении.
Значения Wi^(xi,tj) и Е^х-,-), определенные на интервале АіхАх, могут быть также определены в виде матриц \^Х,-)[и \\Е^х—\.
Показатели Wf(xi,tj) и Е^,-), на наш взгляд, являются альтернативами действующим показателям эффективности процесса анализа и планирования мероприятий для поддержки заданного уровня защищенности линейного объекта, а именно:
• количество мероприятий на участке х1 на момент времени
• количество дней с последнего_мероприятия на участке х1 на момент времени -.
Показатели W1R(xi,tj) и Е^х,-) позволяют более эффективно решать заключительную задачу
процесса оценки рисков - сравнительную оценку риска, за счет следующих преимуществ:
• учитываются различные степени опасности для участков линейного объекта х1 со стороны злоумышленников;
• учитывается различная степень эффективности множества типов контрольных мероприятий {к«}. В качестве примера использования модели рассмотрим результаты расчета показателей Wf(xi,tj) и Е(Х,-) для одного из МН Томской области. При этом будем использовать следующие входные параметры:
• Я - 1 инцидент/год, ^ач - 01.01.2013, ^он -31.12.2013;
• 5 - МН «Александровское-Анжеро-Судженск» (хнач=0,0; Хкон=818,0 км);
• М - данные о мероприятиях, проведенных экипажем службы безопасности (СБ) ОАО «Центр-сибнефтепровод» за смену с 01.01.2013 20:07 по 02.01.2013 2:42 (табл. 10).
В качестве внутренних параметров определим:
• С^=30 млн р.;
• обследование на автомобиле к=0,3; пешее обследование к=0,4; обследование подъездных путей к МН к=0,1;
• Л(=1 мин и Нх=100 м;
• 6=1.
Аппроксимированный график выходнь_1 х параметров модели Wi^(xi,tj) и ЕДХ-,-) для Х-е(394,...440) км и tp02.01.2013 2:42 представлен на рис. 1. На рисунке представлена зависимость влияния информации о проведенных контрольных мероприятиях на выходной параметр модели - величину риска ФБ.
На рис. 2, 3 представлены графики интегральных значений Wt»(Аx',^i) и Е^Ах'для участка
Таблица 10. Мероприятия, проведенные экипажем СБ за смену
Ax' МН «Александровское-Анжеро-Судженск» длиной 400 км при At=1 год. Начальная и конечная координаты участка не указываются, поскольку реальные значения величин с привязкой к конкретному участку составляют коммерческую тайну ОАО «Центрсибнефтепровод».
На рис. 4, 5 представлены трехмерные графики распределения величин W^x^tj)) и E^xt,tj)) для аналогичного участка Ax' МН «Александровское-Анжеро-Судженск» длиной 400 км при tj&At, At=1 год. Графики указывают на слабые и незащищенные места ЛЧ МН и обеспечивают лиц, осуществляющих принятие решений [10], наглядным представлением о защищенности ЛЧ МН во времени и пространстве.
Для построения указанных графиков использовался трехмерный графопостроитель «con-tourf» пакета Matlab, в котором топология поверхности массивов W^x^tj) и E^x^tj) в трехмерном пространстве изображалась в виде цветового набора двумерных линий постоянного уровня. При этом наглядную информацию об изменении функции дает заливка прямоугольника на плоскости XY (участок ЛЧ МН/время) цветом, зависящим от значения функции в точках плоскости (значение риска/уровня риска). Более светлые области графика плоскости соответствуют областям с наибольшим значением риска (уровня риска) и выделяют наиболее опасные места (во времени для участка ЛЧ МН), для которых вероятность появления рискового события наиболее велика. Темные участки выделяют наименее опасные места, где за счет информации, полученной в ходе проведения контрольных мероприятий, вероятность появления рискового события уменьшается.
Время начала мероприятия, t1" Время окончания мероприятия, ф Начальный участок осмотра, x1" Конечный участок осмотра, xf. Тип мероприятия, к«
01.01.2013 20:07 01.01.2013 20:37 394 398,4 Обследование на автомобиле
01.01.2013 20:37 01.01.2013 20:52 398,5 398,4 Пешее обследование
01.01.2013 20:52 01.01.2013 21:19 398,5 400,4 Обследование на автомобиле
01.01.2013 21:19 01.01.2013 21:53 400,2 400,4 Пешее обследование
01.01.2013 21:53 01.01.2013 22:11 408,9 400,4 Обследование на автомобиле
01.01.2013 22:11 01.01.2013 22:43 408,9 409 Пешее обследование
01.01.2013 22:43 01.01.2013 22:54 413,1 409 Обследование на автомобиле
01.01.2013 22:54 01.01.2013 23:19 413,1 413,2 Пешее обследование
01.01.2013 23:19 01.01.2013 23:42 416,4 413,2 Обследование на автомобиле
01.01.2013 23:42 01.01.2013 23:58 414,4 414,6 Пешее обследование
01.01.2013 23:58 02.01.2013 0:12 416,6 418,2 Обследование на автомобиле
02.01.2013 0:12 02.01.2013 0:27 418,3 418,2 Пешее обследование
02.01.2013 0:27 02.01.2013 0:58 418,3 427,9 Обследование на автомобиле
02.01.2013 0:58 02.01.2013 1:18 428 427,9 Пешее обследование
02.01.2013 1:18 02.01.2013 1:34 428 431 Обследование на автомобиле
02.01.2013 1:34 02.01.2013 1:43 431 431 Обследование подъездных путей к МН
02.01.2013 2:24 02.01.2013 2:35 435 435,2 Пешее обследование
02.01.2013 1:52 02.01.2013 2:03 436,2 436,2 Обследование подъездных путей к МН
02.01.2013 2:03 02.01.2013 2:24 435 437,4 Обследование на автомобиле
02.01.2013 2:35 02.01.2013 2:42 437,4 437,5 Пешее обследование
У часток МН «Апександровское-Анжеро-Судженск», км Рис. 1. Значение уровня риска и накопленного риска для участка МН 394-440 км
Рис. 2. Значение уровня риска на участке МН длиной 400 км в течение года
Величина накопленного риска, руб
12 ООО
. 10 ООО
■л>
Г_>
ё
£ 8 ООО
Н
йч
о
о
§ 6 ООО
■г-
| 4 ООО
к
5
щ
2 ООО
°0 50 100 150 200 250 300 350
Время, дней
Рис. 3. Значение накопленного риска на участке МН длиной 400 км в течение года
100 %
9 0
80
70
60
50
500 1000 1500 2000 2500 3000 3500 4000
Участок МН «Апександровское-Анжеро-Сздженск», м
Рис. 4. Распределение уровня риска на участке МН длиной 400 км
В ОАО «Центрсибнефтепровод» происходит внедрение автоматизированной системы (АС), реализующей указанную математическую модель по оценке величины риска физической безопасности. Среди специализированных АС схожей направленности можно выделить разработки зарубежных компаний по управлению рисками эксплуатации МН [23-25].
Первая очередь автоматизированной системы введена в промышленную эксплуатацию и используется во всех подразделениях СБ ОАО «Центрсибнефтепровод». АС позволяет осуществлять централизованный сбор и хранение данных для автоматизированного формирования входных параметров математической модели. На сегодняшний момент в АС зарегистрирована информация о 19 тыс. обследованных участках МН «Самотлор-Александров-ское», « Александровское-Анжеро-Судженск»,
«Игольско-Таловое-Парабель». Эти данные в совокупности образуют основной параметр модели -множество контрольных мероприятий М. Проведено категорирование МН, на основе которого рассчитаны коэффициенты £, произведена оценка коэффициентов к,м.
Работы по внедрению второй очереди АС находятся на стадии разработки технического проекта. Завершение работ по внедрению АС позволит интегрировать разработанную математическую модель в производственные процессы по охране ЛЧ МН ОАО «Центрсибнефтепровод». Использование ко-
личественных показателей Ж(х^-) и Е^х^-) позволит обеспечить руководящий персонал дополнительной информацией для принятия более эффективных управленческих решений при планировании мероприятий по охране ЛЧ МН.
Выводы
На основе методики оценки рисков в иерархических структурах критически важных объектов
[4], с учетом отдельных элементов методик по управлению эксплуатационными рисками МН [9-13], авторами предложена методика оценки риска ФБ для пространственно-протяженных объектов типа ЛЧ МН, реализующая принципы ГОСТ Р ИСО/МЭК 31010-2011 «Менеджмент риска. Методы оценки риска». В рамках первого этапа оценки проведена идентификация рисков ФБ ЛЧ МН, предложен возможный вариант перечня источников риска и соответствующих им рисковых событий. Для проведения дальнейших этапов по анализу и сравнительной оценке риска разработана математическая модель оценки рисков для пространственно-распределенных линейных объектов ТЭК типа ЛЧ МН.
Модель позволяет рассчитать пространственновременное распределение величины риска ФБ, которое основывается на предположении о двумерном дискретном случайном распределении вероятности рискового события. Учитывается разная степень защищенности участков МН путем использо-
вания метода балльной оценки факторов, влияющих на степень опасности участков МН [15]. Результирующее значение величины риска ФБ определяется на основе данных о проведенных мероприятиях по контролю защищенности линейного объекта. Модель учитывает различную степень эффективности множества типов мероприятий.
В качестве модельной иллюстрации произведен расчет выходных параметров модели W^x-,-) и
СПИСОК ЛИТЕРАТУРЫ
1. Федеральный закон Российской Федерации от 21 июля 2011 г. М 256-ФЗ «О безопасности объектов топливно-энергетического комплекса» // Российская газета - издание Правительства Российской Федерации, официальный публикатор документов. 2011. URL: http://www.rg.ru/2011/07/26/tek-dok.html (дата обращения: 28.01.2014).
2. Кононов А.А., Стиславский А.Б., Цыгичко В.Н. Управление рисками нарушения транспортной безопасности. - М.: АС-Траст, 2008. - 210 с.
3. ГОСТ Р ИСО/МЭК 31000-2010. Менеджмент риска. Принципы и руководство. - М.: Стандартинформ, 2011. - 28 с.
4. Оценка рисков в иерархических структурах критически важных объектов / А.А. Кононов, К.В. Черныш, Д.С. Гуревич, А.К. Поликарпов // Труды ИСА РАН / под ред. Д.С. Черешки-на. - М.: ЛЕНАНД, 2010. - Т. 52. - С. 5-15.
5. Управление кибербезопасностью автоматизированных информационных систем / В.Н. Бекетов, А.А. Кононов, С.Н. Осипов, А.А. Фобьянчук // Труды ИСА РАН. - 2011. - Т. 61. - С. 8-13.
6. Цыгичко В.Н., Кононов А.А., Поликарпов А.К. Концепция системы управления рисками нарушения таможенного законодательства // Труды ИСА РАН. - 2011. - Т. 31. - С. 21-41.
7. Стиславский А.Б., Кононов А.А. Автоматизированная система управления безопасностью объектов транспортной инфраструктуры // Управление рисками и безопасностью: Труды ИСА РАН. - М.: ЛЕНАНД, 2009. - Т. 31. - С. 52-73.
8. Безопасность России. Правовые, социально-экономические и научно-технические аспекты. Анализ риска и проблем безопасности. В 4-х частях. Ч. 4. Научно-методическая база анализа риска и безопасности / науч. руковод. К.В. Фролов. - М.: МГОФ «Знание», 2007. - 864 с.
9. D’Atri M.F., Rodriguez D., Garcia-Martinez R. Improving pipeline risk models by using data mining techniques // 24th World Gas Conference. - Buenos Aires, 2009. URL: http://www. igu.org/html/wgc2009/papers/docs/wgcFinal00663.pdf (дата обращения: 28.01.2014).
10. Selva R., Eng C., Mech F.I. Risk Based Inspection (RBI) Best Practice: The Technical Specification for Ensuring Successful Implementation // 13th International Conference on Pressure Vessel & Piping Technology. - London, 2012. URL: http://www.ppsim-tech.com/downloads/risk_based_inspection_best_practice-the_technical_specification_ron_selva_pp_simtech_keynote_pa-per-1_13th_icpvt_2012_london.pdf (дата обращения: 28.01.2014).
11. Selvik J.T., Scarf P., Aven T. An extended methodology for risk based inspection planning // Reliability: Theory & Applications. - 2011. - М 20. URL: http://gnedenko-forum.org/Jour-nal/2011/012011/RTA_1_2011-12.pdf (дата обращения: 28.01.2014).
12. Оценка риска аварий на магистральных нефтепроводах КТК-Р и БТС / Ю.А. Дадонов и др. // Управление промышленной безопасностью. - 2002. URL: http://programs.safety.ru/BTP/ 2002_6/02_06_2-6.pdf (дата обращения: 28.01.2014).
Е(-,-) для различных участков МН «Алексан-дровское-Анжеро-Судженск» при входных параметрах модели, основанных на реальных данных и произвольных значениях внутренних параметров (реальные значения составляют коммерческую тайну ОАО «Центрсибнефтепровод»).
Выражаем благодарность Геннадию Наумовичу Глазову, ведущему научному сотруднику ЗАО «НПФ «Микран», за полезные консультации при подготовке статьи.
13. Анализ риска магистральных нефтепроводов при обосновании проектных решений, компенсирующих отступления от действующих требований безопасности / Лисанов М.В. и др. // Безопасность Труда в Промышленности. Анализ рисков. -2010.- М3. URL: http://www.safety.ru/downloads/pu-blik_s/analiz_riska/4.6.%20% D0%9B%D0%B8%D1%81%D0%B0 %D0%BD%D0 %BE% D0 %B2.PDF (дата обращения: 28.01.2014).
14. ГОСТ Р ИСО/МЭК 31010-2011. Менеджмент риска. Методы оценки риска. - М.: Стандартинформ, 2011. - 74 с.
15. Кукало И.А. Модель угроз системы физической защиты линейной части магистрального нефтепровода // Известия Томского политехнического университета. - 2013. - Т. 322. -М5. - С. 37-41.
16. Тарасенко Ф.П. Прикладной системный анализ. - М.: КНО-РУС, 2010. - 224 с.
17. Постановление Правительства Российской Федерации от 16 февраля 2008 г. М 87 г. Москва // Российская газета - издание Правительства Российской Федерации, официальный публикатор документов. 2014. URL: http://www.rg.ru/ 2011/07/26/tek-dok.html (дата обращения: 28.01.2014).
18. Jiang G., Cybenko G. Fellow Temporal and Spatial Distributed Event Correlation for Network Security // Proceeding of the 2004 American Control Conference. - Boston, 2004. URL: http://www.nt.ntnu.no/users/skoge/prost/proceedings/acc04/ Papers/0173_ WeM11.4.pdf (дата обращения: 28.01.2014).
19. Fink D., Hochachka W.M. Spatiotemporal exploratory models for broad-scale survey data // The Ecological Society of America. -
2010. - М 20. URL: http://faculty.bscb.cornell.edu/~hoo-ker/STEM.pdf (дата обращения: 28.01.2014).
20. Wang X., Brown D.E. The spatio-temporal modeling for criminal incidents // Security Informatics. 2012. URL: http://www.security-informatics.com/content/1/1/2 (дата обращения: 28.01.2014).
21. Безопасность России. Правовые, социально-экономические и научно-технические аспекты. Анализ рисков и управление безопасностью: методические рекомендации. / науч. руковод. К.В. Фролов. - М.: МГФ «Знание», 2008. - 672 с.
22. Querying Uncertain Spatio-Temporal Data / T. Emrich, H.-P. Kriegel, N. Mamoulis, M. Renz, A. Zufle. URL: http://i.cs.hku.hk/~ni-kos/ICDE12_UncertainST.pdf (дата обращения: 28.01.2014).
23. The Synergi Pipeline software. URL: http://www.dnv.com/ser-vices/software/products/synergi/synergi_pipeline/index.asp (дата обращения: 28.01.2014).
24. Pipeline Risk Controller Version 2.0 for Windows. URL: http://www.gulfpub.com/product.asp? ProductID=2428 (дата обращения: 28.01.2014).
25. The Multi-Analytic Risk Visualization. URL: http://www.dnv. com/moreondnv/research_innovation/programmes/materi-als/sensors/marv.asp (дата обращения: 28.01.2014).
Поступила 03.02.2014 г.
UDC 004.492.2
RISK ASSESSMENT OF LINEAR PART PHYSICAL SAFETY IN THE MAIN OIL PIPELINE
Ivan ^ Kukalo,
Tomsk State University of Control Systems and Radioelectronics, 40, Lenin Avenue, Tomsk, 634050, Russia. E-mail: i@kukalo.ru
Sergey N. Grivtsov,
JSC «Tsentrsibnefteprovod», 24, Naberezhnaya Reki Ushayki Street, Tomsk, 634050, Russia. E-mail: GrivtsovSN@tom.transneft.ru
Relevance of the work is caused by need to protect the objects of the fuel and energy complex from illegal intervention.
The main aim of the study is to generalize and to develop a risk assessment technique in hierarchical structures of crucial objects for spatial distributed linear objects of energy industry of the type linear part of the main oil pipeline.
The methods used in the study: Mathematical modeling of risk event probability distribution, expert assessment by the Delfi method, ranging method, direct assessment method, numerical integration of risk two-dimensional function in time and space, interpolation of function of the held event trajectory
The results: The authors have proposed the technique for assessing physical safety risk for linear part of the main oil pipeline, implementing GOST P ISO 31000-2010 concept. Physical safety risks of linear part of the main oil pipeline were identify; the mathematical model for assessing risks for spatial spread linear objects of energy industry like linear part of the main oil pipeline was developed to carry out the analysis and comparative assessment of risk. The model allows calculating existential distribution of physical safety risk dimension which is based on the assumption of two-dimensional discrete casual distribution of risk event probability. Resultant value of size of physical safety risk dimension is defined on the basis of data on the held events on control of linear object security. The authors calculated the output parameters of the model for the main pipeline «Aleksandrovsky-Anzhero-Sudzhensk» considering input parameters provided by JSC «Tsentrsibnefteprovod».
Key words:
Risk assessment, oil pipeline, physical safety, matrix of risk distribution, spatiotemporal distribution, cumulative risk, level of risk.
REFERENCES
1. Federalny zakon Rossiyskoy Federatsii ot 21iyulya 2011 g. № 256-FZ «O bezopasnosti obektov toplivno-energeticheskogo kompleksa» [Federal Law of 21 July 2011 № 256-FZ «On the safety of the fuel and energy complex»]. Available at: http://www.rg.ru/2011/07/26/tek-dok.html (accessed 28 January 2014).
2. Kononov A.A. Upravlenie riskami narusheniya transportnoy bezopasnosti [Risk management for transport safety violations]. Moscow, AS-Trast Publ., 2008. 210 p.
3. GOST RISO/MEK 31000-2010. Menedzhment riska. Printsipy i rukovodstvo [State Standard 31000-2010. Risk management. Principles and Guidelines on Implementation]. Moscow, Standar-tinform, 2010. 28 p.
4. Kononov A.A., Chernysh K.V., Gurevich D.S., Polikarpov A.K. Otsenka riskov v ierarkhicheskikh strukturakh kriticheski vazh-nykh obektov [Risk assessment in hierarchical structures of critical facilities]. Trudy ISA RAN, 2010, no. 52, pp. 5-15.
5. Beketov V.N., Kononov A.A., Osipov S.N., Fobyanchuk A.A. Upravlenie kiberbezopasnostyu avtomatizirovannykh informat-sionnykh system [Cybersecurity management of automated information systems]. Trudy ISA RAN, 2011, no. 61, pp. 8-13.
6. Tsygichko V.N., Kononov A.A., Polikarpov A.K. Kontseptsiya si-stemy upravleniya riskami narusheniya tamozhennogo zakonoda-telstva [The concept of risk management violations of customs legislation]. Trudy ISA RAN, 2011, no. 31, pp. 21-41.
7. Stislavskiy A.B., Kononov A.A. Avtomatizirovannaya sistema upravleniya bezopasnostyu obektov transportnoy infrastruktury [Automated safety management system of transport infrastructure]. Trudy ISA RAN, 2009, no. 31, pp. 52-73.
8. Frolov K.V. Bezopasnost Rossii. Pravovye, sotsialno-ekonomiches-kie i nauchno-tekhnicheskie aspekty. Analiz riska i problem bezo-pasnosti. Nauchno-metodicheskaya baza analiza riska i bezopas-nosti [Russia’s security. Legal, socio-economic and scientific-technical aspects. Risk analysis and security issues. Scientific
and methodological basis of risk analysis and security]. Moscow, Znanie Publ., 2007. P. 4, 864 p.
9. D’Atri M.F., Rodriguez D., Garcia-Martinez R. Improving pipeline risk models by using data mining techniques. Proc. 24th World Gas Conference. Buenos Aires, 2009. Available at: http:// www.igu.org/html/wgc2009/papers/docs/wgcFinal00663.pdf (accessed 28 January 2014).
10. Selva R., Eng C., Mech F.I. Risk Based Inspection (RBI) Best Practice: The Technical Specification for Ensuring Successful Implementation. Proc. 13th International Conference on Pressure Vessel & Piping Technology. London, 2012. Available at: http://www.ppsimtech.com/downloads/risk_based_inspec-tion_best_practice-the_technical_specification_ron_sel-va_pp_simtech_keynote_paper-1_13th_icpvt_2012_london.pdf (accessed 28 January 2014).
11. Selvik J.T., Scarf P., Aven T. An extended methodology for risk based inspection planning. Reliability: Theory & Applications,
2011, no. 20. Available at: http://gnedenko-forum.org/Jour-nal/2011/ 012011/RTA_ 1_2011-12.pdf (accessed 28 January 2014).
12. Dadonov Yu.A. Otsenka riska avariy na magistralnykh neftepro-
vodakh KTK-R i BTS [Assessment of risk of accidents on the main oil KTK-R and BTS]. Upravlenie promyshlennoy bezopasnostyu, 2002. Available at: http://programs.safety.ru/BTP/
2002_6/02_06_2-6.pdf (accessed 28 January 2014).
13. Lisanov M.V. Analiz riska magistralnykh nefteprovodov pri obosnovanii proektnykh resheniy, kompensiruyushchikh otstu-pleniya ot deystvuyushchikh trebovaniy bezopasnosti [Risk Analysis of main oil pipelines in justification of design decisions to compensate deviations from safety requirements]. Be-zopasnost Truda v Promyshlennosti. Analiz riskov, 2010, no. 3. Available at: http://www.safety.ru/downloads/publik_s/ana-liz_riska/4.6.%20%D0%9B%D0%B8%D1% 81%D0%B0%D0 %BD%D0%BE%D0%B2.PDF (accessed 28 January 2014).
14. GOST RISO/MEK 31010-2011. Menedzhment riska. Metody ot-senki riska [State Standard 31010-2011. Risk management. Risk assessment techniques]. Moscow, Standartinform, 2011. 74 p.
15. Kukalo I.A. Model ugroz sistemy fizicheskoy zashchity lineynoy chasti magistralnogo nefteprovoda [The threat model of the physical protection of the linear part of the main pipeline]. Bulletin of the Tomsk Polytechnic University, 2013, vol. 322, no. 5, pp. 37-41.
16. Tarasenko F.P. Prikladnoy sistemny analiz [Applied System Analysis]. Moscow, KNORUS Publ., 2010. 224 p.
17. Postanovlenie Pravitelstva Rossiyskoy Federatsii ot 16 fevralya 2008 g. № 87 g. Moskva [Government Decree of 16 February 2008 № 87 Moscow]. Available at: http://www.rg.ru/ 2011/07/26/tek-dok.html (accessed 28 January 2014).
18. Jiang G., Cybenko G. Fellow Temporal and Spatial Distributed Event Correlation for Network Security. Proc. of the 2004 American Control Conference. Boston, 2004. Available at: http://www.nt.ntnu.no/users/skoge/prost/proceedings/acc04/ Papers/0173_ WeM11.4.pdf (accessed 28 January 2014).
19. Fink D., Hochachka W.M. Spatiotemporal exploratory models for broad-scale survey data. The Ecological Society of America, 2010, no. 20. Available at: http://faculty.bscb.cornell.edu/~hoo-ker/STEM.pdf (accessed 28 January 2014).
20. Wang X., Brown D.E. The spatio-temporal modeling for criminal incidents. Security Informatics. 2012. Available at: http:// www.security-informatics.com/content/1/1/2 (accessed 28 January 2014).
21. Frolov K.V. Bezopasnost Rossii. Pravovye, sotsialno-ekonomiches-kie i nauchno-tekhnicheskie aspekty. Analiz riskov i upravlenie be-zopasnostyu: meAodicheskie rekomendatsii [Russia’s security. Legal, socio-economic and scientific-technical aspects. Risk analysis and safety management: guidelines]. Moscow, Znanie Publ., 2008. 672 p.
22. Emrich T., Kriegel H.-P., Mamoulis N., Renz M., Zufle A. Querying Uncertain Spatio-Temporal Data. Available at: http:// i.cs.hku.hk/~nikos/ICDE12_UncertainST.pdf (accessed 28 January 2014).
23. The Synergi PipeJine software. Available at: http://www.dnv. com/services/software/products/synergi/synergi_pipeline/in-dex.asp (accessed 28 January 2014).
24. Pipeline Risk Controller Version 2.0 for Windows. Available at: http://www.gulfpub.com/product.asp?ProductID=2428 (accessed 28 January 2014).
25. The Multi-Analytic Risk Visualization. Available at: http:// www.dnv.com/moreondnv/research_innovation/program-mes/materials/sensors/marv.asp (accessed 28 January 2014).
