В настоящее время при разработке устройств автоматической пожарной сигнализации все большее внимание уделяется необходимости применения многофакторного детектирования пожара. Извещателей, позволяющих контролировать несколько факторов пожара (повышение температуры, оптической плотности среды, концентрации угарного и углекислого газа), принято называть многокритериальными. Использование нескольких каналов обнаружения в многокритериальных извещателях позволяет снизить вероятность ложного срабатывания системы, а также сократить время детектирования пожара. К ложным срабатываниям извещателей могут приводить такие воздействия, как запыление дымовой камеры, выпадение тумана, наличие технологических аэрозолей, резкое повышение температуры, связанное с особенностями технологического процесса, и т.д. Действие этих факторов на чувствительные элементы пожарных извещателей сходно с действием факторов пожара, однако, поскольку процесс горения характеризуется совокупным проявлением нескольких факторов, многокритериальные извещатели имеют большую возможность его детектирования.
Для надежного детектирования пожара информация, поступающая от каждого канала многокритериальных пожарных извещателей, должна быть соответствующим образом обработана. Большинство современных средств обнаружения пожара производят обработку информации о различных факторах пожара внутри извещателя. Такой подход приводит к существенному удорожанию самого извещателя, но позволяет сократить время, требуемое на обработку информации, поскольку большая часть вычислений производится распределенными средствами. При этом за основу алгоритма все же принимают информацию одного из каналов обнаружения (чаще всего уровень задымления), а информация с остальных каналов лишь корректирует процесс принятия решения. В данной работе, напротив, предлагается использовать многокритериальные извещатели без предварительной обработки информации. Их задача сводится лишь к измерению контролируемых параметров среды и передаче их значений на приемно-контрольный прибор. Обработка информации со всех каналов всех извещателей будет производиться непосредственно в самом приемно-контрольном приборе. При такой схеме можно оперативно изменять настройки алгоритма, адаптируя всю систему к условиям конкретного защищаемого объекта. Основная роль в этом случае отводится алгоритмам обработки информации и детектирования пожара внутри прибора. Целью данной работы является создание такого алгоритма обработки информации от многокритериальных извещателей, позволяющего снизить количество ложных срабатываний систем автоматической противопожарной защиты.
По своей сути любая система обнаружения пожара является системой классификации. На входы системы подается информация о параметрах среды от автоматических пожарных извещателей и система должна отнести данную конкретную ситуацию состояния измеряемой среды к классам «пожар» или «норма». В данной работе предлагается создать экспертную систему классификации, используя вероятностный подход, а конкретнее - вероятностный классификатор на основе теоремы Байеса. Один из примеров применения вероятностного подхода для решения задачи детектирования пожара описывается в [3].
Для простоты будем рассматривать два класса состояния контролируемой среды и полагать, что вероятность пожара (Fire) P(F) и вероятность нормального состояния среды (Normal) P(N) образуют полную группу событий P(F) + P(N) = 1.
Многокритериальный извещатель через равные промежутки времени передает на приемно-контрольный прибор вектор состояния контролируемых параметров среды . Так, например, w1 - температура среды, w2 - оптическая плотность среды и т.д. Для детектирования пожара в этих условиях необходимо принять решение: принадлежит ли данное состояние среды к классу «пожар»? Для ответа на поставленный вопрос в контексте вероятностного подхода нам необходимо вычислить условную вероятность
, т.е. вычислить, какова вероятность того, что в контролируемой зоне возникнет пожар при данных значениях измеренных факторов.
Используя теорему Байеса, можем получить следующее выражение:
(1)
где - вероятность пожара при условии состояния среды w‾; P(F) - вероятность пожара;
- вероятность состояния среды w‾ при условии пожара;
- вероятность состояния среды w‾.
Аналогично можно оценить вероятность отсутствия пожара
(2)
Для удобства классификации чаще всего рассматривают отношение правдоподобия и сравнивают его с порогом классификации h:
(3)
Логарифмируя выражение (3) получим т.н. характеристическую функцию:
(4)
где θ() - функция Хевисайда.
Получив распределение плотности вероятностей ,
и зная отношение априорных вероятностей
, можно определить оптимальное правило классификации. В данной формуле и может быть вычислено на основе известных статистических данных о частоте возникновения пожара в зданиях рассматриваемого типа, используемых при расчете пожарных рисков [2]. Поскольку заранее о характере закона распределения нам ничего не известно, будем считать, что параметры вектора
независимы. Кроме того, если использовать непрерывные величины, нам придется численно интегрировать выражения с интегралами высокой кратности, что проблематично, поэтому в дальнейшем будем оперировать дискретными переменными.
Говоря о восстановлении функции распределения вероятностей в контексте статьи, мы подразумеваем нахождение
где
, при
(5)
где pi(k) - вероятность того, что xi примет значение ci(k).
Для восстановления функции распределения вероятности можно еще раз прибегнуть к формуле Байеса. В [1] показано, что в этом случае для малого объема выборки
(6)
где mi(k) - число векторов выборки, у которых координата принимает значение xi = ci(k), а l - объем выборки, так что
Оптимальным решающим правилом будет:
(7)
где - n-мерный бинарный вектор размерности
полученный из вектора последовательной записью
для всех wk; - вероятность того что xk = 1 для класса «Пожар»;
- вероятность того что xk = 1 для класса «Норма».
Правило (7) при H = 0 оптимально для равноценных ошибок классификации первого и второго рода [4]. В нашем случае ошибки второго рода неприемлемы, т.е. экспертная система должна работать таким образом, чтобы полностью исключить возможность принятия пожара за нормальное состояние. В таких условиях при наличии двух непересекающихся, линейно разделимых классов, Байесовская классификация сводится к поиску значения порога классификации H, доставляющего минимум функции потерь.
Количественную оценку mi(k) для каждого класса можно получить в ходе обучения классификатора, т.е. оценить количество срабатывания канала извещателя wi при условии нахождения в помещении с пожаром или без пожара (при воздействии ложных факторов) соответственно. С этой целью можно провести огневые испытания или использовать математическую модель горения. Поскольку при развитии пожара срабатывание извещателя произойдет в любом случае (вопрос лишь во времени срабатывания от начала горения), то для практической оценки частоты срабатывания канала извещателя необходимо рассматривать ограниченный промежуток времени, сопоставимый с предельным временем реакции извещателя, установленного нормативными документами.
Поскольку горение различных веществ характеризуется различной интенсивностью проявления того или иного фактора пожара, обучение классификатора необходимо проводить на моделях как можно более широкого типа. В данной работе было проведено моделирование шести типовых очагов горения, используемых для огневых испытаний по ГОСТ Р 53325-2009 в программной среде Fire Dynamics Simulator [5]. Примеры полученных зависимостей изображены на рис. 1, 2. Помимо моделей динамики факторов пожара типовых очагов горения, для цели обучения классификатора использовались модели ложных факторов пожара. В работе использовались следующие модели ложных факторов пожара: запыление дымовой камеры, технологические испарения, включение мощного нагревательного прибора (промышленная тепловая пушка), попадание в дымовую камеру мелкого насекомого и грозовой разряд.
Для практической оценки эффективности предложенного алгоритма использовались модели комбинированного извещателя с дымовым и тепловым каналом. С целью повышения характеристик экспертной системы по каждому каналу обнаружения выделялись максимальная и дифференциальная составляющая. Таким образом, для каждого комбинированного извещателя рассчитывались четыре параметра: текущая температура (T, °С), скорость изменения температуры (ΔT, °С/с), текущая оптическая плотность среды (μ, дБ/м), скорость изменения оптической плотности среды (Δμ, дБ/(м∙с)).
Рис. 1. Динамика опасных факторов пожара (температура) для тестовых очагов горения ТП-5 и ТП-6
Рис. 2. Динамика опасных факторов пожара (оптическая плотность среды) для тестовых очагов горения ТП-5 и ТП-6
Как было указано выше, для упрощения процедуры вычислений необходимо перейти от непрерывных к дискретным входным величинам. Переход к дискретным значениям величин можно выполнить следующим образом. Нормативными документами в области пожарной безопасности установлено, что оптический пожарный извещатель должен срабатывать при оптической плотности среды μ = 0,05...0,2 дБ/м. Очевидно, что при оптической плотности среды μ < 0,05 дБ/м состояние извещателя можно назвать «норма», а при оптической плотности среды более μ > 0,2 дБ/м - «пожар». Пограничное состояние при оптической плотности среды 0,05 < μ < 0,2 дБ/м условно назовем «внимание», поскольку нельзя однозначно отнести его ни к состоянию «пожар», ни к состоянию «норма». Таким образом, произведена дискретизация входной величины по уровню. Аналогичный подход применяется и к оставшимся трем параметрам.
Статистические оценки и
по данным математических моделей определяются как частота наблюдения события за ограниченный промежуток времени. Стоит отметить, что обучающий набор данных должен быть достаточно большим и придерживаться соотношения
, т.е. количество обучающих данных для моделей «пожар» и «нормальное состояние» должны соответствовать отношению априорных вероятностей.
Тестирование обученной системы производилось на образах, зарезервированных в генеральной выборке и не предъявляемой системе во время обучения. При тестировании экспертной системы математическое ожидание ошибки первого рода (система приняла ложный фактор пожара за настоящий пожар) составило ≈36%, в то время как при использовании традиционного алгоритма (логическое «ИЛИ») этот показатель находится на уровне ≈66%. Среднее время детектирования тестовых очагов пожара экспертной системой не превышало аналогичного показателя традиционного алгоритма «ИЛИ».
К достоинствам предлагаемой экспертной системы следует отнести существенное снижение количества ложных срабатываний по сравнению с традиционными алгоритмами детектирования пожара. К сожалению, не лишена данная система и недостатков. Полагая опасные факторы пожара независимыми друг от друга характеристиками среды и проводя обучение на ограниченной выборке ситуаций (моделей), Байесовская система классификации теряет свойство оптимальности. В общем случае надежная классификация достигается только в области ситуаций, предъявляемых системе в ходе обучения. Кроме того, алгоритм предназначен для обработки информации поступающей с нескольких каналов одного извещателя. Для получения алгоритма, учитывающего информацию нескольких мультикритериальных извещателей, требуется дополнительно использовать корреляционные методы. В этом случае требуемое быстродействие процессора обработки информации растет экспоненциально.
Рецензенты:
-
Добротин С.А., д.т.н., директор ООО «НТЦ «Безопасность», г. Дзержинск;
-
Луконин В.П., д.т.н., генеральный директор Федерального государственного унитарного предприятия «Научно-исследовательский институт химии и технологии полимеров имени академика В.А. Каргина с опытным заводом», г. Дзержинск.
Работа поступила в редакцию 19.03.2012