Модель оценки устойчивости автоматизированных систем пожаровзрывобезопасности

С. Ю. БУТУЗОВ, д-р техн. наук, доцент, профессор кафедры информационных технологий, Академия ГПС МЧС России

Ю. В. ПРУС, д-р физ.-мат. наук, профессор, заведующий кафедрой управления и экономики Государственной противопожарной службы, Академия ГПС МЧС России

В. Л. СЕМИКОВ, д-р техн. наук, профессор, профессор кафедры управления и экономики Государственной противопожарной службы, Академия ГПС МЧС России

О. В. ЯКОВЛЕВ, д-р техн. наук, доцент, ведущий научный сотрудник, Федеральный исследовательский центр "Информатика и управление" РАН

УДК 681.3.067;614+004

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

Проведен анализ методов оценки надежности и устойчивости функционирования автоматизированных систем управления технологическими процессами (АСУТП) пожаровзрывоопасных предприятий. Показано, что в настоящее время при проектировании данных систем на единой программно-аппаратной базе принимается условие, что физико-технические характеристики элементов систем не меняются во времени. Однако в условиях быстротекущих технологических процессов (взрывы, пожары на промышленных предприятиях) из-за временной деградации, помех, недостатков в сборке и т. д. теоретические показатели эффективности могут снижаться до 25 %. Во избежание данных неточностей в расчетах при проектировании автоматизированных систем пожаровзрывобезопасности (АСПВБ) промышленных предприятий предложена математическая модель оценки устойчивости функционирования этих систем, являющихся частью АСУТП, которая основана на исследованиях в области устойчивости автоматических технических систем Э. Рауса и энтропийных методов А. Дж. Вилсона. Проведено моделирование надежности и устойчивости функционирования каналов передачи информации для различных вариантов построения АСПВБ. Показано, что для оптимизации затрат может быть использован комплексный подход, основанный на использовании единой программно-аппаратной платформы микропроцессорной платформы АСУТП.

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

DOI: 10.18322/PVB.2017.26.11.14-20

Введение

В условиях современного экономического подъема возрастает актуальность надежного функционирования автоматизированных систем управления технологическими процессами (АСУТП) промышленных предприятий, являющихся критически важными для российской экономики.

Автоматизированные системы пожаровзрывобезопасности (АСПВБ) в настоящее время стали играть

ключевую роль в структуре АСУТП. При этом они являются сверхбыстродействующими подсистемами, способными обеспечить детектирование избыточного давления взрывоопасных технологических установок в течение 0,1 мкс. В этом случае возможно осуществить в реальных масштабах времени сброс давления, чтобы защитить технологическое оборудование от разрушения и персонал предприятия от гибели [1].

© Бутузов С. Ю., Пранов Б. М., Прус Ю. В., Семиков В. Л., Яковлев О. В., 2017

Однако соблюдение данного требования увеличивает частоту опроса датчиков давления свыше 107 с-1, частота системной магистрали передачи данных должна составлять более 10 МГц, а для обеспечения надежного детектирования избыточного давления — до 100 МГц. В этом случае актуальным становится вопрос устойчивой работы шины АСПВБ, а именно каналов передачи информации [2].

Целью данной работы является создание методики оценки надежности функционирования каналов передачи информации на основе моделирования величины энтропии АСПВБ. Для достижения данной цели разработана математическая энтропийная модель, основанная на научных работах А. Дж. Вил-сона [3] и Э. Рауса [4].

Теоретический анализ

Анализ научно-технической литературы показал, что в настоящее время отсутствуют системные научные исследования устойчивости функционирования АСПВБ [5-8].

Для расчетов АСУТП используются модели, представленные в [9-11]. Данные расчетные модели базируются на условии, что технические характеристики элементов автоматизированной системы не меняются во времени, вне зависимости от воздействующих факторов. Например, такой параметр, как скорость передачи данных, который является основным лимитирующим показателем системы. В этом случае расчет скорости передачи данных строится на основании произведения (где Р — разрядность системной шины; & — ее тактовая частота).

В реальных же условиях на производительность системной магистрали передачи данных (шина) оказывают влияние всевозможные факторы: проводимость материалов, которая меняется во времени, помехи, недостатки конструкции и сборки, а также многие другие факторы. По данным, взятым из [12,13], разность между теоретической скоростью передачи данных и практической может уменьшаться до 25 % от расчетной величины.

В современных АСУТП, создаваемых на единой программно-аппаратной базе, при возникновении чрезвычайной ситуации возможен конфликт между информационными потоками различных подсистем. Может возникнуть ситуация, при которой первичная информация от высокоскоростных высокопроизводительных источников может быть блокирована менее производительной линией, что приведет к полному отказу системы [14].

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

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

Методика

Для создания модели оценки устойчивости функционирования информационных каналов АСПВБ, основанной на изменении величины информационной энтропии системы [1], первоначально необходимо получить характеристическое уравнение системы

Б

где а0, а1, ..., ап — коэффициенты характеристического уравнения системы;

5 — энтропия системы;

п — количество подсистем АСПВБ.

Известно [4], что необходимым условием устойчивости системы любого порядка является положительность всех коэффициентов характеристического уравнения (1): й0 > 0, й1 > 0, ..., йп >0.

Действительно, уравнение (1) можно представить в виде произведений множителей, содержащих корни 51, 52, ..., Бп:

й„

Если все корни характеристического уравнения будут отрицательны, то все множители выражения (2) будут иметь вид:

йс

где | й11 — значения корней.

Произведя умножение в (3), получим (1), в котором все коэффициенты будут определяться положительными членами | й11 выражения (3), т. е. будут положительными.

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

Для того чтобы оценить устойчивость системы, необходимо использовать алгоритм, основанный на исследованиях английского математика Э. Рауса [4]. Для оценки устойчивости построим таблицу.

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

В данном случае число строк таблицы равно степени характеристического уравнения (1) плюс единица.

ПОЖАРОВЗРЫВОБЕЗОПАСНОСТЬ/ПРЕЛМОЕХРШБЮМБЛРЕТУ 2017 ТОМ 26 №11

15

Матрица коэффициентов характеристического уравнения The coefficient matrix of the characteristic equation

Коэффициент rt Строка Столбец / Column

Coefficient rt Row 1 2 3 4

- 1 a0 = c11 a2 = c21 a4 = c31

- 2 a1 = c12 a3 = c22 a5 = c32

r3 = a0/a1 3 c13 = a2 - r3a3 c23 = a4 - r4C33 c33 = a6 - r3a7

r4 = a1/c13 4 c14 = a3 - r4c23 c24 = a5 - r4c33 c34 a7 r4c43

r5 = c13/c14 5 c15 = c23 - r5c23 c25 = c33 - r5c33 c35 = c43 - r5c43

c1, i- 2 ' = i c1, i = c2, i -2 - ric2, i -1 c2, i = c3, i - 2 - ri c3, i - 1 c3, i = c4, i - 2 - ri c4, i - 1

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

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

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

В этом случае для АСПВБ, построенной по дискретному принципу, величина энтропии 5д определяется по формуле [15]:

Sä = "Ё Pi l0S2 Pi

i=i

где pl — вероятность появления сигнала в i-м канале.

Известно [6], что в АСПВБ всегда задействован хотя бы один информационный канал (например, канал передачи значения давления в технологическом оборудовании). В этом случае вероятность появления сигнала в данном информационном канале равна единице. Если сигналы в других каналах отсутствуют, то, очевидно, вероятность появления сигналов в них равна нулю, т. е. для заранее известного информационного сообщения выражение (4) будет состоять из слагаемых двух видов: либо 1log21, либо 0 log20. Слагаемое первого вида равно нулю, а значит, lim (x  log2 x) = 0 и, следовательно, энтропия за-

x ^ 0

ранее известного сообщения равна нулю.

Следует определить, какое максимальное значение может принимать слагаемое -pl log2pl. Для этого продифференцируем уравнение (4) и приравняем производную нулю:

(d/(dPi)("Pi log2 Pi) =

= "Pi log2 Pi " Pi " Pi (1/Pi )log2 e = 0 т. e. Pi = 1/e.

Так как все слагаемые выражения (4) не превышают значения (1/е) 1о§2 е, то из этого следует, что энтропия АСПВБ есть величина конечная при любом количестве информационных подсистем. Очевидно, при таком распределении вероятностейр1 энтропия принимает максимальное значение, что коррелирует с результатами Э. Рауса.

Теперь воспользуемся методом неопределенных множителей Лагранжа и найдем максимум функции в виде формулы (6), составленной с учетом допол-

п

нительного условия X рг = 1:

г=1

п п

р = Рг1о& Рг - Р1-

г = 1 г = 1

Дифференцируя (6) по pi и приравнивая производную нулю, получаем уравнение (7) или (8):

ХР = - 1оМ2 Рг -1 — | Р11о§2 е - ^ = 0;

" 1оБ2 Р1 = 1о§2 е + ^

Анализируя выражение (6), можно сделать вывод, что вероятность Рг не зависит от переменной суммирования, что может быть только в том случае, если все вероятности равны между собой: р1 = р2 = ... = =Рп = Р = 1/п. Следовательно,

rt max S Д

n i i

= "Ё - l0g2 - = l0g2

i = 1

Известно [3], что энтропия стремится к своему максимальному значению. Основываясь на этом, можно заключить, что энтропия системы достигает максимального значения, когда вероятности появления информационных сигналов равны.

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

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

В случае создания АСПВБ, в которой ряд сигналов зависит друг от друга, величина информационной энтропии определяется по зависимости, отличной от (4).

Тогда систему сигналов (X, У) рассматриваем как один сигнал, а множество состояний такого совокупного сигнала представим состоящими из пар вида (xi, у}), где , = 1, ..., п, } = 1, ..., п. При известных вероятностяхр (х,, у) появления пары (х(, у]) энтропия системы сигналов (X, У) определяется выражением

5

I = 1

Согласно теореме умножения вероятностей имеем:

р

р

где р (х,), р (у]) — вероятности появления сигналов соответственно х, е X, у] е У; р (х, | у) — условная вероятность появления элемента х, е X при условии, что уже имеется сигнал у] е У;

р ( у | х,) — условная вероятность появления сигнала у] е У при условии, что уже имеется сигнал

х, е X.

Подставляя в (9) выражение (10), получаем формулу

5(X, У) = -£ £р(у) рхЬ]) х

=1 ] = 1 (12) х [р(у,0 рС^Ь))].

Учитывая, что 1ОЯ2 (р ( у) р (х,| у)) = 1ОЯ2р (у}) + +1о§2р(х, | у]), последнее выражение запишем в виде:

п п

5(X, У) = -£ р(у) 1ом2 р(у) £ р(х,.^.) -

]= 1 г = 1

-£ р(у.) £ рСх^ у, 1о§2 р( х^о-

]=1 ¡'=1

п

Имея в виду, что £ р(х, \ у.) = 1, получаем выра-

, =1

жение

5(X, У) = -£ р(у/)1оЕ2 р(у.) -

]=1

пп

-£ р(у,0 £ рСх^О ^2 р(щ\у])-]=1 ,=1

второго члена (14) можно рассматривать как энтропию сигнала X при условии, что второй сигнал У получил конкретное значение У,. Эту условную энтропию относительно элемента У, можно назвать частной условной энтропией, которую обозначим как 5 (X | у. ).

Аналогичное выражение можно получить и для 5 (У | х,).

Следует показать, в каких пределах может изменяться энтропия объединения двух информационных сигналов — X и У. Выражение для условной энтропии будет иметь вид:

5(X|У) = -£ £ р(у)р(х1\у]) 1ом2 р(х1\у]). (15) ,=1 ]=1

Известно [9], что для независимых случайных событий условные вероятности равны безусловным. Если предположить, что составляющие X и У сигнала (X, У) независимы, то в выражении (14) возможна замена р (х, | у}) = р (х,), и тогда оно может быть преобразовано к виду:

пп

5

=1

} = 1

В соответствии с (4) первая сумма (14) представляет собой энтропию сигнала У, а вторую сумму

Учитывая, что £ р (у) = 1, имеем выражение

} = 1

5

I = 1

Иными словами, условная энтропия при независимых сигналах равна безусловной энтропии, что приводит к тому, что У) = + 5(У).

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

Пусть сигнал Xпринял значение х1, а сигнал У— значение ук. В этом случае условная вероятность р (ук | х) = 1, а все остальные принимают нулевое значение. Тогда в формулу (15) будут входить слагаемые вида либо р(х,)Ф^21, либо р(х,)0^^20. В любом случае эти слагаемые равны нулю, поэтому 5 (X | У) = 0.

Результаты и их обсуждение

В результате исследований получена формула (14), позволяющая оценивать величину информационной энтропии для оценки надежности функционирования каналов передачи информации.

Тенденцией настоящего времени, связанной с уменьшением финансовых затрат на АСУТП, является создание на единой программно-аппаратной базе автоматизированных систем различного технологического назначения. Однако в данном случае необходимо разбираться в противоречиях, которые

ПОЖАРОВЗРЫВОБЕЗОПАСНОСТЬ/ НРЕЛШ ЕХРШБЮИ БЛРЕТУ 2017 ТОМ 26 №11

17

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

Зависимость (14) должна учитываться проектировщиками АСУТП промышленных предприятий при расчете количества каналов передачи информации в интегрированных системах.

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

Выводы

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

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

Данный подход позволит повысить надежность функционирования АСПВБ, особенно в условиях возможных дестабилизирующих воздействий, при которых объем информации в системе может резко возрастать.

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

1. Бутузов С. Ю. Энтропийный метод определения эффективности функционирования информационных каналов автоматизированных комплексных систем безопасности промышленных предприятий // Вестник Академии Государственной противопожарной службы МЧС России. — 2005.—№4. —C. 6-9.

2. Бородин И. Ф., Кирилин Н. И. Основы автоматики и автоматизации производственных процессов. — М. : Колосс, 1977. — 328 с.

3. Вильсон А. Дж. Энтропийные методы моделирования сложных систем / Пер. с англ. — М.: Наука, 1978.—248 с.

4. Куликовский Л. Ф., Мотов В. В. Теоретические основы информационных процессов. — М. : Высшая школа, 1987. — 248с.

5. Топольский Н. /.Автоматизация систем пожарной безопасности АЭС. — М. : ВИПТШ МВД России, 1994.—200 с.

6. Топольский Н. Г. Основы автоматизированных систем пожаровзрывобезопасности объектов. — М. : МИПБ МВД России, 1997. — 164 с.

7. Ефимов M. В. Теория автоматического управления. — М. : МГУП, 2006. — 420 с.

8. RosenbrokH.H. Computer-aided control system design. — London: Academic Press, 1974.—230р.

9. Вентцель Е. С. Теория вероятностей. — М. : КноРус, 2010. — 664 с.

10. Юркевич Е. В., Азаренкова Н. В., Долкарт В. М.Методология регулирования энтропии в корпоративных информационных системах // Вопросы электромеханики. Труды ВНИИЭМ. — 2013.

— Т. 133, № 2. — С. 39-48.

11. Воронов А. А. Основы теории автоматического управления: Автоматическое регулирование непрерывных линейных систем. — 2-е изд., перераб. — М. : Энергия, 1980. — 312 с.

12. ArocaR. A., Tomkins A., Doi Y., Yamamoto T., Voinigescu S. P. Circuit performance characterization of digital 45-nm CMOS technology for applications around 110 GHz // IEEE Symposium on VLSI Circuits.—2008.—P. 162-163. DOI: 10.1109/VLSIC.2008.4585991.

13. SiljakD. Analyses and syntheses of feedback of control systems in parameters plane //IEEE Journal of Oceanic Engineering. — 2017. — Vol. 42. — P. 253-262.

14. Бутузов С. Ю. Моделирование автоматизированных информационно-управляющих систем интегральной безопасности высокорисковых объектов // Технологии техносферной безопасности.

— 2007. — Вып. 2(12). — 9 с. URL: http://agps-2006.narod.ru/ttb/2007-2/13-02-07.ttb.pdf (дата обращения: 30.08.2017).

15. Шеннон К. Работы по теории информации и кибернетике / Пер. с англ. — М.: Иностранная литература, 1963. — 832 с.

Материал поступил в редакцию 2 сентября 2017 г.

Для цитирования: Бутузов С. Ю., ПрановБ. М., ПрусЮ. В., Семиков В. Л., Яковлев О. В. Модель оценки устойчивости автоматизированных систем пожаровзрывобезопасности // Пожа-ровзрывобезопасность / Fire and Explosion Safety. — 2017. — Т. 26, № 11. — С. 14-20. DOI: 10.18322/PVB.2017.26.11.14-20.

STABILITY EVALUATION MODEL OF COMPUTER-AIDED FIRE-EXPLOSION SAFETY SYSTEMS

BUTUZOV S. Yu., Doctor of Technical Sciences, Associate Professor, Professor of Department of Information Technology, State Fire Academy of Emercom of Russia

PRANOV B. M., Doctor of Technical Sciences, Professor, Professor of Department of Applied Information Technology, Russian Presidential Academy of National Economy and Public Administration

PRUS Yu. V., Doctor of Physical and Mathematical Sciences, Professor, Head of Department of Management and Economics of State Fire Service, State Fire Academy of Emercom of Russia

SEMIKOV V. L., Doctor of Technical Sciences, Professor, Professor of Department of Management and Economics of State Fire Service, State Fire Academy of Emercom of Russia

YAKOVLEV O. V., Doctor of Technical Sciences, Associate Professor, Leading Researcher, Federal Research Centre "Information and Management"

English

ABSTRACT

In the conditions of the modern economic growth, the urgency of the reliable functioning of the computer-aided system of industrial enterprises, which are critical for the Russian economy, is growing.

The computer-aided fire-explosion safety systems now play a key role in the structure of computer-aided system. At the same time, they are a superfast system that can detect overpressure of explosive process units within 0.1 microseconds. In this case, it is possible to carry out a real-time pressure relief to protect the technological equipment from destruction and the personnel of the enterprise from destruction.

Computer-aided fire-explosion safety systems are part of computer-aided systems of technological processes. Value of operating frequency of system bus of these systems is approximately 100 MHz. So, sustainability assessment function of information transfer channel is important scientific and technical problem.

There is mathematical model of assessment of the stability of information transfer channels in the article. This model is based on scientific, theoretical and mathematical models of E. Rouse and on entropy methods of A. Wilson.

Simulation was carried out of reliability of information transfer channels for various options for building these systems. It is shown that an integrated approach based on the use of a single hardware/ software platform of a microprocessor platform can be used to optimize the single of creating an ASP.

Potential destabilizing effects have a negative impact on the sustainability of the work of ASIA. Under these conditions, the amount of information sharply increases in the system. And the model of estimating the entropy will allow the reliability of the automated system as a whole to depend.

Keywords: computer-aided process control systems; system data line; computer-aided fire-explosion safety systems; stability and reliability of technical systems; entropy of technical system; entropy evaluation model of sustainability.

ПOЖAPOBЗPЫBOБEЗOПACHOCTЬ/FIRE AND EXPLOSION SAFETY 2017 TOM 26 №11

1S

REFERENSES

1. Butuzov S. Yu. An entropy method for determining the efficiency of the functioning of information channels of automated integrated security systems for industrial enterprises. Vestnik Akademii Gosu-darstvennoy protivopozharnoy sluzhby MChS Rossii / Bulletin of State Fire Academy of Emercom of Russia, 2005, no. 4, pp. 6-9 (in Russian).

2. Borodin I. F., Kirilin N. I. Osnovy avtomatiki i avtomatizatsiiproizvodstvennykhprotsessov [Fundamentals of automation and automation of production processes]. Moscow, Koloss Publ., 1977. 328 p. (in Russian).

3. Wilson A. G. Entropy in urban and regional modelling. London, Pion Limited, 1970 (Russ. ed.: Wilson A. G. Entropiynyye metody modelirovaniya slozhnykh sistem. Moscow, Nauka Publ., 1978. 248 p.).

4. Kulikovskiy L. F., Motov V. V. Teoreticheskiye osnovy informatsionnykhprotsessov [Theoreticalbasis of information processes]. Moscow, Vysshaya shkola Publ., 1987. 248 p. (in Russian).

5. Topolskiy N. G. Avtomatizatsiya sistem pozharnoy bezopasnosti AES [Automation of fire safety systems of nuclear power plants]. Moscow, High Engineering Fire Technical School of Ministry ofthe Internal Affairs of Russia Publ., 1994. 200 p. (in Russian).

6. Topolskiy N. G. Osnovy avtomatizirovannykh sistem pozharovzryvobezopasnosti obyektov [Basics of automated fire and explosion safety systems]. Moscow, Fire Safety Institute ofthe Ministry of Internal Affairs of Russia Publ., 1997. 164 p. (in Russian).

7. Efimov M. V. Teoriya avtomaticheskogo upravleniya [Theory of automatic control]. Moscow, MGUP Publ., 2006. 420 p. (in Russian).

8. Rosenbrok H. H. Computer-aided control system design. London, Academic Press, 1974. 230 p.

9. VenttselE. S. Teoriya veroyatnostey [Probability theory]. Moscow, KnoRusPubl., 2010. 664 p. (in Russian).

10. YurkevichE. V., AzarenkovaN. V., Dolkart V. M. Methodology of entropy regulation in corporate information systems. Voprosy elektromekhaniki. Trudy VNIIEM/ Electromechanical Matters. VNIIEM Studies, 2013, vol. 133, no. 2, pp. 39-48 (in Russian).

11. Voronov A. A. Osnovy teorii avtomaticheskogo upravleniya: Avtomaticheskoye regulirovaniye nepre-ryvnykh lineynykh sistem [Fundamentals ofthe theory of automatic control: Automatic control of continuous linear systems]. Moscow, Energiya Publ., 1980. 312 p. (in Russian).

12. ArocaR. A., Tomkins A., Doi Y., Yamamoto T., Voinigescu S. P. Circuit performance characterization of digital 45-nm CMOS technology for applications around 110 GHz. In: IEEE Symposium on VLSI Circuits, 2008, pp. 162-163. DOI: 10.1109/VLSIC.2008.4585991.

13. Siljak D. Analyses and syntheses of feedback of control systems in parameters plane. IEEE Journal of Oceanic Engineering, 2017, vol. 42, pp. 253-262.

14. Butuzov S. Yu. Modeling of automated information-control systems ofintegral safety ofhigh-risk objects. Tekhnologii tekhnosfernoy bezopasnosti / Technology ofTechnosphere Safety, 2007, issue 2(12). 9 p. (in Russian). Available at: http://agps-2006.narod.ru/ttb/2007-2/13-02-07.ttb.pdf (Accessed 30 August 2017).

15. Shannon C. Raboty po teorii informatsii i kibernetike [Works on information theory and cybernetics]. Moscow, Inostrannaya literatura Publ., 1963. 832 p. (in Russian).

For citation: Butuzov S. Yu., Pranov B. M., Prus Yu. V., Semikov V. L., Yakovlev O. V. Stability evaluation model of computer-aided fire-explosion safety systems. Pozharovzryvobezopasnost / Fire and

Explosion Safety, 2017, vol. 26, no. 11, pp. 14-20 (in Russian). DOI: 10.18322/PVB.2017.26.11.14-20.